Séparabilité simple de l'équation de Helmholtz ou de l'équation de Laplace pour des systèmes de 
coordonnées de révolution dans l'hyper-espace 


On dénomme x, la hauteur et p le rayon de révolution. Ces deux variables du système de révolution 
sont définies par l'intermédiaire d'une fonction analytique du plan complexe, que l'on va 
dénommée fonction de révolution. ( Voir 1931, EWHobson « The theory of spherical and 
ellipsoidal harmonics » chapitre 10 « Harmonics for spaces bounded by surfaces of revolution », 
1959, L.Robin e Fonctions sphériques de Legendre et Fonctions Spheroidales » Tome III, chapitre 9, 
« Applications des fonctions de Legendre aux surfaces de revolution, avec des systèmes de 
coordonnées curvilignes orthogonales autre que les coordonnées sphériques , 1939, R.Lagrange 
« Les familles de surfaces de révolution qui possèdent des harmoniques », Acta Mathematica 71, 
page 283, repris en note de bas de page dans l'ouvrage de L.Robin , chapitre 9 page 180) 


Le plan complexe correspond à l'un quelconque des hyper-plans où sont définies les 
e méridiennes » des iso-variétés du système de coordonnées orthogonales : Zu +ip=f(é+i n), La 
fonction de révolution possède la propriété « miroir » : KEE Si L'analyticité de la 
fonction f entraîne que toutes les iso-surfaces de révolution n=Cste et Ÿ=Cste sont orthogonales 
l'une de l'autre puisque : 
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Ce qui est la condition d'orthogonalité des iso-courbes sur l'hyper-plan méridien (les vecteurs 
normaux des deux iso-variétés sont perpendiculaires). Les deux nouvelles variables €n et les angles 
restant des rotations de SO(N-1) forment un système de coordonnées orthogonales car elles sont 
définies à l'aide de fonctions analytiques a priori quelconque du plan complexe. L'ensemble des 
rotations de SO(N-1) autour de l'axe x0 des courbes engendrée £&=Cste ou n=Cste forment des iso- 
variétés de dimension N. On a par ailleurs : 
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Tout se passe comme si l'on écrivait les formules de passage vers les coordonnées cartésiennes 
suivantes : 


x, = plé,n)Sin(s, , )Sin(9, - ).Sin(9, )Sin(9) 9 =pel0,2r] 
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En introduisant les rayons successifs : p} = Se et spécialement p = D = p° 
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Nous pouvons exprimer ces rayons en fonction du système d'angle : 
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Comme la métrique de ce système de coordonnées s'écrit : 


Su =dp?° + D d3? 


TA F(E tinŅdé tian) dp —+idx, = f'(£-+in\dé -+idn) 
> dp +dx = f'(E+in)f'(E+in\dE? +an°) 
> dx? + D ds = f'(E+in)f'(E+in dE + dn? be Se (TECH 


Hi 
Le laplacien a la forme suivante dans cette métrique : 
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Ce qui donne finalement : 
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On retrouve l'expression déjà établi dans le cas à 3 dimensions, comme N=2 : 
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L'équation de Helmholtz a également la forme suivante : 
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En introduisant comme auparavant l'opérateur différentielle du moment angulaire sur l'hyper- 
sphère N-1, il vient pour les deux équations aux dérivées partielles : 
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Remarque importante : si la fonction de révolution échange les coordonnées p et x, alors le produit 
suivant reste le même JL +in)f'(€-in), C'est le cas lorsque l'on opère une transformation de la 
fonction de révolution simple de multiplication par i (ou de rotation dans le plan complexe de pet 
Xo). Cette opération est utile pour trouver des systèmes de coordonnées orthogonales 
complémentaires par rotation autour de l'axe p et non de l'axe xo. 


Déterminons des conditions particulières pour lesquelles l'équation de Helmholtz devient séparable 


et par suite l'équation de Laplace l'est pour w=0. Supposons que la solution T soit séparable en 
trois fonctions des coordonnées curvilignes : 
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Le terme angulaire doit être constant dans la séparation, il vient donc : 
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L'équation de Helmholtz devient donc : 
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L'équation de Laplace quant à elle devient : 
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Résumons la recette pour trouver les équations simplement séparées des équations de Laplace 
et Helmholtz d'après la fonction de révolution : 
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Supposons que 


Exemple de simple séparabilité : système hyper-cylindrique 


fE+in)=Éé+in=x+io f(E+in)=/f'(E-in)=1 

PRÉPA) =). ES E Some 
ne a Sg An? zlé) | Së 
J'E+in in S ;J'Etin)+ fin) | Se 
EE Po) 0 


G(E)=G(n)=1 


Il 
© 


LD Joie: ulu +22 -1)Q(0) 


Exemple de simple séparabilité : système hyper-sphérique 
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Exemple de simple séparabilité : système hyper-sphéroïdal allongée 


E f'E-in)=cSinh(é-in) 
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Pour obtenir les équations différentielles usuelles du système hyper-sphéroïdal allongé, il suffit 
d'effectuer les changements de variable suivant : 
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Exemple de simple séparabilité : système hyper-sphéroïdal aplati 
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Pour obtenir les équations différentielles usuelles du système hyper-sphéroïdal aplati, il suffit 
d'effectuer les changements de variable suivant : 


siet Z SEL cosie EE) Tane E (E)E =E) 


E = Sinh) | PE A AEN one. Se, Es os 
a Së, Gell zl El-EzE EEE) 
H'(n)=-Sin(n)H' (Ñ) —> Cotan(n)H'(n)= -7H (i7) 
H" (n)= 4 ()+(-5°)4"(G) 


Fees nes (hf zech A ae 


E?+1 
fm) A Casa mt a -oeh Heil 
Eeler dëse: ëch ër eA aeo 
Gr) arr er ra ae bal HUE D) #n)=0 


Exemple de simple séparabilité: système hyper-parabolique de révolution 
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Nous allons maintenant passer à l'étude de la R-Séparabilité pour un système de coordonnées de 
révolution dans l'hyper-espace. 


R-Séparabilité de l'équation de Laplace pour des systèmes de coordonnées de révolution dans 
l'hyper-espace 


Supposons maintenant que la solution soit R-Séparable de la forme : 
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Le terme angulaire doit être constant dans la séparation: Li oe, H (u ES -10(6)- 0 
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On retrouve l'expression lorsque la fonction R-Séparable est constante égale à 1 : 
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Revenons au cas général d'une fonction de R-Séparation non constante. La première condition pour 


laquelle l'équation de Laplace peut être R-séparable porte sur la séparation des deux termes 
suivants : 
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Cette condition implique donc des relations entre les grandeurs introduites, que l'on va exploiter. 


Comme Ap>0 et SEA 


Pour la suite des calculs sur la séparation potentielle de l'équation de Laplace, intéressons-nous 


maintenant aux deux termes : _Ÿ oi? CRL al ô sf? GL al , il vient : 
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Si l'on revient à une fonction de séparation constante, soit le cas de la simple séparabilité, on sait 
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L'expression : HIEL, (EE) i AR Lal e WIR. redonne bien une valeur nulle : 
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Une autre vérification consiste à se reporter au cas à 3 dimensions d'où N=2. 
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devient FE), (EE) E'n), Lal , AP que l'on avait déjà calculée dans un autre document 
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sur les système de coordonnées orthogonales de révolution. 


Continuons le calcul entrepris pour les deux termes : Q D Toi ô D alem) afin 
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Les équations séparées issues de l'équation de Laplace deviennent alors : 
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On retrouve bien les équations différentielles séparées lorsque la fonction de R-Séparation est égal 
à 1 : 
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Résumons la recette pour trouver les équations R-séparées de l'équation de Laplace d'après la 
fonction de révolution et la fonction de séparation : 
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Lorsque la fonction de séparation est définie comme l'inverse de la racine carrée du rayon p à la 
puissance N-1, alors la forme des équations R-Séparée s'en trouve simplifiée : 
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Chers | 


Lorsque la fonction de séparation est une constante, soit 1, alors les fonctions F1 et F2 sont 
déterminées complètement par la fonction de révolution, on est alors dans le cas de séparabilité 
simple de l'équation de Laplace : 


d f'G+in)-f'E-in) _ z 
MOA me Em) OO peinte 
E(n)=i(n Peu A DË (n) (Fé +in)-flé-in)) 
EE 

H"(n)+ Ein) H'(n)+(4ulu +24 -1)z;(n)+a)H(n)=0 


=x(é)+ x(n) 


Remarque : plusieurs fonctions de séparation sont possibles pour rendre R-séparable l'équation 


de Laplace à N+1 dimensions. La plus belle illustration en est dans le système hyper-toroïdal où 
deux fonctions de séparation peuvent être utilisées, dont l'une d'elle est effectivement l'inverse de 
la racine carrée du rayon de révolution à la puissance N-1. De plus dans tous les exemples suivant 
N-1 sera à chaque remplacé par 24. 


Exemple de R-séparabilité : système hyper-toroidal 


Ce sont essentiellement les mêmes calculs que ceux déjà effectués pour le systèmes toroïdal en 3 
dimensions. La fonction de révolution s'écrit comme suit : 


x +ip = f(+in)= cCoran| ` i es j teotan ET) 
IE aA Ein) E E 
H 


Sin(n) 
SE Coran © + A) Sin(n)+iSinh(£) ER Cosh(£)- Cos(n) 
2 " Cosh(e)- Cos(n) Ss Sinh(£) 
p=c 
Cosh(£)-Cos(n) 


Les éléments de calcul pour la R-Séparabilité sont les suivants : 


=> x +ip = x -ip = f(E+in)= f(E-in)= p = 


Xo =C 


b d É+in d c 
x, +ip = flé +in)=ic SÉ IL i in)= fe in)= - 
2 2Sinh E 2sm $5) 


E ER e E Sinh(£) SS wë Sen 5 Sinh’ (é) 
SENTE e geht" VE A AN 
f'(n+i0)f'(n-i0) ` he 8 1 S 
SE = eet" DEN hl al eut et 2(n)=0 


De plus la fonction de R-Séparation est choisie comme suit : 


N- R? — ANERE) F, (n) X Sinh(£) 
d Ge ` SE © Cosh(E)- Cola! 


= RE (En) = (CoshlE)- Cost)" = Rha, 0)= (CoshlE)- Con) ei lala? seng 
> ME 21> F (n)=0 e55 ® — eil) Sinh (n) 


> faé F(E)=(N -1)Log (Sinh(é))+ (N —1)Log(c)= F(E)=(N -1)Cotanh(é)=> F'(£)=- 


N-1 
Sinh? (£) 


On obtient alors les équations séparées suivantes : 


eeh eidelt sat als: EE EO bad 


ao) du o 0 sai aiser: EN, EI, bai 


Smtëel34 SE 


H"(n)+a H(n)=0 


D Le 5 (u+24-1)\) Aa 
R"(€)+ 22 Cette: RE bs 0 


4A(1-2)-4u(u+24-1) 28  ,4%Cotank? (E) _ Stein 
4Sinh°(E) 2Sinh°(£) À Den 


H"(n)+a H(n)=0 


dëng A e 2 = u(u+22-1) 8 
En posant Ge OG =v(v+22) = F CARRE fa- Sinh’ (£) k Ga 
H'(n)+(& + 2 }H(n)= 0 


Solution  T(E,n,@)=(Cosh(£)- Cos(n)) > 2 - a)4(m)Q(G) 


ä+2 =(v+1) 


Avec le choix de la fonction de séparation comme inverse de la racine carrée du rayon de révolution 


à la puissance N-1, on obtient d'autres équations séparées, il vient la R-Séparabilité-2 : 


Cosh(£)- Seel 
Sinh(£) 


1 


HOT 


ReneP S ago) 


e, A(-2)-u(u+22-1) le 8 E 
> éi l Sinh’ (£) lze) j Solution ren o)-( CAE) Er) 
H"(n)+a H(n)=0 
De plus  u(u+24-1)-2(1-2)=(u+A2Xu+2-1) 
: mep |-( D ai 
H"(n)+a H(n)=0 


Exemple de R-séparabilité : système « Bi-hyper-sphérique » 


— 
em 
le 


KS x:(n)=0 Flël- Eloi 


Dë 
| 2 


LG io ae) 


Dans ce cas les calculs de R-séparabilité donnent d'après la forme de la fonction révolution donne : 


DE sa ip = f(E ternet! SH > Xp +IP = Xo ip = f(&-in)= flé +in) 


= Sinh(£) 
Send Co? - gd Sinh(£)+iSin(n) E ° Cosh(£)-Cos(n) 


a Cosh(£)-Cos(n) =. Sin(n) 
Cosh(£)-Cos(n) 
: à -in 3 . Y , Y 
sa tip = f(E -in)=y Coran ) f'E-in)= — J'(E+in) 
2 asim ( 55) SA 


2 2 


S'E+in)f(E-in)- À —= ` > 
ST EA DR sac (Cosh(£)- Cos(n)) 


Sin’ (n) 
(Cosh(£)-Cos(n)) 


=4E)+%(n)= 1(1)=0 et x(n)= 


(FE +in)-f(E-in) =40° =-47 


f'E+in)f(E-in) ` 1 
(FlE+in)-flE-in)}}  4Sin?(n) 


De plus la fonction de R-Séparation est choisie comme suit : 


1 
| 4Sin” (n) 


N-I R2  LIdéA(é) jdn P(n) 8 Sin(n) 
E Ga S g SZ € Cosh(E)- Cosi) 


= R? (E.n) = (Cosh(E)- rela = R(En)= (CosH(E)-Cos(n) > ei "ell ege" Val 


ell A(6)= Eil 


ll UN —1)\(Log (Sin(n))+ Log(c)= F, (n)=(N -1)Cotan(n)= F,'(n)= SE 


On obtient alors les équations séparées suivantes : 


=) RODO Gulu 21-)-440-4) E)E €), Se ca 


2 4 
F'n), "Na: 0 


H"(n)+ P(o) ei) (sata +24-1)-440- Asbl 4 


S Siet none (rh A-Z 27 E HL 


viel 22 Coran) + la z d 


zrtel-Je-27lnutol-0 


H" (n)+ 24 Cotan(n)H (n)+ fa Gel -0 


En posantä =a-—X communément & =v(v +24) => 


Solution  T(E,n,@)=(Cosh(£)- Cosn)) + SLë in Kiel 


Avec le choix de la fonction de séparation comme inverse de la racine carrée du rayon de révolution 
à la puissance N-1, on obtient d'autres équations séparées, il vient la R-Séparabilité-2 : 


DL als m= RE) (EC) nE)=0 à zl ar F(E)=F(n)-0 


Fe Sin(n) | ASin(n) 
=" (¢)-a E(é)=0 xa 


Exemple de R-séparabilité : coordonnées inverses hyper-sphéroïdales allongées 


Dans ce cas la fonction a la forme : 


TRE y : Cos(n)Cosh(E)+ i Sin(n)Sinh(£) E LE) e ; 
WRS coe * Cosh’ (§)- Sin’ (n) E 
Cosh(E)Cos(n) _. Sinh(£)Sin(n) 


o =Y Cosm (E)- Sira) P” Cosh (E)- Sina) 
Les calculs de R-Séparabilité donnent : 


Sinh(£ +in) 
Cosh’ (é +inņ) 


Sinh(é — in) 
Cosh’ (£ —in) 
2 Cosh’ (&)- Cos’ (n) 
a | E 7 le o1 EE 
i V > Sink’ (é)Sin’ (n f(E+in)-f(E-in)) 4 Sink’ (é)Sin’(n 
Lrio +i0)- f(n -i0)? =-4y EE 


f(£-in)= f'(n-i0)=-y f'(n+i0)=-y 


d 
Cosh(£ —in) 


f'(E+in)f(E-in)=7 


1 1 
AOC) af) 


De plus la fonction de R-Séparation est choisie comme suit : 


N- R2 — LIdER(E) ldnF(n) 5 Sinh(£)Sin(n) 
M g PT Cosh? (E)- Sin (1) 


= R?(E,n)= (Cosh? (E)- Cos?) = RE, n) = (Cosh? (E)- Cos? (m)? => PIRE, ale Sinh” (E)Sin”™>(n) 


= fdg FA (E)=(N— Log Larei Logle) = Giel -Cornil Del: ` set 
= jdn Fale le (Sinln))+ Log(e))= F61) = (N = 1)Coran(n)= F'O) = ta 


On obtient alors les équations séparées suivantes : 


= (nee |lu) -aeh DI DI aae) 
mer nr) naar) DI EO a la- 
E"(E)+22 Cotanh(£ E'(E)+1 2 Her erer 
> HAS En posant ä=a-X communément o =v(v +22) 
rela 2A Cotan) E (+ | 2 ME an) 0 
E" (E)+2A Cotanh(é E'(E)- 1a +” u gi z(é)=0 D 
> Solution TS. gl Woh"). Cos?’ (n) : (EH melO) 
la 24 Cotan()H'(n)+ g- PEA Hat, 0 


Avec le choix de la fonction de séparation comme inverse de la racine carrée du rayon de révolution 
à la puissance N-1, on obtient d'autres équations séparées, il vient la R-Séparabilité-2 : 


rene sen DEST aen 3 ail ei FO eil? 
Si (u+AXu+2-1) | 
Stëll Ell A ee 

| SARE) to reno) LEO) SEH) 


Exemple de R-séparabilité :coordonnées inverses hyper-sphéroïdales aplaties 


Dans ce cas la fonction de révolution a la forme : 


Fr c 8 Cos(n)Sinh(£)+iCosh(£ )Sin(n) _, FAN E NS ; 
nepen n E E X +ip =x, -ip = f(¢-in)= f(E +in) 
des Cos(n)Sinh(£) 27 Cosh(é )Sin(n) 


Cosh?’ (£)- Cos? (n) d Cosh?’ (£)- Cos? (n) 


Les calculs de R-Séparabilité donnent : 


f(£-in)= f'É-in)=-c 


c Cosh(£—in) n E Cosh(E +in) 

Sinh(£ —in) Sinh?’ (E —in) Ein) ” Sinn (E +in) 

E + 2. Cosh°(£)-Sin°(n) 

Ga à ele, eener 
AE E E Sin?’ (n)Coshk’ (£) f(E+in)-f(E-in) 4 Cosh’ (é)Sin’(n 

Limal Gate Con} 


> Së 4(n)=- 


1 
4Cosh°(£) 


De plus la fonction de R-Séparation est choisie comme suit : 


1 
4Sin°(n) 


eelerer! a es Det R(E,n)= (Cosh? (8) Cost) F 

=> DURE ,n)=c""Cosh""(E)Sin" "(n) 

Se e F(E)=(N—1XL0g (Coshlé))+ Log(c) = F(6)=(N —1)rann(é) 
[an F,(n)=(N -1)Log (Sin(n)) = F, (n)= (N -1)Cotan(n) 


=> F'(#)= F,'(n)= 


Co Sa 


On obtient alors les équations séparées suivantes : 


= (greg ader drai DE, EI bg 


Re | 22 42 EN, EN 


eil Tam Sc, e Ia: 
> FAG G=a-X communément & =v(v +24) 
elo 24 Con) a HET D) 0 
In 


SE da? 


Cosh (é) 


N-1 


en 


> Solution T(é,n,®)= ENER Cos?(n)) ? Siet ie) 


tal 22 Cotan) H'(a)+ |a- d D 0 


Avec le choix de la fonction de séparation comme inverse de la racine carrée du rayon de révolution 
à la puissance N-1, on obtient d'autres équations séparées, il vient la R-Séparabilité-2 : 


sie EE Set obt 


l 
Cosh(£ )Sin(n) 4Cosh°(£) S 


1 
Aale De 


ut GE "le (E)=0 Ss 
Cosh? olution A Cosh’ (é)- Cos’ (n) Ke 
pe D, alao soon remo) [ET ere) 
Sin” ( 


Exemple de R-séparabilité : le système Hyper-cardioïde de révolution ou e hyper-parabolique de 
révolution inversé » 


Dans ce cas la fonction a la forme : 
ne c 8 c e BIER -n° )+ 2ién) ae c((e° - oh 2ién) 
HG n) AE -inf JEFTA) 2e +7) => f- in) de Hr 
` E?-7°) 2 én 
2e +7) e "Leef 
Les calculs des fonctions de séparation donnent : 


Eeer ek IN E 


Xo 


En mi PESE iner 


Grimm sani 
zij PUS nr mur Erm Ein) | (EE) 
P Tr] E MG Ou C erry Etim- SE- iny En? 


El 0 


De plus la fonction de R-Séparation est choisie comme suit : 


PR (En) OMR) a poef si REE 
Goal 

3 PYTR (En) = el IEN ant? => PERAG) =e ENT ein) — a 

de F(E)=(N Log ele Log(e)= pl" 

= jan Ft) (UN —Uäeng Lol Alain te e Foi 


Em 


On obtient alors les équations séparées suivantes : 


= (+ REG) 22-0 0-26 DEL El asi 


2 4 
are Eliette dt aleet ED E bat 


GE 


= (a (la EtA] e)no 
Solution T(E,n,6)=(E° + 7° IT: E( Ju DIS ) 


4 E 
Hiel ZC, ONE AH, H(n)=0 


2 


n 


Avec le choix de la fonction de séparation comme inverse de la racine carrée du rayon de révolution 
à la puissance N-1, on obtient d'autres équations séparées, il vient la R-Séparabilité-2 : 


Rene Em EE eh a ee RER 


4n 
agep (Natal 20 


č ; lan 
Solution T(É,n,9) = SEE) (6) 
moq) f- EAA, leid SEH 


Exemple de R-séparabilité : le système « hyper-cylindrique inversé » 


Dans ce cas la fonction de révolution et les calculs des fonctions de séparabilité donnent : 


c = é=+in Ge a í PENE c c 
Ein En Xo +P G in) (Ein) (Erin) 


EME INT D ee re UE +in) FE-inŸ =-40 7 
H 


S (+7) 
EE mear au 


De plus la fonction de R-Séparation est choisie comme suit : 


SE-in)= S'(£-in)= 


2 


N-11 


PTR (E,n)= OMR) et p=c z T > R(En)=(E uni: 


> eln 21 [4E F(E) 1 > F (é) ELE (0) 


S p“ R (En) = en ei éi fdn Pl (DN Log (1)+ Logt) = Ra) ES pl 


On obtient alors les équations séparées suivantes : 


EREE E) Gudue22 dra alte SEI, SIE ae 


2 
fait: 
2 


— 


vie Fra) (ulu 220-4402) 


E"(é)-a E(6)=0 g 
ere mrga) {EEA eet Solution E a D (6 va: siet kaiel 
n 2 


Avec le choix de la fonction de séparation comme inverse de la racine carrée du rayon de révolution 
à la puissance N-1, on obtient d'autres équations séparées, il vient la R-Séparabilité-2 : 


N-1 


Se A ren] nE e = AC) A (0-0 


| Cat", 
H"(n)- re 20 a la) 0 Solution reno) E] Sc lit Kate) 


Exemple de R-séparabilité : le système « hyper-Limaçon de Pascal », contribution de René Lagrange 
« Les familles de surface de révolution qui possèdent des harmoniques, 1939, page 303 formules 
26 


Dans ce cas la fonction de révolution et les calculs de séparabilité donnent : 


BEE E EN >x +ip=x,-ip= LE nl flé +in) 


Ge Sinh°(E)- Sin? (0) 
dee Sinh? (§)- Sin? (n)+ 2iSinh(£)Sin(n) _ T A (Cosh(£)-Cos(n) 
3 e (Cosh()-Cos(n)} 2Sinh(£)Sin(n) 


4 (Cosh(é)- Cos(n)} 


c’(Cosh(£)+ Cos(n) 
4(Cosh(£)-Cos(n)) 


| > f'C-in)f'(E+in)= 


nr POSE) EE 


= 7 = M (E)+ zial 


(Cosh(#)-Cos(n)  (fE+in)-f(é-in)  4Sinh? (&)Sin (n) 
1 


> x(é)=- 45n ln) 


1 
45h) | x)=- 


De plus la fonction de R-Séparation est choisie comme suit : 


N-1 p2 ` Jet) fdnF,(n) Ke 2Sinh(£ )Sin(n) 
TE) 
=> R(E,n)=(Cosh(£)- Cos) = oi R? (é,n)= (1) "2" c~ Sinh” (é )Sin™ (n) 
= eln) — 2y Sinh“ (E) et ei) — alt) 
F(E)=(N -1)Cotanh(£) 
> [a F(E)=(N-XL0g SL) | ët 
! Sinh’ (£) 


> jdn F,(n)=(N -1)Log (Sin(n))= F,(n)=(N -1)Cotan(n)  F,'(n)=- 


(N—1) 
Sos" Lal 


On obtient alors les équations séparées suivantes : 


EREE Galus aa DR, D br: 


ve Emil ere 0-40 dai 
o e ulu+2A-1) Lag. 
E" (6)+ MCeenilélSich di Do S dat: 


+24- 1), 
pi in? 


En posant &=a-X communément à&=v(v+1) 


H"(n}+24Cotan(n 


Sinh’ (é 


rs H(n)=0 
Sin? (n 


Solution T(€,n,©)=(Coshlé)- Col "lz läis) 


mi 
E"(£)+ 24Coanh(£) "(€ fe. u(u+22- h n)=0 
La 


H"(n) + 24Cotan(n 


Les solutions sont identiques à celles du système hyper-sphéroïdal allongé et inverse hyper- 
sphéroïdal allongé. Avec le choix de la fonction de séparation comme inverse de la racine carrée du 
rayon de révolution, on obtient d'autres équations séparées, il vient la R-Séparabilité-2 : 


N-1 


st (Cosh(£) Cos( ) D 8 1 1 8 8 
nes [MO BOT ge aile ein AO AUS 
= | WSA 2e) 0 na 

in (£) | (Cosh(£ )- Cos(n ) Sé 
Solution T(É,n,®)= Së to oe) 
H"(n)+ |- (u + ii -1) e dn -0 | Sinh(£ )Sin(n ) | 


Les hyper-cyclides de révolution 


Pour les systèmes de coordonnées basés sur des cyclides de révolution, je commencerais par le 
système d'E.G.C. Poole conduisant directement à l'équation de Lamé. La fonction de révolution est 
la suivante : 


To a Wei zl 


nz yasn(u,a) a len (v B) 


G=u+iv 


dal. o Maly. B)-acn(u,a) 
BlL-dn ek 8 


Juve hägielkt- vgl Je, geil ei 81 Kla 


ali ` 2 acn(y,a) 
a rh (v, 8) 


1/2 
oll-a*sin? (0)? K(B)= | aoli- p°sin? (o) 


0 


p=ycn(v, p) 


E 
e 
0 
Le calcul de la R-Séparabilité de l'équation de Laplace donne : 


dn(u,a) 
f'(u+iv)f'(u-iv) = me GG dn?’ (ma) _ dn’ (v, B) x(u )= cw J 
Guivy- TTT "3" atA] PS T 


Et le système d'équations séparées en utilisant la fonction de séparation « canonique » Ly ol"! 


re 1 EA Je BÜi-an°(u,a)sn°(v, B) i 
Li Ve(av)* SE teen 
(LH) CE end pt 
Srtale au 3 Une 3 — Usine) C Jil 


H'(v)+ l EA OT E ' BU 


ution v = B(L-dn°(u,a)sn°(v, p)) re 
Sohuion Tuv, O) =- EE 


Et sachant les relations entre les fonctions elliptiques de Jacobi, il vient : 


z" (u)+ f (u+4Xu+4-1)a?sn?(u,a)+E}Elu)=0 


HIER (Ee Clad- 


Posons ZC =(u+AXu+A2-1)-6 = 


Une fois développée l'expression snliv, pl), il vient : 


env, p ) 


Ee 
Pri E -(u+4{u+4-1)as (7, Slükl- 


Si l'on pose t=u+À-1 alors on retrouve exactement les équations de Lamé : 


E"(u)+ E-k +1)a’sn?(u,a)}E(u)= 0 
0 


EES 


Système de coordonnées hyper-sphéroïdales 


Nous allons considérer un plan méridien à deux dimensions de coordonnées « polaires » (p,xo) 
plongé dans un espace à N+1 dimensions, sur ce plan « méridien » une ellipse et une hyperbole 
perpendiculaire forment les deux iso-courbes d'un système de coordonnées elliptique, ces iso- 
courbes engendrent par toutes les rotations complètes possibles (révolutions, groupe des rotations 
SO(N-1) probablement) autour de l'axe x, des N-1 dimensions restantes des hyper-sphéroïdes de 
révolution. 


De même si les iso-courbes génératrices du système de coordonnées sont des cercles et des droites 
quelconques passant par l'origine, par rotation on obtient des iso-variétés hyper-sphères et hyper- 
cônes. 


Commençons par l'exemple du système hyper-sphérique dans un espace à N+1 dimensions défini 
par ses coordonnées cartésiennes (et le Laplacien correspondant) : 


LN Së i=N 92 
= x? RE SE =0 


= 


x = (X1, X353 Xp X0) r= 


Sak 


X 


H 


i=0 


i= i 


La distance à l'axe x, définit la variable p : p= S x et la relation : „2 — p°+x? a lieu. 


i=1 
Introduisons les angles des e hyper-colatitudes » et l'« hyper-longitude », comme ceci : 


ER ln -Ñ =ọ) 


Au = r Cos(9, ) He [0, x] 
x = r Sin(S,, Koslä, Iya E [0,7] 
xX, = r Sin(S,, CHEN Koslä, Iy E DÉI? 


Gë 
r zt 


Syan F Sin(9, CHEN )..Sin(9, )Cos(9,) He Io. z] 
Xy = r Sin(S,, )Sin(9, )..Sin(9, )Sin(9,) KEE [0,27] 
4 =¢ > hyper -longitude HA, — hyper — colatitude 


x, =r Cos(9, ) 


On a bien : IZN GE, Be 
p= De =r Sin(9,) S j 
i=1 


Les deux iso-variétés principales à N dimensions de ce système de coordonnées sont les hyper- 


$ i=N A i=N 
sphère Dé =r’ =Cste €t les hyper-cônes : D x = x, Tan(9,) 2 = Tan, ) 
i=l 


i=0 


ue ` A 5 i=N i=N 
La métrique de ce système de coordonnées s'écrit : S'dx=dr BM, avec r? 2 


i=0 i=l l=i 


Les facteurs d'échelle se calculent comme suit : 


h, =1 
hs, =r 
hg =r Sin(9,) IN 
Ze = hhg hg, hs =r" ëmm (9y Sin” (9y ).--Sin(9,)=r* ] [Sir (9 
he. = r Sin(9, ëialä, Ve rho, loy Ja CT Din E ın ( Ka) in( A r IT ın ER 
I=N l= 
SA J [Sin EN Sin™(9,) 
l=N SL i- 
h, =r [| | Sin(9,) ve SE Se a EWEN 
er 3 El Sin? (9,) El Sin’ (9,) 
l=i+1 l=i+1 
I=N 
hg = WEEN 
1=2 


Les termes du Laplacien se calculent ainsi : 


I=N 

Si 1-1 9 

1 ð Je or |, 1 A ô | Je ôr Jg adn oi Lo 

Aal se et =r I=N Sin ER 
"Je ôr h? ôr TR 20 h,” 08, 


h 2 
S [I Sir (9) 


I=N 
[I Sr" (9,) N 
kr) Ol am ô (sus) ÎI 1 i $ DATES 
in 


1 6 | 
2 2 LN I=N 2 se? il 
de 68,1%, 08) r [Isin (8) Tata) "7 Or (9;) |si ™(8,) 09; 
1=2 l=i+1 
Et le laplacien dans ce système hyper-sphérique à N+1 dimensions a donc la forme suivante : 


1 ôf ôT 1 
Anu? = r` Zu Or J r’ GO) 


ô | 
L opérateur différentiel de moment angulaire E ( Sin™ (9, )— 
y OP iffe g 2 2 Is) SE El EN | 


i ð arì\ li S 1 28 
É (T)= SLL À Sin™ (9,)— 
Sé see | Lo REH Sin (9 jll Sin(9,) al SC? 3 


N J j=i+1 
GN 


E Ee D 


1 CHERS ôT 1 
Sin"”"(9 IER EN UE Lech (Af 
N N N d 


Le N-ième terme de l'opérateur de moment angulaire, dépendant de Oy uniquement s'écrit : 


1 ð ER ua VII 3 ») 
09, 


Exemples de Laplacien hypersphérique à 3 et 4 dimensions : 


Gi [I | E DEEN 


| Sin? (9, ) | Sin (9) 89, 
JSN 


2 
2-D>N=1> D 2L 


1 


2 
AN fleece SEA 


2 
4-D >N =3> L(T)- e E 


REES (0,) 09 Sin” (0) Sin,) 09, 


au 


ôT 
E 


| 


RL ET 
Sin? (9, ) 8, 


cn 


Continuons par l'exemple du système hyper-cylindrique dans un espace à N+1 dimensions défini 


par ses coordonnées cartésiennes (et le Laplacien correspondant) : 


= 


= 
X bis.) p= 


X 


i 
i=0 


i=N ig i=N 2T 
de at: E - 
10 OX; 


D 


— 0 


f 5 ; EN e 
La distance p à l'axe x, est l'une des coordonnées p= D x2, la hauteur xo est la deuxième, les 
D 


i=1 


autres coordonnées sont le système des N-1 angles dans l'hyper-sphère à N-Dimension de rayon p : 


1 ð HOT ôT À 
Ayat erer Ze i J+ 2 t 2 L, (T) 
p On ©] On P 


L, opérateur différentiel de moment angulaire Le (T ) = 


oT 


CA 


| 


Que donne les coordonnées sphéroïdales « allongées » dans un hyper-plan méridien dans un 
espace à N+1 dimensions 


Dans un e hyper-plan a méridien p,xo , introduisons deux variables Zu telles que (formules de 
passages du système hyper-sphéroïdale au système hyper-cylindrique) : 


fe = Ae -1 -n?) = f = Cosh(u) É =c Sinh(u)Sin(v) 


Xo = CËT n = Cos(v) x, = c Cosh(u)Cos(v) 


Les formules de passage du système hyper-sphéroïdal au système hyper-sphérique sont les 
suivantes : 


ZE 1-n° | ¿€ =] Lu" 
GER Tan(3, )= E = sn(3,)= fi ln) 


>> 
clé -11-7 = r Sin(9,) r= Jen +E 107 )= JE +n -1 > Cos(9,)= = 


Transcrit sur le système de coordonnées hyper-sphéroïdales « allongées » : 


X = cén 

X = cfle n'es |Cos(9,.) dy, E [0,x] 

X, = JE n'es Lead 9, . )Cos(9,.,) dy, € [0, x] 

Xy = fle ne )Sin(9,. )Sin(9, ,)..Sin(9, )Sin(9,) à =p € [0,27] 


Les facteurs d'échelle se calculent comme suit : 


2 2 
h; =c S SE 
d E Ni) 2 2 2 2 ne, . 11 
Je AA ho, A = 0E -pe 1-1) : [Sin (a, 
dei Lu" g N-2 l=N-1 


NO a aN as] | 1- Hues Si l-1 9, 
A E) SSC LG erhal av 
el rale, JS 8 Len pe 0-7) à TL sn (9) 


Lu I=N-1 


h, =c Ve 1-17) TI Sina HEWER 
a CN n II in(9,) d e welt Aaf E 1 Sin (9,) 


l=N-1 l=i+1 


hy = EE Abol J Sin(9,) 
1=2 


Les termes du Laplacien se calculent ainsi : 
ol de or |, d o [sar], 1 0e oTr 
Ve El AT 06 Rx D Se | 
EE l l ofe- 
de OER OE N e (em fer Op 
Ve ôr äi 1 Ô HENSEN 
Ve mi on) e 
ô Je ôr d 1 1 IJ Gleeft 
Ve 58, hy, 05, ei i i 


Et le laplacien dans ce système de coordonnées hyper-sphéroïdales « allongées» à N+1 


dimensions a donc la forme suivante : 


ses 27 Sie EZ, 


ce? Ee 2 


= J 


1 1 ; 
+ e? le fe) 0) 


L, opérateur différentiel de moment angulaire 


. 1 d SH oT 
Le. (T)= KS IL 10 ) Sin (9,) 9, (sin DEA 


i=1 FN 


Normalement on retrouve bien le système sphéroïdal « allongé » pour N+1=3, soit N=2 


a G ar) e D} ae AP: 


C 


Séparabilité de l'équation de Laplace et de l'équation des ondes, sans et avec dépendance 
angulaire 


Examinons maintenant la séparabilité de l'équation de Laplace dans le cas où la fonction scalaire 


n'a aucune dépendance angulaire. Dans ce cas, en introduisant la variable de séparation adéquate 
w , la séparation des variables donne : 


ht Ak) We 7 Sie Tr) Elte 


c? é? =n’ OË 06 ôn n 
ler fN He- Ce ENEE DE 
OP ôg ôn ôn 
a © z 2), OH = 
P Dal Zil n ) ôn +oH=0 
2 Saz 2 2 GC — DE = 
e AL 1) z) @E 0 


En posant À=(N-1)/2, le système d'équation est assimilable aux équation différentielles des 
fonctions de Gegenbauer : 


ge £ 2 
t-r)” Lht- l 2 kan EE "ais Set = 


2 GE 22 z dE = = — 22 dE dE JE 
(£2-1) CL 1) SS wo =(E (baue ous — 0 


Pour l'équation des ondes, avec une variable de séparation w également et une variable de 
fréquence, toutes deux bien choisies, il vient : 


AT V0 G SE 1G aye Thr)” Ze? er -0 


DN CE dE ôn ôn 
(8-07 afe- Zu-ef Ziel Zu TEE 
` Laf? Ze? Z ke t-le -r EE -errn ENT 
ô dE 


124 24 GE es Ne 2 CO fs le: 
kd le EE Aen e-Sapr ue e-7( eo 


Supposons maintenant une dépendance angulaire quelconque dans l'équation des ondes qui 
conduise à introduire une seconde variable de séparation, il vient : 


T= E: Ji pel Dës EECH eg z E()A (nE JO Ayal + Lr =0 


SS ee 2 ale? Alte DN sl ES ei 


+ RE Elke (SEL teg Elte OUT, een Meckel 


œ ) Hh ôn ôn y) (g E 
Leif së ie geg 
n -1 
à. dE Da ff GE u(u+22-1) = _ 
L Jp edea , PE DE k- 


L(8)+u(u+22-1)0(6)-0 


Ce qui établit des équations différentielles que je pourrais nommer d'hyper-sphéroïdales. 


Étude de l'équation hyper-sphéroïdale issue de la séparation des variables de l'équation des ondes 
dans un système de coordonnées hyper-sphéroïdales à N+1 dimensions (N=2+1) 


Est-il possible de construire des solutions de l'équation hyper-sphéroïdale à partir des fonctions 
associées de Gegenbauer : 


1— x? 


(—x2) "(x)-(22 +1)x viele y°{1-x2)- nulu 2A "ai: (0) 


Les fonctions de Gegenbauer de première et deuxième espèce sont solutions de l'équation 
différentielle : 


(22 )"(z)-24+Dzy'(z)+vlv +24)y(z)=0 


Les fonctions associées de Gegenbauer de première et deuxième espèce sont solutions de 
l'équation différentielle : 


l-7) e)-C+ Daye) {vy 2a) EE 


Lorsque u est entier, soit u=m, on peut construire la fonction de première espèce comme suit : 


GIE ME 


ee 
2 


On s'assure que pour À=1/2, on retrouve la fonction de Legendre associées : 


MS f atm) ufr pm 
c? (z)=(-1) EE (z)= Pr 


Toutefois cette dernière définition peut en apparence comporter un pôle en z=1 lorsque À>=1/2. En 
réalité il n'en est rien car on démontre que la limite est bornée (voire nulle notamment lorsque 
À>1/2): 


1-24 
TA 1-24 ljm 
iim Cae im (CNE VAT + ar. y2) p7 (x) 


= 0 
PER DER T(v-m+1)(2) 


Aw 
2 


Pour éviter numériquement l'écueil de z=1, on peut aussi définir cette même fonction par la 
formule suivante avec la fonction hypergéométrique : 


SE BEIENEE af rf ; z 


= 1 Sé u,v +2 + u; — +4 + u, —— 
2“T(24 iv uar Z+a+ n) i 


2 


qui cette fois ne présente plus de pôle apparent pour z=1 et Um entier. Et l'on s'assure qu'elle 
redonne bien les valeurs des fonctions de Gegenbauer et les fonctions associées de Legendre : 


1 
als c'e Erste see 22 


5 l-z 
l-z k F|- 1 + ; 
E 2) 2 { v+u,v+I+ul+u; S | 


T(22+v) 1 „l-z 
=0 => C’?! (z)=Cî(z)= F| E 
u => v (z) 7 (2) (22) (v+1) * { Hurt > * 2 | 


Pour la construction des fonctions associées de Gegenbauer de deuxième espèce, je vais me 
restreindre au cas où le paramètre Um est un entier. Dans ce cas la fonction de deuxième espèce 
se déduit d'une formule de dérivation de type Rodrigues. Je donne sans démonstration le résultat 
de cette construction à partir des fonctions associées de Legendre de deuxième espèce : 


cs (= Er a) gr) 


Il se trouve que pour les conditions générales de paramètres u, v telles que u+v est un entier 


u 
A . À DÉI Z sp 
négatif alors la fonction associées de Legendre Q; ( ) n'est pas définie. Dans le contexte de la 
construction des fonctions associées de Gegenbauer de deuxième espèce, lorsque v-m est un entier 


ces dernières ne sont donc pas définies : t À TE WEE TE L em 0- Afin de simplifier la 
2 
suite de l'exposé, je supposerais donc que V 2". 


Les fonctions associées de Gegenbauer de première ou deuxième espèce satisfont à la relation de 
récurrence suivante : 


z Ca (2)= EES EE (22 +m+v -1)C}7 (z) 
í 2(4+v) 

ez th beim, Lelh- DÄ mu UE, el 
EC 2(4+v) 


En itérant cette relation juste sur les fonctions de première espèce, une fois de plus, il vient : 


(m+v+21-2Xm+v+21-1) am (z)4 
(2v +24)2v +22 -2) Wë 
ke 2v(v+22)+(2m-—11-22)-2m° -1 eich 
(2v +24 -22v +24 +2) ‘ 
2 Lon —u — (lan —u 2 c>” (z) 
(2v +24 (2v +24 +2) 


Ch 


Une des formules de dérivation « vers l'arrière » des fonctions associées de Gegenbauer est la 
suivante : 


A,m A,m 
Cäatch- KZ EE (z) 


Une itération supplémentaire de cette formule permet de calculer la dérivée seconde : 


Ci (= v(z(v-1)-1)C2" (2) 2 (2v -3X24 +m +v —1)CŻ4" (z)+ (24 +m +v -1X24 +m+v - 2)C>" (2) 
v Weck 


Mais pour éviter que nous soyons dans le cas V—mMm< 0 pour les fonctions de deuxième espèce, il 
vaut mieux utiliser la formule de dérivation « vers l'avant » suivante : 


came PAGE") mrC 
TT 


Une itération supplémentaire de cette formule permet de calculer la dérivée seconde : 


E (24 +v ](1+ 22(1+ 24 +v)]C2"(z)+z (m—v—1)42 +2v +3)C2"(z)+(m—v —1Xm-v-2)c/2 (2) 
C; ( ) D 


On sait également que la dérivée seconde peut astucieusement être calculée à partir de l'équation 
différentielle , soit : 


Ko 


CAGE = 0 


Ets 


Ce qui permet numériquement d'effectuer une vérification de cohérence. 


Ga+1) TORT 24) ”0"+24 Dle $ 


Toutes les formules précédentes sont identiques avec les fonctions associées de Gegenbauer de 
deuxième espèce. 


Cas particulier des fonctions hyper-sphéroïdales de paramètres À demi-entier 


Pour la construction de ces fonctions, on peut utiliser une autre forme des fonctions associées de 
Gegenbauer qui ne sont définies que pour les valeurs À demi-entières, comme suit : 


V (4) DEEN 
Au 7 3-1 22 2 ee à SA 
cé) NE) OCH 


Cette définition étant valable quelques soient les valeurs de u, et notamment lorsque u=1-À. 


Construction des fonctions hyper-sphéroïdales anqulaires de première espèce de degré et d'ordre 


entier 


Sur l'intervalle [-1,1], les fonctions angulaires de première espèce de degré et d'ordre entier 
peuvent se construire à l'aide d'un développement en fonctions associées de Gegenbauer: 


xe[-11]  Aps/"(x,y2)= ` (1 azr cza) 


Kom 


Fonctions associées de Gegenbauer C7, (x) A e [- 1,1] 
Les coefficients du développement se construisent à l'aide d'une récurrence : 


Adi" (2)+(8,-œ(7°)d/"(7°)+ cd" (72)=0  kef-0,..,0,1...,+0} 
A, Aë (n—m+2k—1{n-—m+2k) 
(2n+4k+22-2)\2n+4k+22-4) 
2 (n+24n+2k+22)+m°-1+(22 EL EEN 
(2n+4k+22-2)2n+4k+22+2) 
2(n+m+2k+22n+m+2k+22 +1) 
(2n+4k+22+2)2n+4k+22 +4) 


B,=(n+2k\n+2k+22)-2y 


C; 


Commentaire sur la récurrence : 


Lorsque le paramètre À est entier alors les coefficients A, Cx et Bx de la récurrence présente des 
pôles. Pour l'instant je traiterais ce cas par la suite. Lorsque le paramètre À n'est pas un entier II n'y 
a plus de pôles pour ces coefficients. Mais il demeure des cas où les coefficients A, C, s'annulent et 
dans ce cas la série est limitée dans l'un des sens des indices. Je dirais que la série est unilatérale 
dans ce cas (vers l'infini ou moins l'infini). 


Exprimons cette récurrence sous forme de deux fractions continues équivalentes sur les rapports de 
deux coefficients successifs : 

d'Zil orda) 

eng y’? 
dch) À, 


> 
di" (y) zalet, Se? 


Notons f, ©) = 


A 
B,- oly? )+ GEI 


k>0 


d'"() C C, 
k<0 a = k Œ k ʻi ñ Go 
< A, GL S EE a) & f;, B, zae RATO fraction continue montante 
í aeh 
d? m(y2) 1=k 8 da” GL 1=k 
k n,k = (-) n,—k z (+) 
> 0 ar y? IT di y? l JS 


fraction continue descendante 


Annulation du terme Ck 


Cette annulation n'est possible que si le paramètre À est un demi-entier. Dans tous les autres cas 
excepté À entier (voir cas à part), la série est bilatérale, soit se développe des indices -°° à +0», 
Lorsque À est un demi-entier, soit À=(2p+1)/2, C, s'annule à l'indice k,-1 alors les coefficients du 
développement s'annulent de même, pour les indices k < kọ. Lorsque C, s'annule en kel la valeur 
de ko est la suivante : 


n+m+2À+1l n+m+2À 
k, —1= T 
C,=0& (n+m+2k+24{n+m+2k+24+1)=0 2 ou 2 
n+m pair n+m impair 
+m+24-1 +m+2À-2 
bss bss n+m+24-1 
© 2 ou 2 S k,=- S 
n+m pair n+m impair 


C'est bien notre cas lorsque l'on prend n >=m, alors : 


n+m+22-] n+m+22-1 
E = > —k, = 
2 2 


Mais il faut prendre garde à ce que certaines des fonctions de Gegenbauer associées sont nulles, 
lorsque n, m sont entier et À demi-entier, car on sait que : 


v,u,peN A C##(z)=0 -u-(24-1)<v<y 


2p+1 
2 
=> Cn (x)=0 -m+(22-1)<n+2k<me-m-n+(22-1)<2k<m-n 
A,m 2 REES = KR 
Ainsi les coefficients 1 (°) / -n-m=+(22-1)<2k<m-n 
la sommation ne commence donc qu'à l'indice : 


ne sont pas utiles. Dans ces conditions 


kọ =—R E —— 5=(n-m)mod 2=(1+m)mod 2=> R =| 77| 


Et le développement se présente sous la forme : 
seji] aps” )= Satro) 
r=- | 
Remarque : en changeant les indices de sommation du développement, il vient : 


Si n-m pair > k= ~ n+2k,=m Pak = Ain art Stréi, (y Sien C” al d Zeck és (x) 


1=0 


n-m +00 


Si n-m impair > k) = 2 — n+2k,=m+1l l=k k, > hps?” (x,y?) = (-1 (Tac, (y les E d Zeck los (x) 


n eh, W n,l m+21+1 
1=0 


Ces deux expressions peuvent être données avec une notation suivante pas toujours très claire : 
+00 
De amà 2)" (x ) Si n-m pair 


1=0,1 "(y 2) (x) Si n-m impair 


Annulation du terme A, 


Et lorsque A, s'annule à l'indice ko+1 alors les coefficients du développement s'annulent pour k > ko. 
Cela donne : 


een kajo 
Ae (n-m+2k-1\n-m+2k)=0<44 ° 2 ou 4 ° 2 
m—n impair m-n pair 
m-n 1 mn 
k, = —— k, = — 1] — —n+6 
© 2 2 ou g 2 EE l< h=” S 1 avec =(m-n)mod2 
m-n impair m-n pair 


Et le développement se présente sous la forme : 


mn) 


LJ 
2 
sel pr) 3 Cargo) cal) 
Ae 
Mais là encore un certain nombre de fonctions associées de Gegenbauer sont uniformément nulles 
lorsque m et n sont entiers et À demi-entier. 


Comme Ee (x)=0 -m-(24-1)<n+2k <m <& -m-n-(24-1)<2k<m-n on a la sommation 


m+n+22+1 


commençant à l'indice Kg SC 2 | soit un développement de la forme : 


k=- 


xebh] psy) D (Calc (x) 


k=—0 


E 


Toutefois ce cas de figure du développement unilatéral vers les indices négatifs est totalement 
équivalent au précédent en raison de la relation de transformation suivante du degré n : 


Sin(xv) Au 
Sin(x(22+v)) "7 


1 
c*" 
Cos(21x) "'" œ) 


p EN > C)" (x)= CH, (x) => hps?” (x,y? )= hps?” (x,y?) 


mat) 


x e [+11] c$” (x)=- (x)= n,m € N > C?” (x)=-— 


_2p+l1 
2 


> À 


pA 
=> hps? (x, 7° ) = KS (- 1) dioak LG ) SRE (x) Comme Css (x) en (x) 
fe 
A k=+00 P 
=> ne (x, y? ) = > (- 1) He (7° ) ee (x) 


TA m-—n+l =| nm 
2 2 


Cette dernière expression donne un développement unilatéral équivalent vers les indices positifs. 


Développement potentiellement fini lorsque -m-2A<n<m, n,m entiers et À demi-entier 


Dans ce cas les coefficients Cx et Ax s'annulent. A,s'annule en k;+1, la valeur de k- est la suivante : 
k, = Ck s'annule en k,-1, la valeur de kı est la suivante : D E Et le 


2 2 
développement se présente sous la forme : 


A m—n—1 


A e [- 1,1] hps?” (x,y? )= 5 GA SE VIER (x) 


p=- ntm | 
2 


Or toutes les fonctions associées de Gegenbauer impliquées dans le développement fini sont 


uniformément nulles pour n,m entiers et À demi-entier : 
Cr” (x)=0 m (24 1)<n+2k<me-m-n-(21-1)<2k<m-n 


n+2k 


Car jam (m+n) (24 1) et J em nr 


Ce n'est donc pas essentiellement le développement qui est utile, mais plutôt la détermination des 
valeurs propres w qui sont alors racines d'un polynôme. La matrice de détermination des valeurs 
propres est une matrice de taille finie entre les indices : 


Knin = — ER Kan = LL 
2 2 


Étant donné que la détermination des valeurs propres est indépendante selon la transformation m- 
>-m+1-22, prenons la valeur maximale de m et de -m+1-2À. Par ailleurs la taille de la matrice est de 
nf) etui 22 —1 

2 2 2 Ily a donc m+(2ì-1)/2 valeurs propres de cette 
matrice. Ces valeurs propres sont indépendantes de n, qui ici ne joue qu'un rôle d'indice, compte 
tenu du fait que le système linéaire est identique par la transformation n->-n-2À. Une quelconque 
des matrices avec les paramètres respectant -m-2A<n<m donne le même déterminant 


De M -old ,,, | =0 
F2 et donc le même polynôme pour la détermination des valeurs propres. 


ne AT EE 
Commençons par la valeur 2 et posons par simplification 2 
AA core | j_24-1,24-1 7 
Dans ce cas y =— 2 2 DES =0 er, 2 2 -o0 et 
mın 2 2 max 2 
1 2 1 2 
M — © Id, = Aenzeien A 
la matrice se réduit à un scalaire : 4 4 


m 
Si nous prenons la valeur : 2 et posons par simplification 2 . La 
matrice est de taille 2x2 : 


|1 
g7% t27 -o —27° 


2 9 2 2 
KEE -2y -øo 
4 D 


M — © Id, — 
2y 


L'annulation du déterminant donne la recherche du polynôme suivant : 


EE —2y? 5 5 9 
Det || 4 9 BEEN }+# +47 + =0 
27? Lo 2 2 16 
4 
Soit les deux racines w : œ, = Zx EN Ss -2-2 ra 
| E GE 
Si nous prenons la valeur : 2 et posons par simplification 2 La 
matrice est de taille 3x3 : 
SEET 67” 
M -o Id, = 4y? Ze A em e 
0 2 25 4 , re 


Construction des fonctions hyper-sphéroïdales anqulaires de première espèce de degré et d'ordre 


non-entier et À quelconque 


Considérons maintenant les deux paramètres n et m comme non-entier et notons les v et u. Sur 
l'intervalle [-1,1], les fonctions angulaires de première espèce de degré v et d'ordre u peuvent se 
construire à l'aide d'un développement en fonctions associées de Gegenbauer de degré v et d'ordre 
u: 
xep] Hs) $ azh er) 
k=—%0 


Fonctions associées de Gegenbauer C7, (x) xE |- 1,1] 
Les coefficients du développement se construisent à l'aide de la même récurrence avec les deux 
paramètres v et u : 
adit g) (B, -o0 azto )+C att 0 )=0  kef-c0,..0,1,.,+00) 

5 (v u+2k 1Xv u+2k) 
r Qvrak+24-2K2v+46+24 4) 
2(V+24)v+2k+22)+u?-1+(22-1{u+2) 

(2v+4k+22-2)2v+4k+22+2) 

 (V+u+2k+24v+u+2k+22+1) 

(2v+4k+22+2)2v+4k+22+4) 


KI 


B, =(v +2k\v +2k+24)-2y 


C= 


La condition pour qu'une telle récurrence puisse être définie quelque soit l'indice k positif ou 
négatif est que v+À ne soit pas un entier . 


Voyons les conditions lorsque -u-2A<v<u pour lesquelles les coeffients C et A, s'annulent 
simultanémement. Si les coefficients C et A, s'annulent, alors A, s'annule en kmas+t1, avec 


k 2 et C, s'annule en Kmin-1 avec x. = | Cela implique donc les deux 
max 2 min 9: 


contraintes Juve 2Au+A)eN 
v+u+24-1EeN 2(v+1)eN 


Et la détermination des valeurs propres w qui sont alors racines d'un polynôme. La matrice de 
détermination des valeurs propres est une matrice de taille finie entre les indices : 


k = u+v+22-1 k = u—v-l1 
mın 2 max 2 


Étant donné que les valeurs propres sont identiques suivant la transformation u->-u+1-2À, prenons 
la valeur maximale de u et de -u+1-2À. À n'étant pas un demi-entier, la taille de la matrice est de 


N= K Ee d k 8 +v+24 = +1 qui est également le nombre de valeurs propres. 
2 


2 


Une quelconque des matrices avec les paramètres respectant -u-2A<v<u donne le même 
TS Det(M -oœ Id, )=0 


déterminant et donc le même 


polynôme pour la détermination des valeurs propres. 


Commençons par les valeurs ,,} 3 - ut, Dans ce cas x; =k oer la matrice se 


S 
5 ` 28 


A7 J >o E Si nous prenons la valeur : ,, + à = 


1 
A 


H=V +1 Ja matrice est de taille 2x2 : 


réduit à un scalaire ` M - œ Id, = et 


2 
2 


M — © Id, = 


1 DH z DH DH e 5 
L'annulation du déterminant donne les deux racines w : on = EE + 1-47? : 


Si nous prenons les valeurs 4 a= u =v +1, la matrice est de taille 3x3 : 
9 2 2 2 ] 
2 + — 47" -o 67 0 
1 16 4 
M — © Id, = 4y? —+—y A -o he 
£ d 4 3 SC 
ó EE Ee 
3 A 3 | 


Évidemment ce sont les même valeurs que l'on obtient par rapport aux valeurs des paramètres u et 
v entiers mais ici dans un cadre plus général lorsque „+ et donc, — 3 sont demi-entier. 


Dans le cas où ,,, 1 est entier non nul , ,; 3-0 positif ou négatif „<o et entier alors la 
matrice est de taille supérieure ou égal à 2x2 et le parcours de la récurrence comporte des indices 
k négatifs ou positifs qui provoque la divergence au pôle du terme en question. La seule manière 
pour que la récurrence ne comporte pas de pôle est : 


l—-2p-1 1-1 1-1 
Knin = = H P = Kna = =0—=7/-1,2 = 0 
min | 2 | | 2 | P max | 2 | P 


Toutefois lorsque , ; 1-0 u+a=1, ON voit plus tard dans ce document que la récurrence se 


simplifie et que les termes C, et A, ne s'annulent pas mais demeurent constant. La récurrence est 
infini dans les deux sens, elle est similaire à celle des fonctions de Mathieu, qui servent à construire 
la solution de première et deuxième espèce. 


Cela confirme donc le fait que lorsque v+ est un entier la récurrence n'est pas défini et ne peut 
servir à la construction des fonctions hyper-sphéroïdales 


Normalisation des fonctions hyper-sphéroïdales de première espèce et de degré et d'ordre entiers 


En l'occurrence la normalisation proposée des fonctions hyper-sphéroïdales de première espèce de 
degré et d'ordre entier est calquée sur la normalisation des fonctions associées de Gegenbauer qui 
sont asymptote pour le paramètre w->0 : 


+l 24-1 
dzl- Des CZ f = 21-22 T 2 T(n+m+22) 
j a | í ©) (GI 24+2n T(n-m+1) 


Cela se traduit par la normalisation suivante : 


ya 2 T 2 m T(n+m+22) 
dxlhps*"” 2 =a! jA 
| aclios; Gr) 22+2n (C(A) T(-m+l) 


Qui donne la condition suivante sur les coefficients du développement : 


n se E 2 ze Tàn+2k+m+24) n 2 m T(n+m+24) 
delt? hps?” 2)_ dr" (y2 f2! 24 = 2124 
f "17 msi”) 2 EE CAY Dia 2t ail 24+2n (TA Th-m+1) 
a ge Ga) 1 T(n+2k+m+24) 1 T(àn+m+24) 
A 2A+2n+4k T(n+2k-m+1) 24+2n T(n-m+1) 


ef) 24+2n I(n+2k+m+24) L(n-m+1) 


=] 
VW 53+on+4# T(n+2k-m+1) T(n+m+2À) 


Évidemment lorsque le développement est unilatéral, la condition de normalisation pour la 
fonction de première espèce s'écrit sous la forme : 


e Lost mm 24+2n Tfn+2k+m+24) T(n-m+1) 


=] 
22+2n+4k T(n+2k-m+1) T(n+m+24) 


Construction des fonctions hyper-sphéroïdales anqulaires de deuxième espèce de degré et d'ordre 


entier 


Sur l'intervalle [-1,1], les fonctions angulaires de deuxième espèce de degré et d'ordre entier 
peuvent se construire à l'aide d'un développement en fonctions associées de Gegenbauer de 
deuxième espèce dont les coefficients suivent également la même règle de récurrence. On part 
d'un développement de la forme : 


xef-11] rgs?” (x,y?)= StG-Ufdiälzkän, tz 


k=—0 


Fonctions associées de Gegenbauer de deuxième espèce CY nax (x) xE |- 1,1] 


Deux cas peuvent alors se distinguer. 


Le premier cas correspond aux valeurs du paramètre À qui ne sont ni entières, ni demi-entières et 
les valeurs de n et m sont entières. Alors les conditions pour lesquelles les fonctions associées de 


Gegenbauer de deuxième espèce ne sont pas définies sont: "MEN  n<m_ Comme on se 


place justement dans le cas n>m alors le développement précédent est valable quelque soit les 
valeur n>=m et il est infini dans les deux directions. 


Le deuxième cas correspond aux valeurs du paramètre À qui sont demi-entières et les valeurs de n 
et m sont entières. Mais si dans le cas des degrés et ordre entiers, la liste s'arrête à partir d'un 
indice négatif ko pour la fonction de première espèce : 


xebh] ` Speis, E Carree) 


n+m+22-—1 
ver JE 


il n'y a aucune raison qu'il en soit de même pour la fonction de deuxième espèce. Aussi doit-on 
prendre un développement infini dans les deux directions. Mais pour les valeurs d'indice pour 
lesquels : 


m+n+2k+21-1eN 
m+n+2k+22-1<0 


Les fonctions associées de Gegenbauer ne sont pas définies. On constate que la limite de définition 
des fonctions de Gegenbauer associées de deuxième espèce coïncide avec la limite des indices de la 
série pour la fonction hyper-sphéroïdale de première espèce. Aussi lorsque l'on sépare le 
développement en deux parties : 


k=-M = +00 
hqs” (x, y’ ) es 5 (- 1) CH GL: fer. (x) + 5 (- 1) diz (7° EE (x) 
k=- M 


Kom 


m =| 25m2 


Le premier développement n'est constitué que de fonctions associées de Gegenbauer de deuxième 
espèce non définies. 


Prenons alors cette première série et réalisons formellement un passage à la limite : 
k=-M =] k=-M -1 
Lin DOY aE O n) De aze o Kë ts 
k=—0 Ae 
=> —v(e)-24 =n+2k +8 >v(e)=-n-2k-24-e€ 


L'astuce consiste à utiliser les formules suivantes de transformation des fonctions associées de 
Gegenbauer de première et deuxième espèce pour les domaines de valeurs des paramètres v, u et 
A. pour z € [-1,+1] : 


2 1 
v,ueR =P p eN bs Ge) 
2p+1 m Cos(u z )Cos(v z )C"(z)- Sin((u +v yr KT, 
„uer 2-222 pen cog (er EU aost ek sav} t 
Silo ve) 
Cela donne : 
cs Gë x Cos(mr)Cos(a(-n-2k -24 SE a Lrl- Sin(z(m-n -2k -24 ell, 1. (x) 
(Ohn+2k+e Y Sin(x(m+n+2k+e+2À)) 
Cos(mr)Cos(r(n+22+e)) nı Sin((m-n-22-er) 
CE" Nä C „m cy ban 
> On) m Sinf(min+24+e}r) Ps tel Sinf(men+24+e}r jE Gg E E 
Cos((m+n+22+e)r) à Sin((m+n+22-£)x) à 
Can Es „m „m 
> Gaane) T Sin((min+24+e}r) EE Sin((m+n+2A +e) sis) 


(O)-n-2#-2à (x) 


H 


Am 
z es (x) ES Sr G) + Cå 9 
Injectons cela dans la série et effectuons un changement de coefficient par passage à la limite : 


k=-M M 
um SC dark Ad Lam "8 CN EU (our, (De can) 
Aux?) zu 2) ; Ke 
CC d,x V => Lim ZC E ajr (y DOSAO be 1) d? "(y Se 2k- aax) 


Supposons Lim 
£e—0 E€ 


L'expression du développement de la fonction sphéroïdale de deuxième espèce devient donc : 


Der SL 24-1 
has") Cali tel Ya ka (x) avec M = 7 tnt | 
em LA 


Il reste à déterminer la méthode de construction des deux jeux de coefficients et la normalisation. 


A,m 2 
Les coefficients di ( ) sont construits de la même manière que pour la fonction de première 
espèce, à l'aide de deux fractions continues descendante et montante : 


R d?” y’ X À 
fi J z SS CG fe ) B — oly te Tak k € D. koax] 
n,k-1 k k+l 
d?” (°) C 
(+) — Zei (+) f k AM 
fe SC y? f B, = oy°)+ AEN | l 


SCHT 2 
Les coefficients di ( ) sont construits à l'aide de la fraction continue montante, en commençant 
à la valeur nulle pour un terme d'indice extrémale suffisamment négatif. 


Ka C, — 
i B, — oly? )+ A E 


Cette fraction continue se connecte à la fraction continue ascendante comme suit : 


ga o(n+m+2k+24n+m+2k+24+1)_ o EE E A mt n+m+2À-1-06 
(2n+4k+22+2)2n+4k+22+4) 2 2 
2(k,-1)=-(n+m+21+1-6) 
n+m 21. 1) 2A=Ee+0-1 
n+m+2(k, -1)+24+1=E€ +6 avec 8 =(n+m+22-1)mod 2 
2n+4(k,-1)+24+2=4(k,-1)+44-24+2= 2€ -2(m-6)-24 
2n+4(k,-1)+24+4=4(k,-1)+44+2-22+2=26-2(m-8)-24+2 
Cun = (e+ -IXe +8) de pl 5-1Ke+5)=-E(-1) 
ANT DE, TR Gén GS) 
0x y? 
+24-1 Ci 
y en Ci ; (-1) Gen SEENEN 
C „u = Lim Me y 3 
eu € (2(m-86)+22)2(m-6)+21-2) NEE 2 
n+m+2-l impair C,,= 
(2m+24-2(2m+24-—4) 
le, ZK Cas 
E let "Benk 
= Lim sat) i =-Lim MA S = — aa) Cu 
Se E aprh) St E RETTEN dé", (y gd Buai- Di )+ A À y eech 
SET DI AE C, 8 
axb’) a T B, -o )+4 7È 


La normalisation des fonctions de deuxième espèce : on utilise le même type de normalisation que 
pour les fonctions de première espèce sur les coefficients des indices k€[-M,+ce)], soit : 


+00 


> 


1= 
© 


- 


+00 


den | DATE ET: 


(a ( dd 2A+2n T(n+m+24+2k) T(n-m+1) 


n,k 


ma AW P IA+2n+4k T(n+m+24) T(n-m+1+2k) 


(n-m)! 


[amra 


at am 24+2n (n+m+24-1+2k)! 
ezte) (n+m+24-1)! 


2 


Que l'on répercute dans le développement : 


Der Ai 
has" (x, 72)= > 1) cle es 22 (x)+ LV dir E avec M == | 
k=—0 


[eme | 
2 


Lo 


amf am 24+2n (n+m+2à), 
(412) 
24+2n+4k (n-m+1),, 


2 


A titre indicatif les indices du premier développement peuvent se calibrer sur un autre intervalle : 


e (x) ke | = 


—n-m—22+1+06 —-2 
2 


| avec Ô=n+m+21-1l mod 2< Ci” se) k'e Del 


Construction des fonctions sphéroïdales radiales de première, deuxième, troisième et quatrième 


espèce 


Les fonctions radiales sont en fait les solutions de la même équation différentielle sphéroïdale mais 
lorsque l'argument réel x est supérieur à 1 : 


(ha) (+) Hale 


(EE: 


On peut transformer très légèrement cette équation différentielle pour tenir compte de l'argument 
x>1: 


e El (2 +1,40) Lo y(x 1)+ Ge D), OS 


Ensuite on transforme cette équation différentielle par le changement de variable suivant : t=y x, 
puis le passage à la limite y->0 : 


702 r DO at ue 0) +(2-2)y(x)=0 


Or la dernière équation l'équation des fonctions de Bessel définie à partir des fonctions de Bessel 
sous la forme : 


azry) Ale Zeniti 06420) 0-0 8 aide cl eich 


Le passage à la limite y->0 des solutions radiales de l'équation hyper-sphéroïdale est donc : 


Lim À t | =a Jusa DO +b Lal 
y t t 


7—0 


On note également que l'équation différentielle hyper-sphéroïdale est invariante par le 
changement : m->-m+1-2À. Cela signifie que les solutions de l'équation différentielle pour m le sont 
également pour -m+1-21. 


Ce comportement limite en y->0 nous suggère d'employer les fonctions de Bessel comme une base 
de développement en série pour les quatre types de fonction radiales : 


veel Hs” (x,y? )= Ji LI Ka behal SEN 


Kom 


Fonctions de Bessel Ian (x) J=J, Y,H0,H® 


Quand on injecte ce développement et compte tenu de l'équation différentielle des fonctions de 
Bessel du développement, les relations de récurrence et de dériviation on arrive à la récurrence 
suivante à trois termes sur les coefficients r : 


Ari (2 )+ (6, -ol än (2 )+ Gr" (7° )= 0 
EE (n+m+2k+22-2{n+m+2k+2-1) 
(2n+4k+22-2)2n+4Kk+22-4) 
o (+24) n+2k+22)+m° -1+(m+2)22 —1) 
(2n+4k+22-2Ÿ2n+4k+22+2) 
(n—m+2k+1){n-m+2k+2) 
(2n+4k+22+2)2n+4k+22+4) 


B,=(n+2k\n+2k+22)-2y 


Cr 


On constate donc que c'est en fait une récurrence similaire à la récurrence obtenue pour les 
fonctions radiales avec les changements suivant m en — m+1-21 : 


Ad?" (y°)+(8, -o0 laien  kef-,0,1...+0) 
= (n—m+2k—1){n-m+2K) 

(2n+4k+22-2\2n+4k+22 —-4) 
2(m+2k\n+2k+22)+m°?-1+(22—1X{m +21) 
(2n+4k+22-2Ÿ2n+4k+22+2) 
2(n+m+2k+22Ïn+m+2k+22+1) 
(2n+4k+22+2)2n+4k+22+4) 


B,=(n+2k{n+2k+22)-27 


C, = 


et les coefficients A et B changeant de signe. Ce dernier changement peut facilement se résoudre 
par une multiplication des coefficients par le facteur LI Ainsi les coefficients de la série radiale 
sont liés aux coefficients de la série angulaire par la relation : 


a} CY aioe) 


La série radiale conduit-elle à des valeurs caractéristiques À identique car : 

- les valeurs caractéristiques ne changent pas avec le changement m en -m+1-2À 

- la valeur caractéristique À est solution d'une équation transcendantale construite à partir des 
fractions continues définies sur les trois coefficients A, B.C équivalente à la fraction continue 


"rm 2) m 2 
construite par transformation des coefficients dif z Mi (=) de la forme 


SK) f(k) our laauelle : = k SK) = SK) = 3 
HD En EE E 


On peut donc à ce stade écrire : 


xell,+0] ` Hëlt", Et) D diem pegaga i=1,2,3,4 
X 


= (yx 
Fonctions de Bessel f z2 (x) J=J, Y, HO, H® 


La normalisation des fonctions sphéroïdales radiales est liée à leur comportement en lorsque yx- 
>20, Les fonctions hyper-sphéroïdales radiales suivent le comportement asymptotique des fonctions 
de Bessel correspondantes, à savoir : 


Lim HS) |Z Lim La (y x) Sin y à - =) 


Pur HSE (x. vk E Lim SS tx = —(y Pia Cos 7 x SS SR = D) 


En injectant cette normalisation dans le développement, il vient : 


7x0 


1 3 SH 
Sachant que Lim |l-— | =Lim -=| =l 
x u 


Lim HSO” (x,y? )=C, Lim EH ET > dizay) y x E re 1) kr) 
EEN x 


Kam 


Lim HS” (x,y? )=-C, Lim hi- =) Ga Ý därs E - S SE kr) 
i (Se? 


GE HS (x, le C ( EE SE vw = a — D) 5 IR SE (2) 


Kom 


y X — co HSO” (x,7°)= -C,(y ie Cl? x = a D) SCH 1} dr (>°) 


Kam 


+00 Lee zb 
Ee EE l Dar 2) 
k=-0 


k=-0 


Soit les fonctions normalisées suivantes : 


1 j 
ja 
x e [l+] ns” (er) fZ _ a qgazmza (y2 )faa y x) i=1,2,3,4 
2 EE 


(x) 
LE 
Fonctions de Bessel f n342 (x) J=J, Y, HO, H® 


Dans le cas de degré et d'ordre entiers n et m tel que n>=m et À demi-entier, la série n'est pas 
infinie dans les deux directions, elle s'arrête les indices k négatif lorsque k < ko et comme la série est 
construite en changeant m en -m+1-2Àà. kim) est égal à la valeur suivante : 


k(m) LH Donc la limite devient:k (-m+1-2A)= ES Donc la série 
2 2 


commence à l'indice x > jun) res avec 8 =(n-m)mod 2: 


On peut donc écrire la série sous toutes ces formes équivalentes lorsque n>=m et À demi-entier, 
compte tenu de la normalisation : 


E BE 3 Get? 
X _ es 
Se >, de géie (92 ont 
2 Sai dip» ASE (y x) 


1 \2 

j= 

E | Z) < 2,-m4-24( 2) Sasa zl x) 
d ,—m n+Àà+ 

2 Š AV J4,-m+1-24 (2 KS E R D x) 

SEH 1) d (y Je- 2 


HSO” (x, z )= p- m8 


132 

iz 

E | A < DEET dE Sasa aly x) 
d à EEN 

2 e "o J2,-m+1-24 (02 EC e (>) (y xÝ 

KS 1) d, ( ) 


2k=m-n+û 


UE? 
(Ez) 
E SEa eag) £ anon zb x) 
n,k 


2k2m-n (y Sek 


2k>m-n 


i =1,2,3,4 Fonctions de Bessel IL (x) f=J, Y,HW,H0 


La limite lorsque le paramètre y s'annule est la suivante : 


Lim HS X y? = m J a(x) Lim HS” X y? — TT Y, (x) 
y>0 g y ? \ 2 x? 70 n y S y 2 x 


Lim Hs6»r| Z y? |= Ja H® (x) Fän HS x seh ES HO (x) 
70 D y 2 x? 50 Se e S SC 


De l'équation différentielle hyper-sphéroïdale, on sait que le Wronskien de deux solutions 
indépendantes est obligatoirement proportionnelle à : 


1 


KEES 


Eë 


Etant donné la condition de normalisation à l'infini, il vient immédiatement : 


Wronskien faso” (x, 7° ) HSE” (x. y? } SS 


Wronskien faso” (x, y? ) HSY” (x, r? } = 


Fonctions hyper-sphéroïdales angulaires de première et deuxième espèce pour x>1 


Les fonctions associées de Gegenbauer de type 1 et 2 sont définies par les formules suivantes pour 
les deux domaines respectifs de valeurs des paramètres : 


tt 2° ARA G) Cat 2 ’ rm D () 
CH, " Dä) T(v-u+1) Le Dä Ifv-u+l) ul 
HEN LA NE 
3 22 de Pn D Zä Ta aiala] 2 2 de Th +u+2A EX Liu 
Co, Ta) Lu 4 Ce (2) Ci, ()=e | ` | Ta) Le $ o A 
v,ueR e 1-24 E J | 24 E SCH E E 
g C =(-1) 2 l P C = ASE 
Zei Zelt: TU) S 3 T a Ta) 4 Je 
de a Zo 
1-24 2 2 1-24 A+ 2 2 1-24 À--+u 
ES? Cäekk) S SCH 8 1 i ge il ) Cie Ser (2-1) | Qi ) 


Les fonctions associées de Gegenbauer de type 1 et 2 sont connectées par les formules suivantes 
pour les deux domaines respectifs de valeurs des paramètres : 


u 
Auch u eN > CŻ(z)=(-1}" G-A Cat 


H 


KE Ve de Ee 
| + cigle "` eecht Ska ee nf) EH 


e + 
ER À- UE N gel (1 ue 2Cosl(A + Da (z -1) 2 + 


l Bee ñ \I-z 
AvVER ueh —-A2-ueN-Cf,()=(-1 Ch est 
re e Ars? [ai EEE] 
1-21 Cu 
Luch zf pence ECH 
(-1) 2 "2 
vu ER 
2p+1 GË e IA 
JCT s Ch (Z = el = CA Bike _ T — Sin À+ ur GC Z 
2 20), (z) WW: 1,(0) (z) 2Cos((1+ ur) Gr" ( i (z) 
peN u+h--4N 
2p+1 G-1p l 
P A e LI Ju D =Z Au 
v,ueR EE peN ECH is enne 


Formellement les deux fonctions hyper-sphéroïdales de première et deuxième espèce peuvent se 
construire à l'aide des développements suivants : 


x e [+1,1] 


hpsż¥ (x,y 


hqs 2% (x,y? 


ps (y2)= Sa o riala) 
Ke 
2) S a UD ee) 
(2) 
Ke 
2)= F Farr Eriala) 
k=—00 


) = ES (- 1} CSC (7° Le (x) 
fe 


hasi” (x.y 


Avec ces définitions, seules les fonctions de première espèce ont une formule de connexion 
similaires avec celles des fonctions associées de Legendre dont elles sont construites : 


H 
Zuch uen fait. LECH, kal È (= nos" (27°) 
(i-z) (z-1)2 
AveR ueN (2) iv E 
Ka UU: 2 (1-2) 2 i . 

h Au 2 Z h ; 2 -5 À h Au ; 2 
uch > qsż!(z,y’)=e AE Siy (zy e j a in(( + ur) e, (z7?) 
1-24 (77 

AveR ueN Len hgs#"(,°)= (- jee hgs*"(2,7°}+ 2 lZ pps? (EM 
J (1-2) 2 2 2 EEN 
1-24 SE 
v,ueR Ja Pl peN > hps! lz, )= C1) e hei? LA ] 
z-1) 2*2 
w R | jé i je 
ja L > hqs” (z,y pe hgst” \z,y? )- z Eee A+ ub) laps!" {z,7° 
g SCH GI ` Si Zei allt: lee 
peN  u+A--#N 
D 
2p+l l BCEE z Vi-z 
HER Ai N +4--EN > hqs," i D Den LIT hps" 
E 3 pe u+4-5 EN > hgs, Gol Gi Gah? SE (:,7°) 


Avec les notations hps, hqs pour les fonctions de type 1 et HPS et HQS pour celles de type 2, cela 


donne : 


u 
AveR uu EN > HS) pp y) 
z-1)2 
Zuch uch (y ( 1-22 u LA, 
| = |r -1) 2 2 T (1-2) 2 ; 
HOS" (zy? )=e >) E] gt y? )+ — Sin((2 + uye) hpsż”(z,7? 
l-a-peN > HQS; (z7?) e a qs, Gi 2Cos((1 + ur) Ce in(( uy) ps (z7?) 
1-24 y 
1_ E auf ,2) (2 a (z-1)2 2 auf, 2 DÉEN Auf 7 
AveR uch 7 A- u EN > HQS’ (zy }=( 1) 2 mer hqs; (z7 biye DM 
j La 
1-24 (1 , 
v,ueR ne p EN > HPS (z, )= (1) U SE 7 Area" 
z-1)2 2 
ne 7 i e, 
pt S 2 „iur Z714) à 2 T SI , 1 2 
= —— Side FAN EEN oe Text — — Sin((à + uye) Ihps+"(z,y 
2 Gr) VE ui nn Sé 
peN HtA-—eN 
D 
_2p+1 D Ju 2] piut (2-1): Ju 2 z l-z Ju 2 
v,ueR RA peN u+à NUE, DM e We hqs; by hs DM 


Transformation des paramètres u->-u+1-2À , v->-v-2À ormules de connexion z->-z pour les 
onctions hyper-sphéroïdales angulaires de première et deuxième espèce et les fonctions radiales 


Transformation v->-v-2À 


Il a été établit pour les fonctions associées de Gegenbauer de type 1 les formules de transformation 
suivantes : 


vAeR meN ou v,ueR = eN 
(SL (z) = Sal t 2) C,” (z) 
Sitz) 


m Sin((u + Als )Sin((v z Alse kickt Stalin Deeg (z) 


Ce Si 
Saa? Sack 


Étant donné la construction par des fonctions hyper-sphéroïdales par le développement formel en 
fonctions associées de Gegenbauer : 


hps?” (r° SC) ait e ee) 
k=—0 

hqs?” (z, y? ) = > e 1% d? (7° jets (z) 
k=—0 


ze[-11| 


Il vient immédiatement les formules suivantes pour les fonctions hyper-sphéroïdales de type 1 : 


vAER meN ou v,ueR 2-2 


oeh 


Sin(x(v+24), ; 7 
h „m | 
Sin(xv) de k g ) 


OT Sin((u + 2)x )Sin((v +4) z) hapsi” (z, y? jz Sin((u +v)x hqs?” (z, y?) 
Sin((u—v)x) 
La transformation v-> -v-2À retranscrite dans les coefficients du développement donne : 


dr (7°)- di (2) 


hps (z, y’ ) = 


bas. (2,7?) 


Pour les fonctions associées de Gegenbauer de type 2, nous avons : 


2p+1 


VER meN ou vueR A: peN 
ree- TEE cira) 
| Sin(xv) i 


z eaa +v}r)C" (2) Sinfr(u +v C4, (2) 
Sin(r(u-v) 


E (z) F 


Il vient immédiatement les formules suivantes pour les fonctions hyper-sphéroïdales de type 2 : 


vAEeR meN ou v,ueR ZE eN 
ie Au  Sin(r(v+24)) D E 
HPS -aa G 5 E Sin(xv EE ç á ) 
eus) 
2 Sin((A +v yr) HPS?” (z, y? z Sin(x(u +v)) HOS?“ (z, 7°) 


N a) Te 
Hen, La Sin(a(u-v)) 


Pour ce qui est des fonctions hyper-sphéroïdales radiales, d'après leur construction : 


er) a RAE 
>, dr (2): v+A+2k j Zë 1,2,3,4 


SE 1} GE bilt (y x 


k=—œ 


Fonctions de Bessel Saala) f=J,Y, HO, HI 


xe[l,+| HSO (x, z )= T 


Prenons l'une quelconque des fonctions de Bessel d'ordre v+À, cette fonction est solution de 
l'équation différentielle : 


2 d' 


ai lt lees 


dx 


ao Mir: 


Donc la transformation  v->-v-2À laisse invariante l'équation différentielle et conduit à des 
fonctions solutions qui sont combinaisons linéaires des fonctions de Bessel initiales. 


Comme ei (x) = Sin(zv)J, (x)+ Cos(zv Y, (x) >J, (x) = Cos(av)J, (x)- Sin(av)Y, (x) il vient : 


ee (x)= Cos(alv AM. (x)- Sin(r(v AIR. (x) 
Y,a (x)= Sinlshe +4), (x)+ Cos(r(v +4), (x) 


Ces relations sont invariantes par translation d'un indice 2k. De plus la relation sur les coefficients 


; dt, ,(7°)=d{# (72) , An A 
du développement “-"-224Ÿ "VW / ne fait que changer l'ordre de la sommation, il vient 
donc : 


HS: (x, MË Cos(x(v + 2))4s 24 (x, y? )- Sin(z(v + A Ups (x, 7?) 
HSO (x, y? ) = Sin(z(v -+ A))4s 24 (x, Ch Cos(r(v + 2))4s 24 (x, DM 
Heil (x, vk etaa ps Eu (x, 7?) 

Ass (x, y? ) St tade prg Oru (x, 7?) 


Transformation u->-u+1-2À 


Si l'on se restreint aux valeurs demi-entières de À, on sait que cette transformation sur les fonctions 
associées de Gegenbauer donne les formules suivantes pour u et v quelconques : 


Ca ME sua) Zeller 


: F(22+v-— | Q- 21- 
Gott ELA) sina- a e) Z Cosa re Cl 


g TEE z) 2 dl ah + n)a) (z)+ rel? + n)a), H 


Gre dl apen l- Zelt pc") 


Il vient donc immédiatement les formules suivantes pour les fonctions hyper-sphéroïdales de 
première et deuxième espèce de type 1 pour u et v quelconques : 


v,Uu ER = N 

wei "back RE sina e eg )-Ž cos(a- aaae Any ) 

hqs -a(z vk Häll at ae )hqs}"" SS WÉ y SE £ Cos((2 — u)x) hps a) (z y ) 
di T(v+u+1) ? 2 í 


su NO Tv uz1 nlla Au 2 A 
hps; Ga T{v+ +24) Sin((à + u)r) hps? Ga +< Cos( + ur) hqs? “(2,7 


z)h 
e kb 
h aaa dÉ Tv uz1 | 2 2) T 
qs, Z, y REESEN Sin((à + u)a) hgst” (z,y "Collins x) hps}" Ga 


La condition de normalisation de la fonction hyper-sphéroïdale de première espèce implique la 
relation suivante : 


D (aze dir 22+2v T(v+2k+u+24) T(v-u+1) =j 
EAN 37 +2v +k T(v+2k-u+1) T(v+u+22) 


Elle doit rester inchangée par la transformation u->-u1+1-2À. Soit : 


Ca bk Mot TL e 


T(v+2k-u+1) v+u+2À) 


Il s'ensuit la relation suivante de transformation retranscrite dans les coefficients du 
développement : 


T(v+u+24+2k) T(v-u+1) 1 


dr #ri-2à EA 
v,k (y ) T(v—u+1+2k) T(v+u+22) à 


Pour les fonctions associées de Gegenbauer de type 2, nous avons : 


_2p+l 


Zuch uch ou vueR À j 


peN 


dite [(v-u+1) . 2 “lui 
C'"""ch TI vn es RE W Cos(r (A+ u) (2) 


ne SE T(v-u+1) 
E 24 (z) =e | A TRA+vty) Ciy 


Il vient donc immédiatement les formules suivantes pour les fonctions hyper-sphéroïdales de 
première et deuxième espèce de type 2 pour u et v quelconques : 


_2p+1 


AveR ueN ou v,ueR À 5 


eN 


jajal 


Tv DUT 1) uer (z, af Sc Ste (a + u)) HOS” (z, r) 
T 


Lo G 4 ) S T(24+v+u) 


-i| +u- Ty- 1 
nosey me CE TEE osilar’) 


Transformation u->-u+1-2À sur les fonctions radiales 


La transformation u->-u+1-2À laissant invariante l'équation différentielle des ondes hyper- 
sphéroïdales, les solutions HSO (y y>) sont nécessairement combinaisons linéaires de deux 


solutions Hëtzt" De plus les 2 solutions HSC (y y2) doivent avoir exactement le même 
comportement asymptotique en x-> que les deux solutions HSY (x,y? Y tout en ayant 


également les mêmes conditions de régularité (ou d'irrégularité en x=0). Il s'ensuit que ces deux 
types de fonctions sont rigoureusement identiques et par delà les 4 types de fonctions radiales : 


HO (x, y? ) = HSO (x, 7? ) 


Transformation z->-z 


Nous savons que : 


2 Cos((v+u)r)C#" (z)+ cŻt,(-2) 


uA e v "i i IAEA 

eil ZE EE + u) a)i (e) (0 
iji OT Cos((v + u) x)?" (z)- có" (- z) EE SE SS To iji 

Cio (z) a Sin i Ss d EA z) Cos((v + u) D e (z) > Sin((v + u) zk: (z) 


Dès lors, il est facile d'obtenir les relations suivantes sur les fonctions hyper-sphéroïdales pour les 
cas v, u non-entiers et entiers : 


nps?" (-z,7°)= Cos((v + u)rhips" (z, cl Za +u)rhas*" (2,7?) mokar Seck 
D v,ueN moe PR 


hgs*” (-z,7°)=-Cosl(v + u) z has" (2,7) lh —nlslbreickai as") (D gs" (er) 


Calcul du facteur de jointure des fonctions d'onde hyper-sphéroïdales entre les fonctions radiales et 
anqulaires 


Reprenons les constructions formelles des fonctions sphéroïdale radiales et angulaires : 


HPS” (x, y? )= D (— 1} d, (7° Jo Se (x) Fonctions associées de Gegenbauer Oi SE (x) x21 
k=—00 


Os (x, Ch SCI d, La 2 T (x) Fonctions associées de Gegenbauer Ces (x) x21 
k=—%0 


HS )= Z E u+i= n SES d: SC "to rat) 


aaa (y2) GE 1) HS (2°) A4 (2 )= SE 1) dôf (2) 


Soit le développement de la fonctions de Bessel de première espèce et application à la fonction 
sphéroïdale radiale de première espèce : 


J. Gol 5 EA RE E Vr D (Si) DG ` 
| AO 2 (x) at ES TU +v +1+4+2k)1!\ 2 


2 2 


Heli (x r°)= IS He E D BEEN (y>) (y x) — Jr (x E > = (- 1) FOR De (>) (z m 
S | 2 


AA Li ) E v.k (y xÝ SE E Ee (y? VE pen CU +v+1+2+ 2k)1 ! 
Comme  Vi,,,z D D ARE SS SC A 
k=—œ 1=0 k=—œ 1=0 


u 
mer RÉ gg Elart) (ray 
> HS) “(x,y ) Ge rbl e T-i+v+l+A+2Æ)I 


Waat, vi. V7 y’ es ne NN Es) 
= HSP (x,y?) 245 E OR (-x PÈ F m y”T(-1- ESA E 


Soit maintenant les développements de la fonction associée de Gegenbauer Q (type 2) et 
l'application à la fonction angulaire : 


1-24 1-24 
SECH 1-24 LI TA 1-24 1, 
eat? SE H et. E den 
Auch WÉI Wb u+ 1) Vor rha) IW u+ 1) Lit 


TL) Ifv-u+1) DL) I(v-u+1) 


2p+] WÉI 1 av WÉI Aw- 
| 2 "7 SC 1 = RE l 
122 2 yr Ei ie 22 Jr EU + 
PEN CA LU S KÉ (- ` i Do ) ot k r(a) l "H ' oal ) 


ra l,, ial 2-4 En lyy 
en. (ey ET) a g oee O ete) ein 


Il n'est pas utile de prendre en compte le premier jeu de valeurs de paramètres UA quant aux 
fonctions associées de Gegenbauer de deuxième espèce, car il se déduit simplement par la 
relation : 


ouer M Ce Gg H rel MA) gi 


IW TE 1) SN IW +u+ 24) Län 
1-24 

22 r IEN rh +u +24) , VA Lay 

AH e = 
AveR uch Gleis "mi e Tý = +) (z Ip €. AN 

1-2 
2p+1 A 22 Jr ER du 

R A EE "ll GE Q Al ) 


WI 
BEE RE "te (x IS p[er DEER DES | 1 ) 


v+l H H 
2 rfv + d 
2 


2 2 DS 
2 en E T(24+u+v) À 1,4 21,4 2A+u+v 24+u+v+l 1 
de" (ne EEN 


2 42 $ 
2,A+v—— Ee? r(a +Y +1) 2 2 
eg 22+u+v+21 22+u+v+1+27 

5 CE (x) jun $ +) t PA tu v) x= (x ; 1) (x PY d 2 x 2 ) 1 1 

2,(0)v i i i 21 

(4) = (H Aera Gen BH 
+v+22-1 
ire da ve 
Or Je S S 
DA  Cutv +244) (u+v+24+1+2/ 
T 22 +21)= T r 
er | 2 ) | 2 ) 
-Lyu e 

Lg E à d x #(x+1)2 (x— UE Deeg 1 

> Co)» (x) = SH 21 
II EIN 3 I(A+v+1+/) 11(2x) 

Hien a. Fae) a D SL dit (yP )RT(u+v+24+21+2k)1 1 

EOT a Ve RÉI Ef SE Gi ZS I(A+v+1+1+2k) 11(2x) 
Comme  Vt,,,z > KE z Ge SCH D z De plus dues o dÉ aż?) 

k=—œ 1=0 k=—œ 1=0 


Hase (e 2)= pe TP egte- D 1 S D T(u+v+24+2k) 1 q>! ( 2) 
dé EA sf" à TEE E EI 


y d Zeck D g u e 
S- D Fre PURE LE 1 R wa Thu-vt2k) 1a 5 
= HOS a(r’) Ta) 2 GA x D SA dech 
1 
Yre azin) | . u rs 1 1=+00 Zi T(u-v+2k) 1 Dä 3 
xh- (-1) aż! 0?) 


> HOS?! (x, Cl 


r(4) El Mer T(-v+2k+1-4-1)1! "= 


Au 2 2, 2 
Or HOSS (7°) tout comme OS; GZ sont solutions de la même équation différentielle des 
ondes hyper-sphéroïdales sur le même domaine. Toutes deux doivent donc être des combinaisons 


Es . i (1),4,u 2 (2),2,u 2 
linéaires des fonctions radiales HS, (xy jus! Gr ). 


HSAH (x, y?) 


D'autre part la fonction radiale 


ne peut avoir les mêmes puissances de x dans son développement de par sa 


. % p ENT 2 (1),2,u 2 ` . 
construction. Il en résulte que les fonctions HOS? (x,y ) et 5, EI: ) sont nécessairement 


proportionnelles car elles ont les mêmes puissance de x dans la série de Laurent. Cela implique que 
SE 2*(- 1) dé Tia (72) 


le rapport : m y”T(-l- rre est indépendant de la valeur k, autant choisir 
E, wi Th(u-v+2k) 
a Lui 
en T(-v+2k+1-2-—1)11 TI" 


Sa (-1) (ai Sa SI 2) 
donc ce rapport à la valeur k=0, soit : kb Arl- SSES 
Pav) RE (ya; (y?) 
ST(-v+1-2-1)1! 
d'introduire formellement le facteur de jointure. Si l'on définit le facteur de jointure comme partie 
du facteur de proportionnalité des deux fonctions comme suit : 


C'est ce lien qui va permettre 


HS (x,y? )= KOA (y? pue). ivre Ja: 3} gosir, (x, y 2) 


Alors le facteur de jointure s'écrit : 


+ (— dÉ uni zt 2) 

T(2)r(1+v-u) y” EE 
Res Ëss D LS yaro) 

A Tr(-v+1-4-/1)1! 


— AA 


1 
KOAH =2 2e" Sin(A x) 


On a établit les relations de transformations suivantes sur les fonctions HPS et HQS : 


Sin(r(v +21) 


HPS?” Lk 
SE e) Sala) 


HPS (z) 
Zuch uch ou v,ueR is peN- 


posis, (SP ES OS) 


Sinlr(u A) 


Pour des valeurs entières d'ordre u , on sait que : 


u=meZ— HOS*",,(x,7°)=HOS/"(z)-x e EE HPS?” (z) 
Sin(v x) 


Dans ce cas la liaison entre les deux fonctions s'écrit donc : 


Seck, Set nee) (à Puos- P 0) rase) 


er kim (2°) 


PE ee Ga Sin(2 x) 


| Sin(v z) Jr mosie (e)s sia cker 


T 


Si maintenant le degré v=n est également entier, alors : 
v=nmenN HSO” (x,y? )= KO (y 2)HPS?" (z) 


Il est possible de montrer la relation suivante sur le facteur de jointure introduit ici : 


e Dais) 
T(2)T(v+u+24) y"  ÆT(v+4+1-/)/! 
Je At) SE 1J d7; ui (2) 
ST(-v+1-2-0)1! 
gear) 


piiks DE giva g tan n(A ei TU ID 4 nl y" ST(-v+1-2-1)1! 
—v-2 


1 
KO HH-24 — 3 2 e7 Sin(A z) 


v—-2À 


Tv +4+1-1)1! 
e k Gi )= dot (7°) 


A? pii- zt AP Ha (y2) 


2 
x e tir z| sinla x) r(4) | 
D 


Sin(z(v + u +24 -—1) 


E Ktiu Gi eet GL: (ËM GL: Le Gi ) S 
En appliquant la transformation v->-v-2À sur la formule initiale de liaison : 


pt) site Ne verre ST oi (ey) 


Sin((u +v +24 -—1)r) EE (2- 
m Sai x) 


= HSE Gr") KO) HF rosz (x,7) 


En usant de la formule de connexion entre les deux facteurs de jointure : 


2 
—i(22-1)r 224-1) Sin( a Cl 
Ha | in( x) y? E 


Sin(r(v + uz 24-21 


KOZH Por LG" Wi (y2 J422 (°)= 


Il vient : 


KA Nr 


E 
SEH bi ja Sin(x(v + u +24 -—1)) EU: FERAN EE 


mT e Jaa Sat? cl WÉI ) 
PURES (7°) 
2 


IL 


SE D e~) 22 Sin(A x 


KOA ara p? E 0>) HOS” (x, y?) 


=> HS (x,y? )= 


On peut encore tourner cette formule de la manière suivante : 
Hs) zu (x, y? ) — HS OAA (x, y?) 


Ge fau Tv u+1 
HQSE MA Ee e ea) A t hosi” (re) 


(1) A,- u+1-24 2 Tr(a) | CF tes 27%- Sin(à x) À,-u+1-24 2 
HS," (xy BESCH KO pp a eap arrg) He S (xy ) 


(all 
Tr (A) S cb 2) ei Ve) 2°% Sin(à x) T(v -u+ 1) 
yida) KO jay agy ayro + u +24) 


En comparant deux résultats obtenus, j'obtiens une autre formule très utile : 


Se HS UAA Gei) HQS?” (x,y?) 


ër gp :) 1 Tr(a) 2 ee 22 Sin(A x) zu 2) 
LES xy E KOZ) yida ape ae) HOS” (x,y 
—i| A—— |m 
Dir Wa, A 1 T(v-u+1) { Tr(a); e Lä eV 221 Sin(1 x) PE 
HS (x, 7° )= HS; ze (x X,Y ) EH EMA 4 (>) HOS," Gr ) 


=> RUE (7° ) = Kä (7° Vera 


y 
T(v-—u+1) SS 


Une des formules précédentes introduit donc un deuxième facteur de jointure : 


Gout d e HW 2241 Sin(A x) T(4) i e HW 2247 Sin(A x) TL T(v-u+1) 
La VJS KO AAA ANA REA) y’ Ndr T KO (y2 az (y AA?) ydr T(v+u+22) 


—i a Tr: 
=> HS» (x, Cl —e z) eY 7 KON (°)4oS°" (x, y?) 
Toujours dans l'hypothèse où v=n est entier et À demi-entier, la formule de connexion suivante est 
applicable : 


v= a Se E HS A (x,y JE 1) Hs "te, y 2) HSO (x x,y Ge Si 1)? HSO (x 7°) 
= 2 HS QE re SS SE 2)= (1) HS7 (x, y? )=-- 1)? HSA (x, y 2) 


Alors si u=m est un entier ou même un réel, n entier et À demi-entier, il vient : 


HS (x,y? )= (1) HSE (x, 7°) 


—i(u+1)r 424-1 Ç; À À S = 
=> HSO” (x, y? ) E emt Gs GC, CG | SH A HE HOS?" (x, y?) 


On a donc trouvé un autre facteur de jointure liant les fonctions radiales et angulaires de deuxième 
espèce. On les note : 


a (-1) 1) Lët SS ch 2) 


-2v- See Tr(a) +v- À TL 1+4)/! 
usala etle asiha) Kele) ea? etes BE EE RTE 
S v+l-2-1)1! 
, , | er 224 Sin(A z) ra) Ÿ gue 
HSO (x, y? )= KOM (7° Juge": Ga Soot DE T(n+u+24) 


Ces facteurs de jointure peuvent également se calculer d'une autre manière. 


Calcul du facteur de jointure entre les fonctions radiales et les fonctions anqulaires dans le cas des 
ordres et degrés entier et de À demi-entier 


On a donc établit pour ces valeurs de paramètres les relations suivantes : 


HSO (x y )= KO (2 )HPS4" (x, 72) 
HS 0)" (x, y? ) GG Kä Li Jusen (x, y? ) 


Dauf 2) (1) 27 Sin(A x) r(4) ` T(n-u+1) 
Wë Siläcl? 424 (y? Li Hi SW di SH SE 


Pour obtenir une autre formule de calcul du facteur de jointure, il suffit d'appliquer la liaison pour 
la valeur d'argument x=0 tantôt sur les fonctions, tantôt sur les dérivées premières. 


Dans le cas des paramètres u=m, v=n entiers et À demi-entier, on a établit que les développements 
étaient unilatéraux vers les indices k positif, commençant par un indice négatif -kO : 


0 -|277| hps” (x,y°)= 2er} 


Li 
HS 0 (x, y? )= EL z) 5 diz VAE af, +) Tegak à 


k=-k6 


AA dÉ SCH 1) di -m+1- SÉ 2) A7" (y dÉ SE 1) aż (y 2) 


(sa? k=—ko 


Mais avant cela il convient de calculer les limites en x=0 des fonctions hyper-sphéroïdales radiales: 


i 


e ET 
Si n-m par —— 


sl} ha (x ui EE à —— a : J Sd”: ,-m+1-2 Char x) k = Es SE z SÉ 1} girafe) k = 2 
È 


k=-k , e hn-m-l 
1 ° Si n-m impair 


-E Lech REES Ze elt ar 
2 


DEED 
PT a nom e Baak ah: 


k=0 


Sr p - 1 yx m+À 
Si n-m pair ue a ER WI Lu il x) HE 


l 
Lim HS>” (x,y? )= Pa) 
AA n ( dë ) 2 J A TEE, y? i m 1 ER 


Si n-m impair Lim ei d SC à hat x) Jui ( x) E T(m + À + 2) 


À,-mH-24[,2 
m d 7 | 


Si n-m pair Lim B x- Ip ——— 
Li Seil dÉ T 1 #0 (2 TEEN 
St ý ER 9 WË Kee 2 À,-m+#-24 
LZ Y y A+m+l S dm 
Si n-m impair Lim x -1 —— 
í E Si Ta +2) 
MT d” -m+l- up) Si n-m pair 0 


Jı Je? K SS Dm 2 À,-m+1-22 
Ge SC —m+l- ay?) T(m Not 1) E s dx b : . Jaye? See al 
POLE EE Li T(m ++ 2) 2 


> Lim HS (x r)- Si n-m pair 
x>0 


Si n-m impair 
Si n-m impair 0 


Appliquons maintenant la proportionnalité en x=0 sur les fonctions ou leurs dérivées premières : 


St nm par Echt? e GE e HPS"(0,y°)=e ? ips"(0,7) 
,À,m 
Hse)" aps ty) E ds) _ 7 desler) 
oem api Kilochit- 2 we 4 "7 SES 
a Lë, x=0 e 
dz E 
, de 2 
Si n-m pair ek Si À 
Ai 2) Tfm+4+1) 422) pps" (0,7?) 
>K” ()e e EE 
Si n-m  impai — g 
i n-m impair a0 TEE 
Z z=0 


Facteur de normalisation radial des fonctions d'onde hyper-sphéroïidales 


Dans le cas des degrés et ordre entier et À demi-entier nous avons introduit le facteur de 


normalisation de la fonction radiale, sous la forme d'une somme unilatérale de coefficient du 
développement : 


A> Smal (y dÉ SL d"(y Ze SCT (bai: TE e EE 1} Ai em WÉI 2) A =(n—m)mod 2 


2k>m-n EE E FS ô 
Act: SU" S (1) "(y dÉ Ÿ (1) d/"(2) a =(n+m+22-—1)mod 2 =(n—m)mod 2 
2k>-m+1-22-n k= [m d BEE? 


Dans le cas des degré et ordre non-entier, alors les sommations sur l'index sont bilatérales, et le 
facteur de normalisation radiale s'écrit : 


aa) Éd Ge) rk É CN dirt) 


Wronskien des solutions radiales et angulaires 


Tout d'abord d'après la forme de l'équation différentielle des ondes sphéroïdales, tout Wronskien 
de deux solutions indépendantes de l'équation différentielle est proportionnelle à une puissance de 
1/(1-z°). Ensuite, à l'aide des formules asymptotiques en z pour les fonctions radiales, il est facile de 
démontrer que le Wronskien des fonctions radiales est : 


1 
22+1 


Wronskien {Hs 024 (x, y’ ) SE (x, y? } = 
y% (e S 1) 5 


A l'aide des facteurs de jointure introduits précédemment dans les deux formules suivantes : 


; fish 
HS9 (x,y? sik EE E 2) HOSZ (x7?) 
HSO (e, cl -e ere an(s (a,r?) 


On peut calculer les Wronskiens pour toutes valeurs du degré et de l'ordre : 


Wronskien {45024 (x, y? } SERRE (x, 7? } 8 


Heli (x, y? ) = KUn (7° per ue CEF gosn, (x, Dal 
m Sin À m 


2 
T(v-u+1) 
T(v+u+2À) 


KO (y) KOA elek 27317 Sin(A x) r(a) 2 Hessi 


Rs DE et 2221 Sin(à x) | r(1) 


KO (yae AAA EE) WRC 


í DE EH WEE 


HSA” (x, y’ ) — AS (r)a gO) Li 405" (x, y’ ) 
aus} ; A,m 2 À qe 2 
a. 2) Sin((2 +v Jr) HPS (x,y )- Sin(z(u +v)) HOS°" (x,y ) 
Hosa (er) Sala ui 


Hell (x, y? ) GE Cos(x(v + 2))4s 24 (x, y? )- Sin(z(v + A Wiselen (x, DM 
Sin(x(v + )) 
22+1 
re 1) 3 
T(v+u+24) PU Gel d j 7 e7 Ain Li e (2 )sin(r(v +1) 
22+1 
TG-u+1) Sin((u -vyr x? SES 


Thv +u+ 22) AC} SS Ae ] e7 A2" (° ler (2) 


Te 


=> Wronskien {Hs 02 (x7?) HSE (x, z )}= 


—v-2 


Wronskien {HOS (x, y’ } HQS” (x, y? } = 


=> Wronskien {Hps? (x, 7} HQS?” (x, z )}= — TG SE 1) EEES] 


L -1) 2 


Finalement le Wronskien des fonctions angulaires hps et hqs, à l'aide des formules de connexion : 


AvEeR uch ou v,ueR Josi 


Hpëicl 7 Ce 
(x-1 


re hqs? “yh Cxhps*"(r, 7°) 


ZP 
2 


HOS "(7° }= AE S 


Il vient : 


Wronskien {HPS?" (x, y 2) HQS" (x, D 2} 


[G+u+22) l =} EICH BE "eil Ain DË He SÉ 2) 
T(v-u+1) CG) e- j 


2 a. 
Pronstien psi (er prasie lor Mie ëch? 


e 


Solutions exactes de l'équation différentielle des ondes hyper-sphéroïdales pour u+1=3/2 et u+À=- 


1/2 


Modifions l'équation différentielle de la manière suivante, et trouvons deux solutions exactes de 
l'équation hyper-sphéroïdale : 


1- 


(1 ll ortt orh x?) dar: D), (x)= 0 


est eye) (a at (ur së SIE 


1— x? 


22 +1 
B=- 


dx? 1— x’ 


e u(x) , sis dura) h IRE) 


EK 


A 
Nr 
2 2 

ufu+27)-(1+7)- 
A(+2)-- 

e {u+22)-(1+2)+2(1+2)=0 8 au +22 )+ 2 -1=0 6 (7 +u +13 +u-1)=0 


= (AO lesch ele 


Simplification ulu + St D + DI + x27)= (1 Ge Cl | 


2 


d d'u) 


Simplification © =a = 2 + 7)= +y ?u(x)= (0) 


TE E E 


UI 22+1 
=> y(x)= (1 -x7 ya (A Sin(y x)+ B Cos(y x)) 
~ 3 1 2 
G A+u-l=0> u S EE EE Z À 


22+1 
=> y(x)= (1 — x? Ja (A Sin(y x)+ B Cos(y x)) 
On a vu que ces solutions correspondent effectivement aux cas où la valeur propre est déterminée 


o=(5-1)5+4) 
par la résolution d'un système matriciel réduit à un scalaire : 2 2 ; 


DI u=-1>0=a=0=> y(x) A Sin(y x)+ B Cos(y x) 


2 


On retrouve les solutions pour le cas À=1/2 : | SR | 
(2) u=1=w=a=0— y(x) A Sin(y x)+ B Cos(y x) 
LA? 
DI u ‘ > ©=@ É y(x) A Sin(y x)+ BE x) 
Pour À=1, il vient : | D j 
(2) u 1 Hi 3 SE A Sin(y x)+ B Cost x) 
2 4 nf 


(1) A 2>0=4 2 y(x) A Sin(y x)+ B Cos(y x) 
Pour À=3/2, il vient : | S SS Sp 
A Sin(y x)+ B Cos(y x 
2 0 _2 
( ) A © CO y(x) a — 


bras 5 15 G- A Sin(y x)+ B osl x) 


2 4 (-x} 
15 jaje A Sin(y x)+ B Cos(y x) 


Pour À=2, il vient: 


(2) u=0 panei uf 


(1) u Jsou 6 Dë À Sin(y x)+ B Cos(y x) 


Pour À=5/2, il vient : bës x} 
(2) u=-1=w=-a=-6= y(x)- A Sin(y x)+ B Cos(y x) 


(1 — x? E 


Comme dans le cas des fonctions sphéroiïdales, il n'est pas facile de construire les fonctions d'ondes 
sphéroïdales liés aux valeurs propres w solution d'un système matriciel fini de dimension 
supérieure. 


Les solutions exactes pour u+À -Ż 


Selon les propriétés de parité des fonctions hps et hqs, de finitude sur l'intervalle de [-1,+1], et 
notamment en x=+1 et x=-1, cela impose que Si"&7)=0 ou Cos(y)=0, On a donc : 


2 zc I(n+u+22) 
24+2n(r(1)) T(n-u+1) 


Comme Jasos" x, =2 
R 


e à DES 2 e (EES 277 5 S 
2=2 > fax hos D =2 = we KAS 
u+ a d ps? (er ) patan CA dn en) 0 fe S 
2 


Les fonctions HPS et HQS peuvent se définir à partir des formules de connexion établies pour la 
valeur 1+1=3/2 quelque soit l'argument z complexe en dehors des coupures respectives des 
fonctions de type 1 et 2. Selon la construction des fonctions hyper-sphéroïdales, cela implique que 


HI 
BIE jat peN>ueN 
À soit un demi-entier et u un entier : 2 2 . 


5 a Sin Pa à a Cos? Pa 
nr 1 EIER IEN 2 1 2(2n+21WT(2)) 2 
GTa SE SE à SE? Ÿ 272 7(2n+22+12 Seel di il 
lo;) E wl 2770» n 3 (q-e)? 
ER ` 1 oœ: (-1¥} Sin(y,x)= @;Sin(y,(x-1)) n pair <-n-1 impair 
hps, (x. Y )= z RA n+l 
(1-x°) + Ier t- UU 2 Cos(y,x)=o$Sin(y,(x-1)) n impair n-1 pair 
n+2 
At @(-1) 2 Cos(y,x)=-@:Cos(y,(x—-1)) n pair <-n-1 impair 
hgs, (x Y )= Zr nAi n+l 
(1-x°) 4 Ier tU? Sin(y,x)=-w{Cos(y,(x-1)) n impair <n-1 pair 


Les formules de connexions se présentent comme suit : 


ds in KA de =. aA 
ver pen QE re) mos o ge EE pp 
(2-1): GA 2 ÿz-1 


Mais ces formules sont valables pour une construction à l'aide de fonctions associées de 
Gegenbauer et aboutissent à des fonctions à valeur complexe. Or ici la construction des solutions 
exactes est différente et surtout directe, on peut donc proposer une formule de connexion 
simplifiée en éliminant la partie purement imaginaire : 


3 # 3 3 g 
a" E nos EE 


GÉIE 


Les fonctions HPS et HQS pour x>1 sont donc les suivantes : 


Hé 1 ; 
HPS,” (x,y? )= el ege ee 
ee LG: = 1) 3 (O; Sin(y, (x E 1)) n impair 
Fe hs Gs 
EE 1 œo\Cos(y,(x-1 n pair 
HOS, © kache) 2 Ee 
Li L 1) ` Le Cos(y,(x— 1)) n impair 


Pour ce qui est des fonctions radiales, la condition de normalisation à l'infini est la suivante : 


3 
Lim HS? (x,y? )= (y x) SE SS 77 D) 
y x> 


3 
Lim HSZ (x,y? )= (7x) #3 Col: x ES SS S gd 
EE 


Or les fonctions sinus et cosinus de la condition de normalisation sont des combinaisons linéaires 
des fonctions solutions obtenues plus haut. Les fonctions radiales sont donc des fonctions 
comportant un déphasage et respectant la normalisation à l'infini : 


E (x-1)- En) 


EES) 2A+1 
nm (x ES 1) 4 Yn ? 
C (24-1) 
3 ol e 7240) 
rä 2) 1 4 
HS >Yn = EES 24H 


LG: —1) 4 Ya ? 


Comme À est demi-entier, il vient : 


Sin. (x-1)- pz) 


(22 1 
LA HS, F2 (x, MÉ 2441 241 

ie, (x? —1) + Yu 2 

2p+1 D pr 
ER uo. PA 

- SE ` i Cor, (x 1) 2 ) 
HS, X, Yn J57 24+1 EES 
Li Sat) 4 Ka S 


Le Wronskien de ces deux dernières solutions est donc le suivant : 


OT CT re 


Cela correspond à la valeur usuelle du Wronskien des fonctions radiales. 


Les deux constantes de normalisation sont liées aux intégrales suivantes : 


e i kam o PAF) a (F) 
hag QU = l ni m >l, = [a 27 
Le STE EE =l hf -1 D 2 
Col TT; Cos?| DZ x 
1 ar) H KE ) V | 
aie n=2p+1 (n+ H ) [ax — 1,=[& i 
2 
2 


Parmi les intégrales à calculer commençons par les suivantes déjà calculées avec les solutions 
exacte de l'équation des ondes sphéroïdales : 


n=2p peN 2y=2pr 
=x e 1+x 


fas Sin ie eh, (Ja GE x) lan 2) Sa nm 
0 


+1 27 
t=2y(1-x) faroera) falt) 


0 


+1 4y 2y 
reais fal Ur) fat cost) LEES 


0 


+1 D dy e: +1 Sin Kä +00 
> fax a x) _ l fa! Col >n pair jé 2 2-7 y Euler + Log(2nx)+ Í a oeh) 
=x À t 2 t 


Zus 


n impair 


À : Cos? 2 x) 
De même on obtient : jl | 
] t [dx — s ee epech [ dt Coste Ea 


Zu 


Puis les suivantes : 


ja Sin’ (y x) ` [a 1-Cos(2yx)_1 (Ja 1-Cos(2y x) f 1-Cos(2y x) a [a 1-Cos(2y x) A [a 1-Cos(2y x) 


(BEA l-2} A Daf TT ep, à I-x o} l+x 
We n=2p peN 2y=2pr 
+1 2y 
äu PACE fat [a 1-Cos(2y x) -2 Ji: Ge 
0 We? 0 (x) 
d 1-Cos(2y x) 1- Cost) l- Costy l- Cos(2y x) ` H 1- Cost) H 1- Cost) 
t=2y(1+x) ja en -fa ja [a WE 2y [a S [ai S 


Jär, TI ski? Cos(t) Au Gel 


Par intégration par partie fa 
0 


cl a cirai. É [et + r =o] | Í E SC l Di À sit) | Ï en) 


3 -ef 4 4y d 4 


Z 1-Cos(t) "€. Cos(t) H Sak 1 a ) 
r fja =y Euler + Log(4y)+ [a >n pair [i= S A Ever + Laon +27 Î ma Y dE dt d 
0 4y -1 


dës 


1- Cos(t) l Gelb Ja Sin(t) 5 Cos(x)-1 + [at Sin(t) 
‘à t ve 


HIT 


3 ; Cosi 
De même on obtient : a | ) 
n impair [a—— —- ` y Euler + Log(2nx) +27 Í dt —= Sin(e ), fa dt CH 


Si D — x? f 


Où l'on voit introduites les deux fonctions connues Sinus intégrale et Cosinus intégrale : 


Zu 


Cosint(x ja et Cos(e ) Sinint(x) = [at Se? 
0 


3 
Solution exactes Dour, 1 = _ 1 provenant de la transformation |” +4 = 5 
S u—-u+1-22 


Avec les propriétés de transformation du paramètre u en -u+1-2À, il vient : 


v ueR EI peN A+u=> 
1 
E win 
DS, Zjem FR APS, ZG PS, Z 
ER (2v +24 +1{2v +24 -1) 
1 
pa") E (a) : vw" 
pq, "kel hpu" Let ps,” Le 
dun (2v +24 +1{2v +24 —1) 


Il vient immédiatement : 


. 2 nr 
nz 1  (rG})(2n+24-1on+24N2n+224) 4, | 2 d 
teacher 8 Ja 
el a fa)? 


1 (TA) (Qn+24-1)2n+24{2n+24+1) Î k Co") 


Wl Ze , Laf 
te y 1 &Sin(y,(x-1)) n pair <n-1 impair 
hps, ? (x, Y )= 7 2H | eco: ` S S 
(i S DES &Sin(y, (x- 1) n impair<n-1l pair 


SE 2) 1 o 1)) n pairen-1 impair 
DN 
) 


HES 


Laf" 


Les fonctions HPS et HQS sont les suivantes : 


CC 


O, Cos(y, (x 1 


n impairr<n-l pair 


1 


aee a a re por 
Po (x —1) 3 (O, Sin(y, (x-1)) n impair 

= > 

1 ©:Cos(y,(x—1)) n pair 
2 |a Cos(y,(x— 1)) n impair 


e=" 


nos" "kack 


Les fonctions radiales sont les suivantes et respectent la normalisation à l'infini (elles sont 
identiques à celle pour le cas 1+À=3/2) : 


EE (2 2) 
See E EA? 1 4 
HS, X Yn 2441 EES 
nm (x? = 1) 4 Fr. 2 
MD (22 -1) 
ete 1 4 
ZE BK ee 224 224 
(x? —1) + Ya? 


Comme À est demi-entier, il vient : 


j um. ZZ 
Dun ( d 1 Sin. (x 1) > ) 
= i zl, J= 2441 2441 
Van 2 (x — 1) 4 D S 
Saai p 
SE SC Gu -t-a i Cosy, (=-1)- 22) 
HS NM (x, Ya ) et 


22+1 22+1 
2 ER En 
(x gei 1) a Va S 


Construction des solutions exactes de l'équation différentielle des ondes hyper-sphéroïdales lorsque 
u+-1/2 est un entier supérieur à 1 


On peut calquer l'approche de l'ouvrage de J.Meixner & FEWSchafke « Sphaeroid_funktionen » 
dédié aux fonctions sphéroïdales n°3.534 page 274, il est écrit que les solutions sont de la forme : 


Poe ENS SC) Li Die ap (x) pour x> 1. Où P est un polynôme de degré au plus m-1. 


Fort de ce résultat, on peut proposer l'« anzatz » suivant à intégrer dans l'équation différentielle 
des ondes hyper-sphéroïdales lorsque x£[-1,1] : 


D list a: e va) 0+7 x?) OCT Mai 


1— x? 


Pour meN? u+=m +. y(x) =o AT e Pr (x) = (1 -xJ d: EEP A (x) 
1 


1 2A-1+2m 
—x2) 2 
P Lellig, x"? +... +ax+a 

m-—1 R m—2 se 1 0 


en fixant à chaque m, et déterminer les conditions de l'annulation des divers termes d'un polynôme 
résultant afin de connaître les termes dudit polynôme. Il vient pour m=2 : 


24+3 SES? 
v()= (1-22) + erfarit te] tq 160° +80(44 —5)+164* -402 +647? +9=0 
y 


On retrouve les deux racines du polynôme déjà établi auparavant : 


5-42 —-41-47° 
Oh = 


4 
160° +8@(44 —5)+164* — 402? +647? +9=0 => 
Ge 5-42 +4,/1-47° 
2 
4 


Et les solutions réelles sont à prendre parmi les parties réelle ou imaginaire de la fonction y(x) ainsi 
construite par injection des seules racines réelles. Il vient pour m=3 : 


(x) en 2 j4 -1+40@ _ -9-1287* +40 À -164 +40 0-32 À œ-160° 
y 8y 1287? 


225—1036 À +560 4 — 64 2° -256y° (17 +4 4?)-4(259 +2567? -280 2 +48 2° )@ + (560-192 2) ©° -64% =0 


Il vient pour m=4 : 


pyt l+4o -9 192y° +404 -164 +400 -32X 0-160 
eea ET RS 1287? 
»x)=(1-x) WR 4 
i 225+8128 y -1036 4° -1792 7°% +560 4 —64 2° -1036 @-17927°0 +1120 40-1924 © +560 0° -192 4 © -640° 


30727? 
-11025-368647 +51664 4? -31584 2 +5376 45 -256 7° -128 7? (3149-10962 +801*)- 
16(-3229+ 39484? -10084* + 64 4° +647? (-137+202))o 
32 (087 + 320y° -5044 +484 )o° + 
+256 (1-4 Zu 
-2560*=0 


La relation de récurrence sur l'« anzatz » : St en est à 4 termes. Les 
pG)=f-x) + e Pala) 


valeurs propres À racines du polynôme s'obtiennent en remplaçant dans l'équation différentielle 
hyper-sphéroïdale un polynôme de la forme : a +a,x+a,x? Pour ne retenir que le terme de degré 


0 qui doit s'annuler, soit la condition sur les trois coefficients : 
la? -14m0 -m))+2 (iy a +a,)=0 


Ce qui constitue un polynôme de degré m en À. En prenant le motif suivant intégré dans l'équation 
différentielle des ondes hyper-sphéroïdales : 


H 141 1+2 1+3 1+4 1+5 
4X +a X + di2 X + d;,3X + d;,4X +d,,5X 


Puis en recherchant le terme de puissance de x comportant 4 termes successifs, il y en a trois mais 
qui se déduisent l'un de l'autre par translation d'indice puis en annulant l'un des trois termes, on 
obtient la relation de récurrence suivante : 


.40+42 -1 
=1 a AE 
8y 


-2iy(l+1-m)a, Je en +3m-—m° +1(2m-3)-X% +o) an +(1+2)2iy ap +(1+3)a,3)=0 

Ou = do Le) polynôme de degré m- 1 en © 

da =a Le) polynôme de degré m —2 en o 

=> afore EE ya +a,)=0 polynôme de degré m en © 
A2 = Amp Le) polynôme de degré 1 en o 


m—2 


a, 1 =1l 


Pour l'exemple les premiers termes du développement s'écrivent : 


404421 2 [m1 Wer äi Uer Ai 0) ms 
Fer A 2 1287? S 
STEET EEGENEN 
 LJäe 437 -140 +42 —äläe +44 -25 Sp 
e? K ) r 
30727° 
y(x)= D CEET ei" In "lan 3 Mm -3) + 
l, =287+408(8 +0)-112(2 +o} -16m (40 +44? —13)+3m4o+42 äere? -9) | ns 
7687” 
, o+42 äer: A0 -94o +42 -25f40 +47 -49) 
983047 * 


Solutions exactes régulières en x=1 et x=-1 


Dans cette section on va supposer que le paramètre À demi-entier est supérieur ou égal à %. 


Pour m=1 : la solution s'écrit : 


deet) nid es vii: E EE 
l—=x°) 4 1—x2) 4 


o / P(ow,2,7°)=4œ©+4%# -1=0 


La régularité des fonctions y1 et y2 est obtenue lorsque qu'elles sont finies toutes deux en x=1. Or 
ces deux fonctions présentent un développement suivant autour de x=1 : 


»(1)=0 2 Cos(y)=0= y(x) ~ = Sint) 
2 4 (=x) + 

ll en Sin(y)=0= y(x) = = Cosi) 
2 4 (=x) 


Il suffit donc que y soit solution de l'une des équations transcendantales pour l'une des fonctions y1 


ou y2 EE que w soit racine du polynôme P1 et que À<=3/2 pour que y1 et y2 soient 
Sin(y )= 0 


régulières (elles s'annulent toutes deux en x=1). Il y a une infinité de solutions de ces équations 
transcendantales, indicé par n=1,2,.... 


Pour m=2 : la solution s'écrit : 


1 E E E | 
Ge; Lette) (ei nla) 
D 4 Sp 


a 8y x Cos(y x)-(42 -1+40) Sin(y x) 


2243 y x)= =y, (x) 
8y (1-x°) + 
e? (x) a (42 —1+ Ae) Cos(y x)+8y xX Sin(y x) Ta à Si 6 
8y (=x) 4 


en / P(w,2,7°)=16&° +8o(44 —5)+164* -402 +647? +9=0 


La régularité des fonctions y1 et y2 est obtenue lorsque qu'elles sont finies toutes deux en x=1. Or 
ces deux fonctions présentent un développement suivant autour de x=1 avec h(x) une fonction 


régulière en x=1, mais ne s'annulant pas en x=1 : 


2 E 
D 0 8 Cos(y)=(42 -1+40)Sin(y)= »(x)= A le Sin(y)+(1-x) + Sin(y}a(x) 
1287 2 4 (1x) 4 
P,(©, À,7° 5-24 
v,(1)=0— (44 -1+40@)Cos(y)= -8y Sin(y)= y, (x)= ke, 4 le Cos(y)-(1-x) + Cos(y)h(x) 


12872 4 (=x) 4 
Il suffit donc que y soit solution de l'une des équations transcendantales pour l'une des fonctions y1 


= 2 — 1 D D A 

ou y2 : JSY Cos(y)= (42 -1 +40) Sin(y) que w soit racine du polynôme P2 et que À<=5/2 pour que 
(42 -1+4w)Cos(y)=-87 Sin(y) 

y1 et y2 soient régulières (elles s'annulent toutes deux en x=1). Il y a une infinité de solutions de ces 

équations transcendantales, indicé par n=1,2,.... 


Pour m=3 : 


sa-l- ) + er 3y E 


o 1 P(o,2,7°)=225-1036 4 +560 4° —64 4° -2567° (-17+4 4°)-4(259+256y° - 280 2 +48 2 leit 190 2) @ -640° =0 
zsl-nlskbindsl ntx)=nG) ` kal 

»(1)=0= 167 Sin(y)= -40 +42 -9)Cos(y) 

Es 0 = 167 Cos(y)= Met +42 -9)Sin(y) 


Et les développements autour de x=1 sont les suivants avec h1(x) et h2(x) deux fonctions régulières 
en x=1 : 


us 2 A4 -1+4@  —9-128y° +40 4 —-164 +40 © -32 À @-16 0° 
X +1 X + 


Alle Bear hA edel DZ cos) 


22+1 
1—x) 4 


P,(©,2,7° a(x) <. 7—24 2 
yga- enr MO nj À A 
(1 — x)= 
Là encore les fonctions y1 et y2 sont bien régulières dans la mesure où y est solution de l'une des 
équations transcendantales et w soit racine du polynôme P3 et que À<=7/2 . 


Pour m=4 : 
pat E 1927? +404 -164 +400-327 @-160° 
a pn ` Eis | 
S Laf , 225+8128y -1036 4° -17927 À +560" 642 -10360-1792 7w +1120 0-1924 ©+560 © -192 2 © -640° 


3072y° 
o | P(0,2,7°)= 
-11025-36864 y* + 5166427 -31584 4" + 5376 4 — 256 4° —128 y° (3149-10964? +80 1°) 
16(-3229+ 3048 42 —1008 4 + 64 2° + 64° (-137+204))o 
32 (9087+320? - 504 4 +484 )o* + 
-256(21-4 Ale 
-2560*=0 
=) () nex- ` rä 
-225-7744 y? +1036? +2567? Z SO + 647 + 10360 +256 °w -120w +1920- S600? +1920? +640" ofp) 
24y(9+647° -40% +164 -400+32 Lo +160) 
- 225- 7744y” +10364 + 2567747 — 5604" + 64° + 10360 + 2567 "0-112020 +1924 © - 5600° +1924/©° + 640° Co) 
DEES NEE EEN 


»()=0= Cos(y)= 


»()=0= Sin(y)= 


Et les développements autour de x=1 sont les suivants avec h1(x) et h2(x) deux fonctions régulières 
en x=1 : 


P,(©,2,7° }h . 9-24 
ne CALME ant. sinyal) 
(1 — xy a 
P,(©,2,7° h 9-24 
aeea ebe A AË Schéi 
1x) 4 
Là encore les fonctions y1 et y2 sont bien régulières dans la mesure où y est solution de l'une des 
équations transcendantales et w soit racine du polynôme P4 et que A<=9/2. 


On peut donc proposer un résultat général qui assure la régularité des solutions en x=1 et -1, mais 
qui reste à démontrer : 


iy x l=m-1 


"Lal LE Dax naire nial- Bella zs m) 
l-x | 4 


_;40+4% -1 
87y 


m-2 


ENEE EN a, +(1+2\2iy a +(1+3)a,3)=0 


a =a\æ) polynôme de degré m- l en œ 1 

| d ) Pa i e P,(v.2,7 =a or Ai ml ml +2(iy a +a,)=0 
ET? lo) polynôme de degré m- 2 en © 4 
— | Polynôme caractéristique de degré m en © 


=, eil polynôme de degré let o | æ racine de P, (o,2,7°)= 0 


a n-2 = m-2 
a, 1 =l 

À, À, m- m 
oca de E). aaa eA ET 


Af Laf? 
2 m+l- 
AT D je BAACH. vr? 


(1 = x) 4 
i j Mod(m,2) 1-Mod(m,2)| Pn OR (x) 2m+1-24 
y tel que u(x =1,4,7 )= 0> rx) (Sin(y)) (Cos(y)) pe "d r seid 
1-x) 4 


gil) 


h (x), & (x), h, (x), D (x) fonctions régulières en x=1 
Jm zl 


=> Lim y(x)=0 Lim y,(x)=0 Lorsque 4< 


xl 


Il existe donc bien des fonctions angulaires déterminées soit par l'annulation en x=0 soit par 
l'annulation de la dérivée première en x=0 et l'annulation en x=1, pour lesquelles w est racine du 
polynôme caractéristique, y racines d'une équation transcendantale en nombre infini et 
dénombrable (on peut donc les indexer par l'entier n=1,2,...) et que A<=(2m+1)/2. Les fonctions 
angulaires ont été construites par la récurrence des coefficients du développement définies plus 
haut, le polynôme caractéristique s'en déduit également. Les parties réelle et imaginaire 
conduisent à deux fonctions de parité différente. 


Je n'aborde pas la question de la normalisation de ces fonctions angulaires, sachant qu'il s'agit 
d'appliquer sur une des solutions la norme préconisée dans ce document. 


Transformation de l'équation différentielle hyper-sphéroïdale : forme de Liouville et lien avec 
l'équation différentielle de Mathieu 


Par le changement de variable suivant, l'équation différentielle prend la forme dîtes de Liouville : 


er El ell lech x?) CE 


1— x? 


x = Cos(3) y(x) = (1 -x y BEN 8 (sin? (9) ` g(a) 
z ADO +y25in2(9)- Hi om) (QE o 


Sin’ (9) 
dene, e uusii) 


Lorsque u =1-À ou u=- À, l'équation différentielle devient une équation de Mathieu : 


2 2 
eigtl n LG Cos(289)]g(9)- 0 
2 


2 
gl  a=ø+7% +7 => g"(9)+(a-2qCos(29))g(9)=0 


Lorsque u =1-À , la normalisation des fonctions hyper-sphéroïdales se simplifie : 


Ÿ kcal d 24+2v Ifv+2k+u+24) T(v-u+1) = 
keete REESEN ENEE 


ais Dai) =1 
k=-0 
De ce fait lorsque u =1-A, il s'agit également de fonctions hyper-sphéroïdales dont les solutions 


peuvent être construites à partir des fonctions de Mathieu, ce que nous allons voir dans le point 
suivant. 


Les solutions exactes lorsque u=1-À 


Partons tout d'abord de la relation de récurrence initiale : 


adit oe) (B, -o0 azto )+C aż )=0  ke{-,..,0,1...,+} 
A 2 (v u+2k 1Xv u+2k) 
= Y Qvr4k+24-2)2v+ 444244) 
(v+24Xv+24+22)+u?-1+(22-1{u +2) 
(2v+4k+22-2)2v+4k+22+2) 
 (V+u+2k+24v+u+2k+22+1) 


(2v+4k+22+2)2v+4k+22+4) 


B, =(v +2k\v +2k+24)-27° 


C; 


Lorsque u =1-À cette relation est particulièrement simplifiée : 


2 
A=G = B, =(V +2k)\y +2k +24)- 7- 


2 


= a H)-(@+ 2x +22- 0e) ) aig) aii) 


De la transformation de l'équation hyper-sphéroïdale, on peut penser que de nouvelles solutions 
exactes sont obtenues lorsque u =1-À à l'aide des fonctions de Mathieu. Comme dans le cas où À est 
demi-entier, on peut utiliser la définition des fonctions associées de Gegenbauer : 


1-24 1-24 
1922 de E sn 1922 x 124 AL 
ci” (e)= C17 SaN e Coh e)=C17 mao A i 2. (z) 
Il vient : 
À eg 22 Ve ( dE A 1,1-4 Ee dE j 
=1-4 > Co) TZ) 4 P ed ET EI e ae 
pet GG) EG) € Al GG) CN Led e A 


+ 


Or l'on sait que : pz CE 2 Cos((2 EE Q? je = Sin((A EE 
Zus mT Dh Zus 2 Dk 


D'où: y=1-2 — 


hps” pe)= E di he Jehu) 
Ce qui fait que les deux développements : keng s'écrivent ici 
hqs} ” (x, 7° ) = £ (- 1) de LG ) Cr (x) 
k=-—0 


HI Dee (x, y? ) =( Ir: D (1 X (E D (-1) d 7 (7° }Cos((2 +v)ArcCos(x)) 
sous la forme : ec 


bas" (x, y’ ) = ( ee mT Ga (1 x? 1: De 1) ae (7° }sin((a + v)ArcCos(x)) 


Cas particulier lorsque v+ est entier 


Dans ce cas la relation de récurrence est invariante par le changement d'indice suivant : 


2k->-2v-2À -2k : 


L arnih) (v+ oa ant ae) aih) A 2)=0 


(v + À) k => 7 a 4 ee? LE USER +24+ 2k Ẹ K Sai 2 = (ai (2)+ as (7°) 0 


Dès lors la situation est en fait plus complexe qu'il n'y paraît. En effet lorsque l'ordre v+À est un 
entier la construction des fonctions de Mathieu est en générale double. Il existe deux valeurs 
propres caractéristiques de cette équation notées a et b. Comme l'équation des ondes hyper- 
sphéroïdales est liée à celle de Mathieu lorsque u=1- (voir forme de Liouville), il en est de même. 
Et les solutions hyper-sphéroïdales se construisent à l'identique avec les fonctions de Mathieu 
angulaire « ce » ou « se » d'ordre entier pair ou impair. 


De plus l'invariance de la relation de récurrence par le changement d'indice implique tant la 
symétrie que l'anti-symétrie des coefficients du développement. En effet en posant la 
proportionnalité : a (7°)= o dii" (°)= æ? =1 le coefficient a est égal à 1 ou -1 : 


OM 


Ja") 


A, bk a) où dit, bh 


—d° 
On note l'indice kọ: = 74> dirait) où 42 


Dans ce qui suit on remarque également que : 


Cos((v ++ 2k)ArcCos(x)) = Cos((v KEIER EREM 2k)ArcCos(x)) 
Sin((v + + 2k)ArcCos(x)) _ —Sin((v +A1—-2v-22— 2k)ArcCos(x)) 


Comme dans le cas des fonctions d'ondes hyper-sphéroïdales, il faut tenir compte du changement 
de signe du facteur (-1)* dans la transformation : k->k-k lorsque v+ est un entier : 


v+A pair GOU = 1 


k,=v+A—(-1) —(-1)°(-1) | 


v+A impair (— 1)? EE 


Lorsque v+À est un entier, il y a toujours deux valeurs caractéristiques w qui permettent de 
construire deux fonctions hyper-sphéroïdales dont les parités sont différentes. Ces deux valeurs 
caractéristiques sont choisies parmi les deux premières valeurs propres ordonnées de façon 
croissante et nous pouvons conventionnellement les associer au degré v et l'ordre u=1-À . Nous les 


ENNEN 


Ainsi les résultats suivants sont constatés : 


Lorsque v+À est pair alors pour la première valeur propre të (y? ), les coefficients sont anti- 
symétriques : di" (°)=0 er de (7°)=-d{/(72), et seul le développement de la fonction 
2 


hyper-sphéroïdales de première espèce est identiquement nul. Il ne reste que celui de la fonction de 
deuxième espèce : 


hps? (x,y 2)= ( y= AL x)? XCD EE 


bas dÉ y 2)= —(—1 Yarl- x "7 EN di (e )sin2lp + H)reCos(x)) 0 


WE 


v+=2p— 


Lorsque v+À est pair alors pour la deuxième valeur propre ohy) les coefficients sont 
symétriques : d1? (y?)=0 et di’, (7?)=aż1(7?), et seul le développement de la fonction 
2 


hyper-sphéroïdales de deuxième espèce est identiquement nul : 


hps” (x, y? )=( ee 27 (1 x)? 26 1) 42 (7° JCos(2(p + k)ArcCos(x)) + 0 


v+=2p E i 
a CN sëch X dE LE di G) snp + HareCos(x))= 0 


Lorsque v+À est impair alors pour la première valeur propre féie CAE les coefficients sont 
symétriques: Ko aw et d2()=dait(,2) mais avec le facteur (-1)® seul le 
2 Va: v, 


développement de la fonction hyper-sphéroïdales de premier espèce est identiquement nul : 


hps"? (x, y? ) = ( TE CAL (1 x’ y$ SS (- 1) SE (7° ) Cos((2(p + k)+ 1)4rcCos(x)) = 0 


v+=2p+l— p 
hqs? (x, °)= —— (et: Se SET D >: C fe di 0) sinp + k)+ 1)ArcCos(x)) #0 


Lorsque v+À est impair alors pour la deuxième valeur de œ ab) les coefficients sont anti- 


v ,ko—k 


GEN - ko 
symétriques : Den et dil2(y)=-d{i"(,°), mais avec le facteur GT seul le 


développement de la fonction hyper-sphéroïdales de deuxième espèce est identiquement nul : 


hps?" (x,y 2)= ( De CAL (1 x? 1: 21 1) di (7° ) Cos((2(p + k)+ 1)4rcCos(x)) #0 


v+1=2p+1— Si GE 
gt) "ur Ae Gel e? EC ai last Elie DreCos()) = 0 
LE 


Compte tenu des parités évoquées ci-dessus, les développements se simplifient. Lorsque v+À est 
pair pour la première valeur propre œ #t-4(y2) il vient : 


DEZ? 2 BEZ? 2 ae 2'7 2 ZE k 74,1- 2 H 
v+2=2p— hps; (x,y )=0 hqs; (x,y )= ( 1 "eer x ) 2 KEE d: (y }Sin(2(p + k)ArcCos(x) 


k=ko 


Lorsque v+À est pair pour la deuxième valeur propre œ Sta vient : 


A (- 1)°° Ce D ) Cos(2(p +k JArcCos(x)) + 


1-4 
2=2 hps ™ (x,y )= yi2 (l-2)? EE À 
ET E 
Datt 


Lorsque v+À est impair pour la première valeur propre w Ce Fil il vient : 


1 LA À k=+% 


v+i=2p+1— hps"” (x, y? ) =0 bas"? (x, y? ) = ( IW m TO) D x° J2 > (- 1} di’ Li ) Sin((2(p + k)+ 1)4rcCos(x)) 


Lorsque v+À est impair pour la deuxième valeur de œ Sg Fil il vient : 


` EE a (- 1}° dý GL ) Cos((2(p + k,)+ 1)4rcCos(x))+ 
v+4=2p+1- hps? (x, Y )= (-1) 2 AL -x ) 2 +2 Soy dili (2) Cos((2(p +k)+1)4rcCos(x)) 


k=ko+l 


hqs?"* (x, 7°) 0 


Voyons maintenant de quelle manière ces deux développements sont effectivement liés aux 
développements des fonctions de Mathieu de première et deuxième espèce. 


Liens directs avec les fonctions de Mathieu 


Les fonctions de Mathieu d'ordre v quelconque sont solutions de l'équation différentielle de 
Mathieu et possèdent un développement en fonctions sinusoïdales suivant : 


g''(x)+(a—2qCos(2x))g(x) = 0 


S ce, (x,g)= Yci,(g) Cos((v +2k)x) 

me, (x,g)= De (a)ett#x > EE 

SR se,(x,g)= X eielgl Sin((v +24)x) 
k=—00 


Avec un tel développement alors les coefficients du développement suivent la relation de 
récurrence suivante : 


q ed, eil 28 ale, Ga S0 )= 0 


Sachant que la normalisation usuelle pour les fonctions de Mathieu de première et deuxième 
espèce est la suivante : Sc, Golf =1: 


Keen 


Or lorsque u=1-A, la relation de récurrence des deux fonctions hyper-sphéroïdales de première et 
deuxième espèce redonnent bien la même relation de récurrence que celle des fonctions de 
Mathieu par les substitutions suivantes : 


r aaah BI E EE +24)- r - ol) dé (y }+ die) = 0 


Posons v=r—-À a = oly?) LA? + E d cO d d GL 2k) läit? (7° )+ q ae D 


Posons (-1} d4 (y)= c5,(g)— q ca ala)+ (+ 2k} — ak, (g)+ a chala) =0 


2 
Avec plus précisément a— a, (a) © — ©” GL a, (a) = À +o (7°)+ a 


d 2 
De plus la normalisation des coefficients du développement est identique. 


Cas v+À non entier 


Les fonctions hyper-sphéroïdales sont strictement proportionnelles aux fonctions de Mathieu, et 
comme la normalisation est identique l'on peut écrire directement lorsque v+À n'est pas un entier : 


pps a eye) (AE 2 (1-2) ce, sel 


r (4) 
SA À 
DEA AE —(-1 rl) sena) 
2 2 2 
a=% ai aile: Ai -5 = el Ai 


Remarque les valeurs caractéristiques a et b de l'équation de Mathieu sont identiques lorsque 
l'ordre v+A est non entier et les deux fonctions de Mathieu d'ordre v+À non entier forment bien 
deux solutions indépendantes de l'équation de Mathieu. 


0 


Cas v+À entier 


Par contre lorsque l'ordre v+À est entier la correspondance est plus complexe, car l'étude des 
valeurs propres de l'équation hyper-sphéroïdales prouve que l'on atteint à chaque fois les deux 
valeurs caractéristiques a et b de l'équation de Mathieu : 


2 Bi (a) < ds (a) 


a) a. E 


La méthode de construction des fonctions hyper-sphéroïdales est donc plus complexe car il faut 
tenir compte de la parité des solutions et du fait que pour un même indice v+À il y a deux solutions 
avec deux valeurs caractéristiques distinctes. 


Lorsque v+À est pair alors pour la première valeur propre œ a1-A(2) : 


= À —À À 
v+4=2p>0 > hps?” (x, vk 0 ` bosci" (x, MË ( (Ze SC D A P: se, (x,q4)=-( Jr o D x? E se, (x,q) 


Lorsque v+À est pair alors pour la deuxième valeur propre œ sm Call 5 


1-1 À 1-1 À 
v+4=2p > hps ™’ (x, Cl ( TE SC seCh SC (x, q) = ( he? SC x ye ce, (x, al bas" (x, y’ )= 0 


Lorsque v+À est impair alors pour la première valeur propre œ ee Tal ; 


SA ER —À À 
v+4=2p+1—> hps"” (x, vk 0 has?" ? (x, MË ( e Sr D x? P se à (x,g)= 1Ÿ 3x SC D x? Je ess (x,q) 


2 


Lorsque v+À est impair alors pour la deuxième valeur de œ ie Tal g 


1- 1- 


2 2 à 2 
v+4=2p+1— hps; (x.y?)=( (e SE x]2 de Aa (e Se D E ce,,1(x,q) hqs; ™ (x, y’)= 0 


Pour déterminer les fonctions hyper-sphéroïdales correspondantes dans le cas UA il suffit 
d'appliquer la transformation u->-u+1-2A. Or pour l'instant la formule de transformation n'est 
connue que dans le cas où À est demi-entier sous la forme suivante : 


HER hp) A an + n)a) mosi” (E)Z Cos(( +) pa) 
GR T(v+u+24) T 
— 
ES e PE Deeg E FED (sQ + u)r) hpg}"(2)-T Cos((à + u)r) hps?" H 
ù T(v+u+24) 2 g 
E Ne (EEN Sage 
veR u=1-1 hps; (z)= T (3 hpa? (z) 
Avec u=1-A, il vient : 2p+1 4 (+ g ) 
e EE hpqż™>(z)= Z NA, Säit) 


2 T(v+1+1) E 


Pour déterminer maintenant les fonctions de type 2 (pour z>1), appliquons les formules de 
connexion avec les fonctions de type 1, sachant que Sin(n(A+u))=0 : 


H 
VER u,AeN> hps" (z)=(-1} k-i} hps" (z)= LE ln? (z) 
DP (2-1) 


£ u 
vueR À GC Dech 3 hps! (e)= Des (2) i GE hps*(z) 
2 


Ds 


16 E de 2 É S 
TS pare E ige ne Een Ziel 


Construction des fonctions hyper-sphéroïdales dans le cas énéral, non demi-entier et v demi- 


entier 


On a vu qu'avec un développement sous forme de fonctions associées de Gegenbauer, la 
récurrence comportait des pôles pour les valeurs de v+À entières. La récurrence habituelle était de 
cette forme : 


Adtt(7°)+(8, -oœ(y°))d#(72)+cd!#(7°)=0 kef-,.,0,1,.,+00)} 
(v—u+2k—1Xv-u+2k) 
T Gy+4r+22— 2)2v+4k+22-4) 
(v +2kv +2k+24)+ u? —1+(24 -1u +24) 
(2v +4k +24 -22v +4k+2A=+2) 
o (v+ u+2k+2A)\v+u+2k+24+1) 
(2v+4Kk+22+2)2v+4K+22+4) 
1 1 1 
2v +4k +24 -4° 2v+4k+24-2/2v+4k+22+2/2v+4k+21+4 


A; y 


B, =(v +2k\v +2k +24)-27° 


C, 


Pôles possibles sur 


Il y a un procédé très simple permettant d'éliminer ce type de pôles, il suffit de proposer un 
développement utilisant directement des fonctions de Gegenbauer. 


hps” (x,y°)= -xY 


1727 1 +o 


KE 2)c Cå a(x) avec A=1-2-u 


lorsque l'on injecte le développement proposé dans l'équation des ondes hypersphéroïdales, à 
savoir : 


(=) a") 22) ra) 0+7 -)- ge SE y(x)= (0) 


1 


y(x) = (1 se Ke 20) A=l-i-u— a 2A)x g'(x)+ (2+ CD g(x)=0 


gl) = 3 ei, (2) 


Alors cela justifie de développer la fonction g(x) sous la forme : E . Et par 
l'application successive des relations de récurrence sur les fonctions de Gegenbauer : 


est. EDEN (+ GA +v HE 

” 2(A+v) 

(v+1)v +2) ^ (2) v?+2Av+A-1 dë (2A +v -12A +v -2) na G 
4(A+v)(A+v +1) © 4(A+v-1)(A+v=+1) ” 4(A+v)(A+v-1) "' 


= x? C^(z)= 


On obtient la récurrence suivante : 


Aert, Oli, rg =0  kEf-00,..,-1,..,0,1,..,+00} 
A=1l-2-u 
Q = © + Ai N 
ei (v+2k-—1)(v +24) 
4 (v+A+2k-2Y{v+A+2k4-—1) 


+2kYv+24+2A)+(A-1X1+2A) 
(v+A+2k-1{v+A+24+1) 


_y? (+2A+2k+1){v+2A+2k) 
4 (vV+A+2k+1{v+A+2k+2) 


A, 


2 
B, =(v+2k\v +2k +24) S br 


C; 


Remarques : les conditions pour lesquelles la récurrence ne présente plus de pôle sont liées au fait 
que la somme v+A ne soit pas un entier. L'intérêt d'avoir introduit cette récurrence réside dans le 
fait que la somme v+À soit un entier. Dans le cas contraire le développement avec les fonctions 


associées de Gegenbauer est suffisant. Or À=1-1-u—v+A= (~ +4)+1-4-(4+ u), Lorsque 
u+ est un entier_ alors il suffirait que 1-À ne soit pas un entier, mais même si la récurrence est 
définie, nous verrons rapidement que les fonctions de Gegenbauer de première espèce ne sont pas 
définies dans ce cas. Donc le développement introduit ne résout pas le problème de la construction 
lorsque u+À est entier et v+À entier. En revanche pour le cas u+À non entier et v+À entier, on 
constate que le développement portent sur tous les indices k de -2° à +20 car les coefficients A, et Ck 
ne s'annulent jamais. Cela permet de compléter une partie de la construction des fonctions hyper- 
sphéroïdales. Le deuxième cas où la récurrence comporte encore des pôles est lorsque 1-2À-u est 


un entier, soit u+2À -1 est un entier. 


Exception faîte des remarques précédentes, il n'y a aucun souci particulier pour la construction des 
secondes solutions à l'aide des développements dont les coefficients suivent la même relation de 
récurrence : 


u+22-1 +œ 


basis" 2 Sables.) avec A=1l-2-u 
Ärem 


Les fonctions de Gegenbauer de première espèce (polynôme notamment lorsque de degré est 
entier sont définis ainsi : 

1-2A 

CC I-2A 1l 
SE Zelbrääain 2 eh 


SOS ERNEST E 


C'est bien un polynôme car le terme non polynomiale de la fonction associée de Legendre 
1-2A 
ek 
compense le facteur : 


hypergéométrique : 


. Cette fonction est également donnée à l'aide de la fonction 


c*(e)= TV +2A) F(-vvs2aasie) Kë 
T(v+1)7(2A) SE 2 


Par extension, les fonctions de Gegenbauer de deuxième espèce sont définies à l'aide des mêmes 
fonctions de Legendre de deuxième espèce : 


1-2A 


27 JrI(v+2A)f 222 La : 
rv +1) (A) ( z) Piw ze] II 


El 


Mais attention dans le cas où A est un entier négatif ou nul alors la fonction de Gegenbauer est 
identiquement nulle et ne peut être utilisée dans le développement. Cela correspond au cas où u+À 


est un entier strictement positif 


Pour z > 1, les définitions sont : 


1-2A 


ap 22 VrI(v+2A),, JE a 

Ch T(v +1)r (A) E-1)+ 7 Sa Dan 
A 22 Jr I(v+2A , YA La 

Co (z)= En ){ zi) 4 o? Jl SSC 


Et avec les fonctions hypergéométriques : 


Ota T(v+u+1) LG dä vrsi 3d 


3 S e E - v+uz-1l,-2,-3 
zarva) E 


Ou encore : 
1 
A 2 (Ah de T(v+2A) v+2A+1 v+2A SR 
Cop E Za ltbe +A +1) z^ 25 | 2 3 2 ;v+A+l]; 7? 1<z<o 


ou bien encore : 


(z^) Var (v +24) Ve fie 1 
cA o) dl ) e mli- -) A A av) 1<z<o 


Même si l'on connaît la normalisation des fonctions de Gegenbauer et la normalisation des 
fonctions d'onde hyper-sphéroïdales, on ne peut déduire une relation de normalisation sur les 
coefficients du développement. 


Les expressions des dérivées premières et secondes de ce type de développement sont les 
suivantes : 


DEET 


(1-x°) S hps?" (x,y°)= S gt 0) ctal) 
k=—00 
CH fa- E ps D = 52A Zer (cial) 
k=—%0 


d Aaf? HD Saala gio) 
k=—00 


ou bien encore : 


+24-1 Ae 
ri 2 Sgh) Cia) 
(ess 


u+2A-1., +o 


£ Ire" (x, y }= (1 =x? (SC Zei (2) Al ER hs (x)+ x(u +24- 1) er (x)} 


hps?” (x, y?)= (1 x?) 


2 u+22-1 , + 


FS TA E }= l- Zei (2) hı -x 4A(A HINCA elt delt 2) (u +22 0) A CA (x)+(u +241 


Il reste maintenant à construire les solutions pour les deux cas : 


Cas n°1 : u+À est entier et v+À entier 
Cas n°1 : u+2À -1 est un entier et v+ entier 


On a vu que pour u+À=1, le cas v+À entier se construit avec les fonctions de Mathieu périodiques 
(ou anti-périodiques) de période n et 2 n. C'est ce que nous allons essayer dans le prochain chapitre 
de ce document où je vais étendre cette construction à l'aide des fonctions de Mathieu périodiques 
associées. 


Cas limite y->0 pour un degré v quelconque et un ordre u+À entier : petite étude des fonctions de 
Gegenbauer associées sur z€[-1,1] de paramètres _(v,À,1- À) 


On sait que dans le cas où y=0 l'équation des ondes hypersphéroïdales se ramène à l'équation 
différentielle des fonctions de Gegenbauer associées : 


(-22)»"()-(1+2)z »@ (4 +24)- Ss 2) y(e)=0 


1+22-1 
y(z)= (1 — 7° F 2 g(z) A=l-4-u y(z)æ C?” (z) si première espèce 
(1 -z )g"(z)- (1+2A)zg'(z)+(v+2-AXv+2+A)g(2)=0  gļlz)æ CA, a (z)= CD, zl si première espèce 


v+u+22-1 


=> < ou bien 
(1-22 )g"(z)-(8-22-2u)z g'(z)+(v — u+1\v +u+22-1)g(z)= 0 


En poussant plus loin avec une forme de Liouville et spécifiquement pour u=1-À : 


z= Cos(t) 
| 2-1 >u" (t)-2(u +2-1) Cotan(t}u'(e)+(v + u +22 E DÉI 1)u(s)= 0 


g(2)= (S) ` (Aule) 


u"'(t +(v Alz =0 
u=1-2> u(t)= Cos((v +2)r) ou Sin((v +2}) 


a(2)= (sn (0) (Duo) 


Cela suggère donc que la fonction de Gegenbauer associée en question a la forme suivante : 
À 
lee D — 7 IS Cos((v + 2) ArcCos(z)) 


Or d'après une formule valable lorsque À est un demi-entier et pour le paramètre u quelconque : 


1-22 


= 1-24 , 1, 
huer A-Z pen CH) EE) EO 


ESEE a= 2e 


Comme „5 O SE 1 2 26 +1)ArcCos(z)), il vient : 
SÉ E E SS 


ds Wë peN C#(2)-(- TE 274 Cos((v +2)A4rcCos(z)) 


BORNE 


22-1 
Par mesure de simplification on peut supprimer le terme ks 1) ` car par la suite nous n'allons que 
regarder la forme des fonctions de Gegenbauer associés des ordres voisins : -À, 2-À,3-À, etc. 
ce qui je l'espère n'aura aucune influence sur les résultats obtenus. Partons donc de l'expression : 


aia (2) z Eet Cos((v + 2)ArcCos(z)) 

WÉI D ay E 

(2) Sin((v + à)AreCos(z) 
Zb 


Pour obtenir l'expression pour la fonction de deuxième espèce, je tire profit de l'expression du 
Wronskien des fonctions de Gegenbauer associées: 


Vå 


CO (z) =Á 


H ENT 27 T(v+24+u) a 2 
vk (z Lee )}= Tr(a} TOFI) D ) 
CC CS EIER 


Il vient alors : 


whch (z), Coi GH = T) í : 7} W(Cos((v + 1)A4rcCos(z)) Sin((v + 2)4rcCos(z)))= - 


m (v +4) BE 212 40) a 
AT l-z°j ? =-4— — (pv +14)> AU 
LOT rG) 
L'expression de la fonction de deuxième espèce est donc : 
gek ae o o) 
2 IEN D w e k 


114 (2) nor Cos((v + 2)ArcCos(z)) 


Résumons les expressions obtenues : ui e H Se 
CA (z)= 8 T Sin((v +2)A4rcCos(z)) 
SH 2/T() ab 


; ; e x , , Ciz) cez). A 
Avant d'établir un récurrence à trois termes sur les fonctions Y ) (Gi ) , il faut mettre à 
profit une formule de liaison par la transformation de l'ordre u->-u+1-2À, pour débuter la 
récurrence par la connaissance d'au moins deux de ces fonctions : 


CE (z2) z T( TS } ET + u)r)C}* DÉI S Cos((4 + uja), H 


Thinz 34 


AA a vlek (Ga rc") 


BEEM 2 T(v-u+1) 2 
EE Coell/ Co 
v (z) e TVv+u+24) os(( + u)r) g) 
a T(v-u+1) 
2 T(v+u+21) 


CA (z) A 


Si A+ueN> 


Cu Cos((4+ u)r) C> (z) 


Formules que l'on peut mettre en relation lorsque à=1/2 et avec celles-ci : 


P-"(z)= 2 T(v-u+1) “() 
seb) z Sin(ux)T(v+u+1) _ 2n+1 JEN 
| QG) $ Stol To az Hetz? Ge 2 
Á 2  I(v+u+l)” 


On peut utiliser ces dernières relation pour calculer les fonctions sur le jeu de paramètre suivant : 


eent 2 TEA) ee à nr) ct G) ut, ACEN 
rx T(v+u+2À) E v A YLA 
À,-u+1-24 o T T(v — L +1) ji R D Chtätz 
Coy e a PARC SÉ EMA 
Crit 217 Cos((v+2)4rcCos(z)) CA (7)= 217  Sin((v+2)4rcCos(z)) 
| v IEN E v (v +1) (4) Ab 
VAER 
S 114 x  Sin((v+2)4rcCos(z)) 4 p T _, Cos((v +4 )ArcCos(z)) 
Ce ` Gleck 2 S 
2 r(1) l-2} (v +1) (4) IG 


Tirons maintenant une loi de récurrence à partir de la récurrence sur les ordres pour les fonctions 
de Gegenbauer associées. Utilisons pour cela une correspondance entre les fonctions de 
Gegenbauer associées et les fonctions de Gegenbauer (même si cette correspondance n'est valable 
que pour les ordre u entier) : 


ÉT(A+u) 


VUeN C#(z)=(-1)"2# (1-22) C4 (z) 


SE" 


Les fonctions de Gegenbauer suivent la récurrence suivante : 


Ai -z X2 + ua + u+1)C 4 2(2)-2 (A+ u(i +24 + uee e) +v = uv -u+ 24+ u)! (z)=0 


v-u-2 v-u-l v-u 
Alors les fonctions de Gegenbauer associées suivent la récurrence suivante sur les ordres u voisins : 


1+2(4+u)) 


Cae (z)+ 2 ( SE (z)+(v-u\v+u+22)C*(2)=0 
l-z 


A,n—À 
d S zZ ` z $ 
Selon cette formule alors les fonctions ~” ( ) suivent la récurrence suivante : 


EC EH Co (z)+(v+2-nÂv+2+n)c"()=0 Cl -Z(+ 2n) us 2)-(v+2-nfv+1+n)C?(z 
leede E E 


S EH ce (z)+(v +2-nv +1 +n) Cy (2)= 0 E (z)= - car G-A DE D +4+n) C (z) 


Ou 0 l g 
DÉEN 1-7? 


E (z ) + 


Lorsque À=1/2 cette formule est à mettre en relation avec la formule de récurrence pour les 
fonctions de Legendre associées d'ordre demi-entier : 


1 2 p(t) p 1 RE oS 
u=n a Dr P, 2(z)+ Fast GE n |P, 2(z)=0 


S x e Z 
Cette relation va permettre de calculer les expressions des fonctions ~” sur les 


ordres successifs. Par exemple pour u=3/2 : 


ETAT SCH Sin((v + 2)ArcCos(z)) AA 217 Cos((v +4)ArcCos(z)) 
Gel hu ATI 7 e 8 Sei TU) GEN e 
21 m -a Cos((v +2)A4rcCos(z OT EN m . Sin((v +1)ArcCos(z 

Cor (z)= Ge ara 6 SE ý Cio) (z) XTA) b Wë D 


+ AW- 2 Sin((v + 2)4reCos(2)) 
+ AW- 2 Cos((v + 1)4reCos(z)) 


+ À )Sin(ArcCos(z))Sin((v Bä A)ArcCos(z)) 


+ )Sin(AreCos(z))Cos((v + 4)ArcCos(z)) 


v 2 Cos((v + 1+ 2)4reCos(e))] 


v-l+4 Sin((v +1+ aJareCos(=)) 


) E EEEE E E A) 


Trsa 


m=l+1 


m=n-l 


[Ir +2+m 


rbl bes CA) AF) 
l-z 
Ciel Ova aC) 
DAA 
CE u o L z Cos((v +4)4rcCos(z))-(v 
= riz) 
Chr (2)= Te sv +A)arecos(z))-(v 
2rGK-2)2 
z= Cos(ArcCos(z)) EN Sin(ArcCos(z)) 
212 
en (z) = —— GE (- Cos(ArcCos(z)) Cos((v +A)ArcCos(z))-(v 
riz)? 
=> 
Ci (z)= = SSES (Cos(ArcCos(z)) Sin((v + A)ArcCos(z))- (v 
TG 2)? 
1-2 
Cr A(;)= S SS | Re Cos((v -1+2)A4rcCos(z))- 
z riz) S 
Con, (z)= 4 T (== Sin((v 1+2)4reCos(z))-" 
2 r(a)1- 22)? 
On peut donner une forme très précise à ces fonctions d'ordre demi-entiers quelconque, sous la 
forme : 
WEE EN I=n-1 m=l m=n— 
c?" (z)= A+n-1 D ( de E j< E (Te ; À WE „| 
RÉI (BA 1=0 (n- m=1 = 
e Kee e l=n-1 vg: 
Cr ()= e 1-7)! s nm] 


GD — z?) 2 


| Silbe + 1 + 21-(n—1))4reCos(z) 


m=l+1 


Remarque: comme n=À+u, alors dans l'hypothèse où v+À est un nombre entier, on a v+À-n+1= v- 
u+1. Dans ces conditions le développement de la fonction de première espèce est en cosinus impair 


si v-u est pair 


Voyons si l'on ne peut pas trouver un développement de la forme 


et en cosinus pair si v-u est impair. 


l'injectant dans l'équation différentielle : 


u' '(t)- Au +21) Cotan(t}u'(t)+(v + u +22 -1v 


u=n-&u' '(t)-2(n-1) Cotan(t)u'(t) H (v+4 +n 


u(t i (v+4+21-n+1)t en 


Lag 


— u+1)u(r)=0 


1)(v+2-(n-1))u(r)=0 


Ce développement conduit à une récurrence à deux termes : 


1=0—>(n—1Xn-1-v-2)4" D zwet Al A7 
1=1—(n-2X{n-2-v-1)4; =2(2+v+2) A4! 


(n—1-1Xn-1-1-v-2)4%, =(+1{7+1+v+24) 42, 


Cette relation de récurrence s'arrête lorsque l=n-1 et tous les coefficients d'indice supérieurs 
s'annulent. Le développement s'arrête également pour 1=0 du fait l'annulation du membre droit de 


1=+0 


la récurrence pour l=-1. Cela justifie donc d'écrire le développement fini „(t)= D A 


i(v+4+21-n+1) , Le 
d 


1=0 
As _(n—-1-1{n-1-1-v-2) 
A (2+1X72+1+v+2) 


rapport entre deux coefficients d'indice qui se suivent donne : 


Le même rapport dans la formule donnée se calcule ainsi et donne le même résultat : 


REÇUS CG, (fiv DEE "(Te +A+ m) 


m=1 m=l+1 


B, =C ec AT A-n: ell 


m=1 m=1+2 


Ba n-1-lv+2-n+1+1 
B, (EA v+2+1+1 


La récurrence établie est donc celle qui construit les fonctions de Gegenbauer associées d'ordre 
u=n-À, de première et deuxième espèce en prenant les parties réelle et imaginaire, à une constante 
près bien évidemment. 


Valeurs particulières des fonctions associées de Gegenbauer telles que u+À est un entier pour 


certaines valeurs de l'ordre v 


Pour certaines valeurs de l'ordre v, les expressions établies précédemment se simplifient : 


Ca ()= Re Kat E EE Ti A-n4 "(Te +A+ m) Cos((v +4 +21-(n -1))ArcCos(z)) 


r(ai-z) 2 1=0 I(n-1-/)! m=l m=l+1 
2-22!" m l=n-1 (n sch 1)! m=l m=n-1 

CE SDA : (fi A-n: m| [[v+4+ m) Sin((v +4 +21- (n —1))ArcCos(z)) 
{ais}: HR AL 


Pour v+À=1-n, le développement ne contient qu'un seul terme quelque soit la valeur de n>=0 : 


= CE (en lettre DCH ar 


X cir- (z) i 

p= =0 A A+n- 

HP riz) + r ra) D! 
Ce ,(2)=0 


Pour v+À=2-n, le développement ne contient également qu'un seul terme quelque soit la valeur de 
n>=0 (deux termes puis un seul) : 


m=n-l 


Į[@-n+m)=0 Ee ZER AA 
d à arni gi E nt 2A vw (27-2)! 
m=l+1 Qarz = f= 2 2 2(-1 1 2 Za pe 1- 2 2 1 1 
à eneh- 2200 Tir m) (ont 
(ET ER ECH (z) =0 
l-2 l-n 
chae): 
=, 25 
Pour n=1l ONE 
Cr (ee Tee 


r(a-2} 


Pour v+À=3-n, le développement ne contient que deux termes quelque soit la valeur de p>0 (trois 
termes se réduisant à deux) : 


m=n-l1 


Il E SET m) =0 Ae) 27 (- me? 

m=l+1 GAA = 1 2 SC 2 C 24 C 8 ve 

l£n-1 = n (z) ( Z ) EN Ga ll n) os( rc os(x)) (n )) zai 
a CZ (z)=0 

lzn-3 


2 1-2 


cè (z)= (L -=z E a Cos(2ArcCos(z)) 


Pour nz 


Iag 
Cora (z) T BI -7 E T. Sin(2 ArcCos(z)) 


T(1) 


Pour v+À=4-n, le développement ne contient que deux termes quelque soit la valeur de p>0 (4 
termes se réduisant à deux) : 


m=n-| 
El (4 -n+ m) =0 
m=l+1 
l#n-1 = 2 y (27-5)! 
ee Jecke? Sat? E li-n) Cosl3AreCosl))-6ln-1)Coslarecoshi)) 
l#n-3 (OA (z) =0 
(ET EE) 
i A 1-2 
C(2)=(1-2)7 SO Cos(34rcCos(z)) 
Pour n=1 
214 2 na" e 
C2 (2)=-(1-z ) S ra) Sin(3ArcCos(z)) 
A+ AA 
Crat-Ä-a2fz races Arcok) Cos(4rcCos(x) 
Pour n=2 
A1 T A1 


Cell": o (Cos(84reCos(x))-3 Cos(4reCos(x) 


On constate donc que quelque soit la valeur de n>0 pour les fonctions de première espèce : 

- pour v+À=1-n et v+A=2-n, le développement ne contient qu'un seul terme non nul. Pour v+À=1-n 
c'est un développement en cosinus pair et pour v+À=2-n un développement en cosinus impair 

- pour v+1=3-n et v+1=4-n, le développement contient deux termes non nuls. Pour v+À=3-n c'est un 
développement en cosinus pair et pour v+À=4-n un développement en cosinus impair 

- pour v+À=5-n et v+À=6-n, le développement contient trois termes non nuls. Pour v+À=5-n c'est un 
développement en cosinus pair et pour v+À=6-n un développement en cosinus impair 

- ainsi de suite … 

- et pour n=0 et n=1, le développement n'a qu'un terme. 


Deux remarques sur ce petit chapitre 


La première : il se trouve qu'une série d'articles de 1949 et 1950 abordent d'une certaine manière 
ce type de développement avec les fonctions de Gegenbauer dans le cadre de l'équation des ondes 
sphéroïdales. Il s'agit alors d'étudier les fonctions dies « associées de Mathieu a Voici au moins 
trois références : 

-1949, R.Campbell, « Sur une expression remarquable des solutions de période 2 nr Pi de l'équation 
de Mathieu associée » 

- 1950, R.Campbell « COMPORTEMENT DES FONCTIONS DE MATHIEU ASSOCIÉES POUR LES 
GRANDES VALEURS DES PARAMÈTRES » 

-1950, R.Campbell « CONTRIBUTION A L'ÉTUDE DES SOLUTIONS DE L'ÉQUATION DE MATHIEU 
ASSOCIÉE » 


Il me semble qu'après un examen plus attentif de l'article « CONTRIBUTION A L'ÉTUDE DES 
SOLUTIONS DE L'ÉQUATION DE MATHIEU ASSOCIÉE », ce dernier comporte « peut-être » des 
erreurs, notamment sur la relation de récurrence des fonctions de Gegenbauer, ce qui peut 
conduire à une erreur sur la relation de récurrence des coefficients du développement. A titre de 
comparaison, les auteurs L.J.CHU, J.A.STRATTON proposent dans l'article de 1941, « Elliptic and 
Spheroidal Wave functions » également un développement en fonctions de Gegenbauer (cas de 
polynôme), mais ces derniers n'ont pas la même normalisation, il est moins trivial de vérifier la 
récurrence (Construction des « V-Functions a formule 83). Toutefois la récurrence de L.J.CHU, 
J.A.STRATTON est non seulement plus cohérente que celle obtenue par R.Campbell, mais par un 
changement de paramètres est identique à celle établie ici. 


La deuxième : telle que présentée par R.Campbell, il est clair que les fonctions associées de 
Mathieu sont les fonctions sphéroïdales avec un argument Cos(t) ou Sin(t). C'est ce que nous allons 
voir, car pour le cas des fonctions de Mathieu associées périodiques la construction la plus 
« cohérente » est celle réalisée par R.L.Ince conduisant à une récurrence à 4 termes. 


Dor extension à l'étude de l'équation des ondes hyper-sphéroïdales, il n'y a pas d'intérêt à 
rechercher des solutions non-périodiques des fonctions de Mathieu associées sous forme de série 


sinusoidales puisque ce sont celles que j'ai abordé par le développement en fonction de 
Gegenbauer associées 


Construction des fonctions sphéroïdales dans le cas u+À et v+À entiers, lien avec les fonctions de 


Mathieu associées 


Une des formes de Liouville de l'équation des ondes hypersphéroïdales est obtenue avec le 
changement de variable x=Cos(t) : 


q- araa) oyla) AR CH Main x = Cosh) 


UE 


SE > wielt 21233 u) Cotan(t) w'(t)+ lo- u(u+22)+ 7° Sin?(t))w(s)= 0 
y(x) (Sin? (O) Z w(t)— w" (t)+2(1—2 -— u) Cotan(s) w'(r)- lo H(u+2A-—1)(1— u)+ y? Sin?(t))w(e) = 0 


Une forme équivalente est obtenue avec y=Sin(t) : 


(1 ve ale ll eech Si CH x = Sin(t) 


1— x? 


Si KE (Cos? (0) aa w''(t)—-2(4 + u)Tan(t)w'(t)+ lo- u(u+22)+ y°Cos” (1) w(t)= 0 
y(x)= (Cost) = 2 WOS w" (t)+2(4 + u—1)Tan(t) w (A+ (o+ (u +22 -1X1 - u)+y°Cos? (t))w(t)= 0 


Selon un article de 1923 de E.L.Ince, « Associated Mathieu Functions a l'auteur transforme de 
l'équation sphéroïdale en deux substitutions que je résume ainsi : 


D x2)y"(x)-2 x va) x?) 


u? 1 
ES z Mad x = Cos(t) D 


y(x) = CHOLLO >q= d a = © + : + A => Po +2q Cos(2t) A Je) = 0 
gO = (sr fu el sl = (sin (0) e) = (sir OO) + 20 (9 = (sir OP lat) 
=> Hi "(t)+ 2A Cotan(t) u,'(t)+ (a — AT +2q Cos(2r)}u ()= (0) 


sic) - (sr) > mA h) =h) ee) E mE) (sin (0) Zait 


> um," (t)+2(1— A) Cotan(t)u,'(t)+ (a D — AP + 2q Cos(2t))u,(t)=0 


Si nous appliquons une démarche similaire à l'équation hyper-sphéroïdale : 


q- EA aa) Boayl) AD aen w= Cos) 


UE 
= y(x)= = (Sin? (©) al >u" (t)+ 2(4 + ul Cotan(t)u", (t)+ lo- u(u + 2A)+ 7° Sin?(t))u (t)= = 0 
y(x)= (Sin? (0) = g u, (£) >u", (t)+ 2(1— A u) Cotan(t)u' 5 (t)+ lo + (u +22 — 1X1 z u)+ y?°Sin?(£))u BO = 0 
Par correspondance des paramètres nous écrirons alors : 


2 2 2 2 
den A=u+À a-N=o-p{u+24)+È = +7 re ann 
WOP 2A Cotan(t)u,'(t)+ (a — AT +2q Cos(2r)}u, (£) = 0 
1-u=1-A+2 
u+22-1=A+2-1 


u," (t)+2(1- A) Cotan(t)u,'(t)+ (a — D — AY + 2q Cos(2r)}u,() = 0 


= (u+22-1)1-u)=2-(1-A) 


Comme l 


Ainsi que l'équation intermédiaire : 


Q- h" 6)-0+22)x stall rech x?) Olai x= Cos(t) A=u+À2 


>= (sir D) 2 ORTER 


2 
J a=0@+% + 5 => AUDE +2q Cos(2t) Tat 0 


Sin’ (t) 
Les deux transformations intermédiaires introduites par E.L.Ince s'écrivent ici sous la forme : 


3 D Aru A 
(Sin? (c)) 2 g(t) = (Sin? O} u(t) > g(t) = (Sin? (0) 2 u(t) = (Sin? ()}2 u(t) 
EI _u+2A-1 l-u—-2 EA, 
(Sin? (0) 2 g(t) = (Sin? (0) 2 H: (t) <> g(t) = (Sin? (e)) 2 = (Sin? (0) 2 
Nous retrouvons donc toutes les équations données par E.L.Ince mais cette fois-ci dans le contexte 
de l'équation hyper-sphéroidale. Si l'on pose À = 3-y+1-22-1- 4-2, le système d'équations 


est identique car l'équation est invariante par le changement À 1-A&u—-u+1-22, 


Résumons donc le résultat obtenu après ce dernier changement : 


2 2 
D 1 y 

A —l À = EELER 
à S A 4 2 


u+22-1 


CONCHO 


y) = (Sin (O) u) = u" (€) + 21-A) Cotan(t)u,'(t)+ (a -0 -AF +24 Cos(2t))u,(t)=0 


u (£) => Hi dat 2A Cotan(t) u'(t)+ (a — AT +2q Cos(2t))u, (£) =0 


Remarque : sur cette généralisation de l'équation de Mathieu on remarquera que le paramètre q 
est de signe inverse dans l'article de E.L.Ince. En effet communément l'équation de Mathieu se 


note : 
y'""()+(a—2q Cos(2r)) y(r)= 0 
Et par extension celle l'équation de Mathieu associée : 
y""(#)+2A Cotan(s) y'(t)+ (a — À — 20 Cos(2r)) y(r)= 0 


Partons maintenant de ces deux équations proposées par E.L.Ince en recherchant la solution sous 
la forme de fonctions périodiques ou anti-périodiques de période x par des développements 
identiques en la forme à ceux employés pour les fonctions de Mathieu. Prenons par exemple une 
solution paire et anti-périodique en t de période x : 


1=+00 
Di (£) = ` A7. Cos((21 + 1) t) 
1=0 


Injectons cette forme dans l'équation différentielle : 
Sin(t) fu, "(t)+2A Cotan(t)u,'(t)+ (a — AT +2q Cos(2r)) u ()}= 0 


AX,,Cos((21 +1) 


Donne sur le terme d ‘expression suivante : 


A (ag Sinell -1)t)+ (q+ (+27 AY -a) Sin t)- (q +(+ 21+ AF — a) Sin(2(t +1)r)+ q Sin(2(1+2)r) 


Qui doit s'annuler une fois les termes de l'expression réparti sur les termes voisins du 
développement (ils sont tous en fonction Sinus dans l'équation différentielle de départ). 


Cela conduit à la relation de récurrence à 4 termes suivante : 


1=0 > (a-(a +1} )4? +(-a+g+(A-3Ÿ) 4-9 4? 20 
1=1 >q A7 + (a (A +3} allan +(-a+(A-5Ÿ +4)4 +4 4 =0 
1=2 >q 4 Ale (asf -q)4ż +(-a+(A-7F +q) 42 +q 44 =0 


o Ai, nl DAT 4q) A | ( a+(A-21-3) + elAi A}, =0 


Sous forme matricielle et en changeant de signe pour l'ensemble du système d'équations linéaires 
on obtient : 


DAAT -g-G-A) q 0 S 0 -1 0 0 0 
-q  q+G+A)Ÿ -qg-(5-A) q n De 0 -1 0 … 
0 -q g+(5+AŸ -4-01-47 q ss -1 0 … 
0 -q -q-(21+1-Ař 0 eg, 0 EE e eet 
0 -q qg+(21+1+AŸ … ma e 2070 
| 0 | {0 0 0 …| 
[1 111.. 1] 
CE CAE 1l 
= 10 0 1 1 b,,=0 i>j 
S(M-aB)A=0S(B'M-al}A=0 B'- = [b] | 
0 0. 1 bone 
… 0 0 
0 0 0 


On voit donc qu'il s'agit d'un problème de détermination des valeurs et de vecteurs propres de la 
matrice B'M Les coefficients du développement sont les vecteurs propres calculées à partir des 
valeurs propres déterminées par diagonalisation. 


Prenons un développement de la forme : u (t)= ces (t)= KS ACos(21 t), Par le même procédé, on 
1=0 
parvient à la relation de récurrence à 4 termes suivantes : 
1=0—(a-A° -q) 4 +(-a+qg+(A-2Y) A? -q 4 =0 


1=1— A? +(a-(a+2Y -q) Garasa +4) 42 +9 A =0 
1=2 >q di +(a-(A +4} -g) 4 +(-a+(A—6Ÿ +4) A4! +q 4 =0 


q A24 (a (21+A) q)4ż +( a+(A -21-274 elAi -a A2, =0 


Aya 


Sous forme matricielle et encore une fois en changeant de signe pour l'ensemble du système 


d'équations linéaires on obtient : 


lq+4 -2q-(2-AŸ q 0 
-q  g+(2+A) -g-(4-Aÿ q 
0 -q g+(4+A) (6-A) 
0 -q 
e 0 -=q 
| 0 
S(M-aB)A=0@(B'M-a1)A=0 B'-|° 


0 


O oo me me 


0 1 -1 0 0 ONE) 
sé Së 0 1 -1 0 … À, 
0 Se 0 0 1 -1 O0 …|| 4, 
=d =0 
-(21-AŸ 0 „n0 0 ol 
qa+(21+A) o ép 00 0: 1 À 
| {0 0 o n] 
lolo 1] 
DE O 
De 
Ds 1 taib y= 
0 0 1 
0 


On voit donc qu'il s'agit également d'un problème de détermination des valeurs et de vecteurs 
propres de la matrice B*M. Les coefficients du développement sont les vecteurs propres calculées à 
partir des valeurs propres déterminées par diagonalisation. 


Dans l'article E.L.Ince introduit les notations des fonctions associées de Mathieu sous la forme 


suivante : 


WOP (a +2q Cos(2t) 


a) Se? 


g0) = ccen (ra) (sn (PS 45 
el ccena) (sin) Z5 Aia (<a) Cos((21 +09 
g(t)= sceta (tq)= (Si (7 


glt) = sce, a (t-q) = (SiP (O) 2 


= q) Cos(21 t) 


ee q) Cos((21 T 1) t) 


1—A ten 


2 A; E q) Cos((21) t) 


Nota bene : je change légèrement cette notation suivant le changement d'argument (« c » pour 
x=Cos(t) ou « s » pour x=Sin(t), et inversion du signe du paramètre q) 


Évidemment lorsque A=0, les développements sont effectivement d'une forme équivalente aux 4 
types de solutions périodiques de l'équation de Mathieu : 


1=+0 


cce, (t;=q)= eer 


q) Cos(21 t) 


e Anl 


1=+0 


q)Cos((21+1)r) 


0 
CO (t, q)= cena lt 


A=0 
e EWDÉ (a +2q Cos(2t)) elt) =0 Si 


et Aral 


(e Leon 


q) Cos(( 21+1)r 


sces, (t, q)= Sen, 


q)= Sin(c (JS At 


1=0 


D 


q) Sin((21+2)r) 


ee BL: 
l=+0 


)Cos(( 


0 
SCE nyi (t, q)= Sesat 


= Sin(t DAS 


q) Sin((21+1)t) 


> Bal 


Parallèlement étudions l'autre forme de Liouville obtenue avec y=Sin(t), à savoir : 


er E (422) 20 07h x?) HUE D) DE x = Sin(t) 


2 
1— x 
2 


ee 
a 4 


2 
A=l-2-u gesot 


wel Jaen? aa 29 Tan (9+(a- À -2qCos(2r)}u, (#)= 0 


u 


y(x)= (Cos (t) u, (t) => u" (t)-21- A) Tan(t)u,'(t)+ (la D — AT -24Cos(2t))u,(t)=0 


On note qu'il suffit de réaliser le changement d'argument t->n/2-t dans les équations précédentes 


u+24-1 


y) = (SiO "` 


y) = (Sin (O) u) = u" (€)+ 21-A) Coran(t}u, (e)+ (a D — AP +24 Cos(2r))u, (1) = 0 


u (£) => Hi dot 2A Cotan(t) u '(t)+ (a — AT +2q Cos(2t))u, (£) =0 


pour obtenir la même chose. On peut donc facilement trouver la récurrence d'une solution de la 


kon 
u(r)= X Bin Sin((21+1)t) 
1=0 


forme Partant du développement : 


E +1) E = )) ? 


1=0 


SA A^ Sin((21+1)r), il suffit d'alterner les signes des coefficients de la matrice. Soit sous forme 
1=0 


matricielle on obtient le système suivant : 


(+AŸ g+(3-A) q 0 0 1 0 0 olf 4 
q g+G+AŸ g+(5-A) q | 1 1 0 À, 
0 q g+(5+AŸ g+(7-AŸ 0 |, 0 1 1 0 |) 4 |, 

0 q g+(21+1-AŸ 0 0 O … oi … 
D 0 q g+(21+1+A) „m 0 0 1 Ara 
0 | 0 0 
[1 -1 1 -1 1 
0O 1 -i 1 1 
e(M-aB)A=0@(B'M-arla=0 pi IT" ea] E SE 
… 0 0 Se: WT ET Ee 
m 0 0 1 
o o o 


Par une linéarisation trigonométrique, on arrive au même résultat. 


De la même manière la récurrence d'une solution de la forme : u(t) 


= SC 4kCos(21 t) solution de 
1=0 


l'équation :“"()—2A Tan(t)u'(t)+ (a-~ —2qCos(2r))u(s) = 0 est construite comme suit : 


g+M 2q+(2-AŸ q 0 0 
q g+(2+AŸ g+(4-A) q SS 
0 q +(4+A) (6-47 0 e 
0 q SEGA 0: 
si 0 q g+(2i+AŸ … 
0 | 
[1 —1 1 -1l 
0 1 -i 1 
e (M-aB)A=05(B"M-a1]ja=0 B=’ | S SS 
ao 0 0 1 
[o o 0 


Donc pour l'équation différentielle intermédiaire —— te). ( a—2q Cos(2t) 


construire 4 types de solutions : 


A(A—1) 


WOP (a —2q Cos(2t) 


Cos’ (t) 


k0-0- 


g()= ene, (t,q)= (Cos 


1 1 0 0 0|| 4, 
O 1 1 O0 … À, 
0 O 1 0 || 4, 
= 0 
D O … O || … 
… … 0 0 1 …||24,, 
|o o Ge es | 
o 
1 
b,,=0 i>j 
dés = |b, | SÉ H 
1 DI P. =(-1}) j 
sal 
A(A=1) g(r)= 0; on peut 
Cos’ (t) 


gle) = cseh, (tq) = (Cos D} E 4) cosl t) 
gle) = wel, Aral (Cos (0): EE 
gle) mei, Jr alzclieerill? OÙ 


B5 —A 


21+1 


(q) Sin((21 +1) 


D)? D A5" (a) Cos(21)1) 


Là encore lorsque À=0, les développements sont effectivement d'une forme équivalente aux 4 types 
de solutions périodiques de l'équation de Mathieu : 


A=0 


a g"'(t)+(a -2q Cos(2t)) g(r)= 0 


> 


=+00 


cse, (t, q) SE? Ce (t, q) be ES (4) Cos(21 t) 


E) 


1=+ 


LUE (t, q) = Sezny (tq) = 


Ce 


=0 


Rei (t, q) = Më? (t, q)= Cos(t (e Aala 


Je 


Gë (tq) Ce (t,q)= Cos(t GES, Alq 


1=0 


BJ Sin((21 + 1) t) 


) Sin((21+1)t)= È sala) Sin((21+2)t) 


l=+%0 


RB. (4) Cos((21 + 1) t) 


1=0 


)Cos(21 t)= 


Quant aux solutions de l'équation AUDE +2q Cos(2r)- AA po- o, la correspondance 


Sin’ (t) 


avec les fonctions de Mathieu est également obtenue lorsque A=0, sous la forme matricielle qui est 
exactement celle des fonctions de Mathieu , notamment pour les ordres impairs : 


1 1 1 1 
O 1 1 1 
a 0 0 1 1 
B M= 
UU … 
E 0 0 
[o 0 
A=0 
l-q -q 
0 = 
= B'M- S 
0 
| 0 
Et pour les ordres pairs : 
[1 1 1 1 
O 1 1 1 
í 0 0 1 1 
B M= 
UU … 
Ve 0 0 
[Oo o0 
A=0 
[0 -2q 
0 = 
= B'M- 1 
0 
| © 


NEA es GA) à 
1|] -4 g+(B+A) -g-(5-Aÿ 
Se 0 -q qg+(5+A) 
1 0 -q 
; | 0 
Jl o0 
0 o] 
0 
q 0 D 

q 0 
-q (+1) … 

e 

Ilan —äe-D-AT q 
(lie g+(@+AŸ -a-(4-AŸ 
SG 0 -q g+(4+A) 
1 0 -q 
. : 0 
J 0 
0 o] 
0 
— q 0 D 

-q 0 
-q DI 


—q-(21+1-A) 
g+(21+1+A) 


0 


-g-(21-A) 0 


g+(21+A) 


0 


Voyons si l'on peut directement construire pour l'équation différentielle intermédiaire : 


u''(t)+ 2A Cotan(t) u'(t)+ (a — AT +2q Cos(2r))u(r) = 0 


des développements possibles à la manière des fonctions de Mathieu, sous la forme 
SE Bh Sin((21+1)r) ou sous la forme 5 BÅ Sin((21)r). Pour t=0, ces deux développements 
1=0 1=0 


s'annulent, ainsi bien évidemment que leurs dérivées secondes. Les dérivées premières sont a priori 
non nulle. Si nous injectons ce résultat dans l'équation différentielle alors : 


u''(t)+ 2A Cotan(t)u'(t)+ (a — AT +2q Cos(2r))u(r). = ZA Cotan(t) NO = 0 


1=0 


Comme A<>0 (sauf pour les fonctions de Mathieu) alors la solution ne peut qu'être identiquement 
nulle. 


De même avec l'équation différentielle #"()-2A Tan(r)u'()+ (a-~ -2q Cos(2r))u(r)= 0 par 
simple GE d'argument, on en déduit qu n'est pas possible de construire des solutions de 
la forme 2 AA „Cos((21+1) £) OU de la forme SS BA Sin((21)t) excepté dans le cas des fonctions de 


1=0 


Mathieu, soit lorsque Ac 


Cas limite y->0 pour u+À entier 


On sait que dans ce cas l'équation des ondes hyper-sphéroïdales se ramène à l'équation 
différentielle des fonctions de Gegenbauer associées. La forme de Liouville permet d'appréhender 
les solutions lorsque u+À est entier : 


D eg 2 (1 24)x ae) Di H 24) de CHA A=l-u-2 x = Cos(t) 


2 
e D 1—x 


preii 1 DEET 


OEM EEE DE 0 210220) > ule) 
S el 29 Cotan(t)u'(t)+ (v 2) -x }u(s)= 0 ce A=l-u-A=1-n 
u"(t)+2(1- u 2) Cotan(t)u'(t)+ (v +1- uv + u+22-1)u(s)=0 a=(v+1) 


On a donc à notre disposition toutes les solutions de cette équation, sous la forme de fonctions de 
Gegenbauer associées de paramètre À+u=n entier et son développement sous une forme finie 
= Sa Cos((v +4 +21-(n-1))t)- Dans notre cas les paramètres v et À sont également reliés par la 


relation v+À est un entier, aussi des développements du type r B2,,Cos((27 +1)r) ou 


T=+oo Le Së e 
u(r)= 2 Bi, Cos(21 t) sont équivalents en posant : 
l =—œ 
2I +1=v+4+2l-n+1 21 =v+A+21-n+1 
Impair 5- v+i-n ou pait ~ v+A-n+l 


l =l+—— l =1+ 
2 2 


Lorsque q=0, La récurrence des coefficients du développement devient : 


A dÉ (7 +1+4) DÉI d a+(A 27 af +q) As a Ais =0 
EENG GA la) a+(A 27 2) +q) Ana =q Ana =0 
Sch (a-(r+1+a)) 44 + a+(A-27 3) Aa =0 
a=(v+4) = 
Katen Lef Allan d a+(A 27 2) ) 44a =0 

(6 +4) -(7+2-n)) Ars (6 +4) -(n + 27 +2)) SE 

(6 +a} -(7+1-n)) A} =(¢ AY (n +21 + PES 


Alors la récurrence des coefficients du développement impair Ae 5. Aż „Cosl2T +1) 1) donne : 
LE 


(wa Biel läis ie Af (+27 +2} ) 45 ia TË EECHER 


Et la récurrence pour un développement pair de la forme „(t)= Ka A? Cos(21 t) donne : 


l=—0 


(Af -QT +1-nf ) 4 =(0 +4} (042741) ) 45 ei 1—/)(v-(n-1)+1+4)4Ż =-(1+1\v +1 +4 +1) 4}, 


Ce sont donc les mêmes récurrences d'avec les coefficients du développement des fonctions de 
B n-1-lv+A-(n-1)+1 
[+1 v+A+l+i 


Gegenbaurer associées de la forme : c:-1(,), à savoir : 


"(t)+2A Cotan(t)u'(t)+ (a MO =0 


D'un point de vue matriciel l'équation “ correspond pour 
> u" (t)+ 2(1-n) Cotan(t)u'(r)+(a-(1=n) Jult)=0 
un développement impair à ce système linéaire : 
KN et es at CP 0 0 8 0 | À 
ob, za? 0 af Aaf 0 e S ER 
2 
0 0 1 I 0 0 (6-n)  … 0 p alear 4 |o 
LE SE 0 0 … —(2+n) 0 | 
de a 0 E A. Saul 2 e 0 O (27+2-n) … Ara 
VE Ae, e Sit, 0 e d i D 8 BE 
IG af -a -12(n-1) -20(n-1) -28(n—1) F IA 
0 (4-n) -a —20(n-1) -28(n-1) s H A 
e 0 0 (6-n) -a -28(n-1) e -A 4 |20 
0 0 T D ll 
Se 0 0 (21+2-nf -a || 4 
L 0 L 


Les valeurs propres de ce système sont 47 Daf ä-alfle-af 


Et pour un développement impair au système linéaire : 


1 1 1 1 — Fifa ts 0 0 0 A 
OS ES GE NEE 0 Daf -(G3+n) 0 G e ER 
0 0 1 1 … …| o 0 (5-nf 0 0 EE 
E eee 2 Ss 0 0 … —(21-1+n) 0 À 
sis 20/08 ZE get 45 S 0 O (21+1-n) … À; 
0 dÉ e 22,041] 0 E e e e 2 st? 
[(i-nŸ a -8(n-1) -16(7-1) -24(n—1) d UA 
0 (3-n) -a -16(n-1) -24(n-1) E DER 
Gi 0 0 (5-nf -a -24(n-1) d OU 
0 0 Si eg 
0 0 (21+1-n) -a …|| Au 


2 
Les valeurs propres de ce système sont 47 (-n),(8-n) (5-7). 


2 2 

A chaque fois les premières valeurs propres (1-1) er al correspond à un développement à un 
2 2 

seul terme, les deuxièmes valeurs propres G-a} et (4-1) à deux termes, les troisièmes valeurs 


propres (S-a? et (6-7? à trois termes et ainsi de suite. Or on a vu dans la méthode de 
construction des fonctions de Gegenbauer associées telles que u+À entier que quelque soit la valeur 
de n>0 : 

- pour v+À=1-n et v+À=2-n, le développement ne contient qu'un seul terme non nul. Pour v+À=1-n 
c'est un développement en cosinus pair et pour v+À=2-n un développement en cosinus impair 

- pour v+À=3-n et v+À=4-n, le développement contient deux termes non nuls. Pour v+À=3-n c'est un 
développement en cosinus pair et pour v+À=4-n un développement en cosinus impair 

- pour v+À=5-n et v+À=6-n, le développement contient trois termes non nuls. Pour v+À=5-n c'est un 
développement en cosinus pair et pour v+À=6-n un développement en cosinus impair 

- et ainsi de suite … 

- pour n=1 et n=0, le développement n'a qu'un terme. 


Il est clair que cette séquence correspond à celle que nous trouvons ici. On peut donc en conclure 
l'identification des valeurs du degré v suivant la succession des valeurs propres établies dans le cas 
des développements pair et impair lorsque y =0 selon la correspondance suivante : 


( v+i=l-n—v+u=1-22- développement en cosinus pair 
a = (2 = nf v+i=2-n—v 
ETA v+A=3-n—v+ 


22 — développement en cosinus impair 


22 — développement en cosinus pair 


u+A=n a = 


v+i=5S-n—v+ 


22 — développement en cosinus pair 


u=2 
OUER 
4- Di v+i=4-n—v+u=4-2À —> développement en cosinus impair 
u=5 
u=6 


6- n) v+i=6-n—v 22 — développement en cosinus impair 


Pour la normalisation des deux types de fonctions de Mathieu associées, 


on peut par exemple étendre la règle de normalisation des fonctions de Mathieu (pour A=0), à 
savoir que la norme du seul développement de Fourier est normalisé à nr : 


f dt (Eata )Cos(21 à) - =r > 24 (4) + 


eier? Cos(2r)- ll glt)=0 Er 


(Ca = 
[af Zar, )Cos((21+1)r Wl State Aao =1 


far (si (Q7 (ecet, (all = far (Si? (Q^ (cceh,.(t,-g)) = far (Sin? (t) IN Jace! / (t- a) = fa (Sin?( M (scet (t Alf =m 
0 0 0 0 
Et de même : 


far Ei A (g)Cos(21 o) SAN 
Side eil -1> Et =: 


Í dt (Cos? OÀ (esef, (t, a) = Í dt (Cos? DIR (ssef (t, a) 


0 0 


© 
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Expression des dérivées premières des fonctions intermédiaires de Mathieu associées 


J'appelle les deux fonctions périodiques comportant les seuls développements en cosinus, les 


fonctions périodiques intermédiaires de Mathieu associées : 
u'"'(t)+ 2A Cotan(t) u'(t)+ (a — AT +2q Cos(2r)}u(r) = 0 
u(t)= ceiz (t, -q)= Za — q)Cos(21 t) 


DOE Sg EE (t, -q)= > Aia = q)Cos((21 + 1) t) 
Les dérivées premières et secondes sont alors triviales : 


1=+00 R | d’ w 
Se AÏ(- q}Sin(21 t) keet, (t.-q)}= 


£ feit (t-q )}= 


Ikon 
— 5 AUT Az (—q)Cos(21 t) 
1=0 


-Fenyan id frein (ra) 


£ let?) G—al- 


Ye E 4%, g)Cos((21+1)r) 


Expression des solutions périodiques en Cos(t) de l'équation des ondes hyper-sphéroïdales avec un 


développement en fonctions de Mathieu intermédiaires associées 


Les fonctions d'onde hyper-sphéroïdales sont donc construites comme suit en 


Al (sr) + DA" (a) Cos(2rn) 
ONCOM Se) coler 
v= (SO Z4 a) Cos((21+ 09 
EN 


Et en argument x : 


+22-1 l=+% 
hal" 
2 ES 
of al (1+24)x nl Luc x) e y(x)=(1-x°) 2 
2 D = 
x Cos(t) A=l-u-2 q= azott +} 
ve Ze 


argument t : 


DETTE 


21+1 


Sa 


L Gë ((21+1) ) ArcCos(x ) 


N 


Y 


7i Ju (21 ArcCos(x )) 


y(x)=(1-x) Ga GE (21 ArcCos(x ) 


2 


A 


a Col (21+1) ) ArcCos(x ) 


Soit avec les fonctions de Mathieu associées intermédiaires : 


+22-1 
y(x) z D -x i 2 ce | recos) 
_u+22- 
stell" eitzi" | äere 


y(x)= D -x E cei" Ange | 


CS 
A 

E 
4 


y(x)= D EEN E e Ange 


Ce qui donne les dérivées premières : 


Et les dérivées secondes : 


NRA ES 


+24 2 +24+1 2 
y'(x)= (x e 2 SIE LC) +(u+22-1)x (1-x° y“ 2 äere GC 
di K 1=ArcCos(x) S 
+22 2 +22+1 2 
y'(x)= Bi -x j 2 — Las , £) +(u+24-1)x (1 -x E SE E HR £) 
d À 1=ArcCos(x) S 
z 2 2 
y'(x)= “(1 -x E SS et? LC) -u d =x’ E) eet: JH Ange £) 
dt à t=ArcCos(x) 4 
| 2 2 
y'(x)= Bn -x E Za? fa gi =H d -x Je? Ee Ange GC 
t=ArcCos(x) 


+14 72 2 +212 2 NIT 2 
il = ( -r e 2 e: W a =y (D u+41- 1) í -x e Ei f- £) + (u +24- 1) (e (u + 24) + d Í -x ES cei" ED a 
t=4AreCos(x) t=4rcCos(x) 
+14 72 H +212 2 +2143 2 
KORI x P 2 - cei" A aus AA? À ceini" ie + (u+24-D{r(u+24)+) LA? ein Arts 
di a t=AreCos(x) di 4 t=4rcCos(x) 
u A 2 -3 . 2 u, 2 
viel x Fësch) al 2ull 2 KE S Wa +u (lu -1)- l-r E ce" reich 
S 4 t=ArcCos(x) di 4 t=ArcCos(x) 4 
2 2 2 2 3 2 
yh) -7 Fäll al 2ull 2 KÉ äi +u DND BE ci dE 
g 4 t=ArcCos(x) di 4 i=4reCos(x) 4 


Expression des solutions périodiques en Cos(t) de l'équation des ondes hyper-sphéroïdales avec un 
développement en polynômes de Chebychev en x 


En introduisant dans les expressions précédentes les polynômes de Chebychev de première espèce, 
on obtient un développement particulier des solutions de l'équation des ondes sphéroïdales. Ils 
sont définis comme suit : 


DATA ern y? 
v(x)=(1- x ) S DA dÉ Glas 


E Est 8 v(x)=(1- Te e Z) T, (x) 
T, (x) = Cos(nt) = Cos(n ArcCos(x)) ru e% Sc. S ) E 


Alto 2 
dich rar Zahl 
1=0 


Avec ces expressions les dérivées premières sont alors les suivantes : 


y'(x) z D = 2 KS SA {- eil {27 T DA 7 x? )U(x)+ X Lu +22 — 1) Taal) 


l=0 4 
yale ER RE MO ere been 


2 


OT Set. | a+), (x à E) 


DO A CT laf 
K 
4 


Et les dérivées secondes : 


u+24+3 kä 1l 


yl- 2 SA LC (u+ 24-1 (u+ 24, Je 3 alt u 24 ui 0, ()+ (211) elt, (x) 


1=0 


4 
CORRECTE Al u+22-1)/(u+24)+1)7,(x)-2 x äh vi {2k-24-44+1)U eh Ah il, ell 


r-t- Ee- AUX ue (u-1)-1)7,0)-2x 11-72 +2k-D 0, el Af li, All 


Def PSE Je te-n)-2)n 6-2 delen tege eeh Al 


Équivalence de l'équation des ondes sphéroïdales avec l'équation des ondes hyper-sphéroïdales 


Un court chapitre qui concerne l'équivalence entre ces deux types d'équation différentielle. On a vu 
que la solution de l'équation hyper-sphéroïdale pouvait se construire avec un développement de la 
forme : 


Ga)" 6)- (1422) stil) HUE D) DEE 


1— x? 
_u+22-1 
y(x)= (i-x°) 2 g(x) = (1 -x J)e") -68 —22- 2u) x g'(x)+ (o+ (1- uY22 +u —1)+ CD -x° )) g(x) =0 
Par le changement de paramètre le développement de la solution l'équation hyper-sphéroïdale 
peut également s'écrire ainsi : 


u>-u+1-27> (1-2) y) -1+ Sale etlech x?) RL) jo 


y(x)= Det — (x )g''(x)-(1+22 +24) x e'(x)+ lo- u (22 + u)+ v°(1-x2)g(x)= 0 


On peur également construire une solution de l'équation des ondes sphéroïdales à l'aide du 
développement suivant : 


Es) -2s 7b) of) ) 9-0 
wll (he) = (1-2 )8"(x)-2(1+ 0) x stellen — ua +1)+ 771 a Detail = 0 
Or la dernière équation associée hyper-sphéroïdale sur g(x) : 
(2x2 )g"(x)-(1+22 +24) x stelle - ulu +24)+ 7° (1-22) gx) = 0 
est équivalente à l'équation associée sphéroïdale : 
(2 lertsl- 203 A) x g'(x) +o- 219 +1)+7 "(x ) 2) = 0 
Par le changement de paramètre : 


(x? )g"(x)-(1+22+2u) x g'(x)+(o-u(u+22)+7°(1-x2))g(x)= 0 
TS E EGF EO) — 21011 
de 2 ` j — D 7 )g"(x)-2(7 +1)» g'(x)+(& aln H 1) H zh xX ))e(&)=0 


Saor SE Ee 


Le développement sur les fonctions sphéroïdales : 
7 u + To 
pst xy? )= l-2} Det) (x) 
Zeen 
conduit à la récurrence à trois termes suivante : 


Aen, bech, -oet )+ Cgf )=0 kef-o,..,-1,..,0,1...4+0} 
(v +2k-1\v+2k) 

SEET 3X2v +2u+4k-—1) 
2(v+24)v+2u+2k+1)+(2u-1)1+w) 
(2v+2u+4k-—1)2v +2u+4k+3) 

2 (V+2u+2k+2)v+2u+2k+1) 
E E ado 


A; = 


=(v+2kv+2u+2k+1)-2y 


C, = 


Et le développement sur les fonctions hyper-sphéroïdales : 
ist (er) (ef Dent) (re) S ei) CA) 


conduit à la récurrence à trois termes suivante pour laquelle est appliquer la substitution des 
paramètres w et u par les paramètres w et y : 


Agina B, -a)na  kef-00,..,-1,..,0,1..,+00} 
Di 5=0+2 -1 2ŭ=2u+2A-—1 


(v +2k-1\v+2k) 
O a Ak- 32v +20 +4k-1) 
2y >(v+2kÏv+2%2+2k+1)+(22-1N1+ 7%) 
(2v +20 +4k-—12v +20 +4k+3) 
(v +2ŭ+2k+2)v+2ŭ+2k+1) 
f (2v+2%+4k+3)2v+2%+4k+5) 


(v+2k—1)(v +24) 
n OFINA IAr Ak- 42v +2u+24+4k-2) 
2 2v+24{v+2u+22+2k)+(2u+22-2)2u+22+1) 
(2v+2u+22+4k-2Y2v+2u+22+4k+2) 
(v+2u+22+2k+1v+2u+22+2k) 
(2v+2u+22+4k+2)2v+2u+22+4k+4) 


A, =y 


=(v+2k)(v +2 +24+1)- 


C, = 


A; = 


4B, =(v+2k\v+2u+24-1+2k+1)-2y 


C, =7° 


Développement en série de la valeur caractéristiques w de l'équation des ondes hyper- 
sphéroïdales allongées lorsque l'hyper-sphéroïdicité y est proche de 0 


D lartsl- e aisect et: SE 2), (x)=0 


1— x° 


La récurrence des coefficients du développement d'une solution de la forme : 


x e [-1.1] hps” #(x, °)= 5C 1) d? (2 eii (x) 
Zeen 


Fonctions associées de Gegenbauer C4, (x) xe [- 1,1] 


—> hps?” (x,0)= dż (0)C2” (x) 


se présente originellement comme suit : 


adit o) (B, -o0 azto? )+C aze G? )=0 bel: Dt) 

y (v—u+2k—1)v-u+2k) 

dE TEEN EES TEST 2) 

2 (V+24k)v+24+22)+u?-1+(22-1Xu+2) 
(2v +4k+2A-—2\2v +4k+24+2) 

a (v +u+2k+2A)\(v+u+2k+2A+1) 

(2v +4k +24 +2)\(2v +4k+24=+4) 


KI 


B, =(v +2k\v +2k+24)-2y 


C; 


Elle peut aussi s'écrire de la façon suivante : 


Adi (°)+ H + S 7 SE E Seed atly dÉ Cd zëck? )- 0 


y? 


CO (v-u+2k-1{v-u+2k) 

T (2v+4k+22-2)2v+4k+22-4) 

à D V +2k)(v +2k+24)+ u’ —1+(22-1Xu+2) 
À (2v+4k+22-2)2v+4k+22+2) 

a __ (v+u+2k+24\v+u+2k+24=+1) 

i (2v +4k+22 +2)2v +4k +24 +4) 


Introduisons un développement de la caractéristique À et des coefficients du développement sous 
la forme : 


j=% ; d?” \y j= 
ol? )= Zur ne Dar 


ya 


Nous connaissons la valeur de la caractéristique et la valeur de la solution pour y=0, il vient alors : 


dł (0 1 si k=0 
œo(0)=1,=v(v+22) O due = LG re 
v,0 F 


Comme par définition pour le coefficient k=0, le rapport est identiquement égal à 1, il vient : 


Zä" =1=>d =l d;=0 j>0 


Écrivons la récurrence pour k=0 : 


o+ B, + Cid 0 +4 HAS d ANA ECH Se SE tel 


j=1 j=1 
Comme d'5 = 40 =0 

SC seg 
5 oi +B -1+5 (4d; +Cd -l,a )y” =0 


J= 


Les termes en puissance de y” s'annulant identiquement, il vient : 


1, =v(v+22) 
1=B 
l'a = 4d; + Cod: ; 


Cela donne une récurrence permettant de calculer les coefficients La dès lors que l'on connaît les 
d d d 


—1,j? 


coefficients Lj, Comment maintenant déterminer les coefficients ad, pour cela 
retournons à la récurrence de départ et insérons les développements dans la récurrence du terme 
d'indice k quelconque : 


j=% 


(Hd + Cydia; + Bid, fe +2kv+2k+24 jS d A Sora = 
j=0 j=0 


j=0 
SES DEE (EE EE AT a)y + ee die SE Zä 
S 


Calculons le premier terme en JA ` We 2 222,0 Sie 0H =v(+24) il 


k#0—d,,=0 
confirme le résultat précédemment établi. 


Calculons le terme en y} il vient : 


(4d,,,+ Cdi, + Bd, ,+((v+2kŸv+2k+22)-1,)d, a)y” SN „7? Si 
j=0 m=0 
m=% i=00 m=% j=% 4 
KI SE Du = y, KSC m Lay? utm) = E dj lu <— j=i+m 
m=0 i=0 m=0 j= i=0 
a (4d,,,+ Cds, + Bd, +(v+2kXv+2k+22)-1)d,,1)7° nu Sais d, Auf 
j=0 j=0 i=0 
> 5 Adya; + Cdi j TE EE E EIERE AE SE y” =0 
j=0 i=0 
> Adra; +Cidyn + Bd, +v +2k\v +2k+24)-h)d, ja SCH ala =0 
i=j i=j 
disk (Ad, tut Bu, ) du An At Eau t Bde ) 
<d = i=0 


SE: (v +2k)v+24+2k)-h E 2k(2k+2v +21) 


=0 


De nouveau il s'agit bien d'une formule récursive, car dès lors que l'on connaît les coefficients 
d'indice : 


Wade diodes EE = Le EC re 4da; EC) 


1,194 j+ 


Ceux d'indice j+1 sont parfaitement définis. L'existence de la récurrence : 


i=j 


d,a- (4d, +C,d 


kri 


+B,d, .) 


k+1,j 


SE 2k(2v+2k+24) 
permet d'établir le résultat important suivant : ***0 di; =0 pour j< W En effet partons de la 
condition initiale "3. . —0 pour k=+1,+2,,+3, qui est vérifiée comme on l'a vu 


précédemment et conforme à l'hypothèse de récurrence. Supposons la propriété vrai pour j : 


| | Si k>0 MS 0 | | 
Vk#0 j<lkl= ell du, =0> SE et d,;=0 te Dei Ae 
k-l,j 
Si k>0 d Adi 
bi (+24 Nv+2k+24)-2 
5 
Si k<0 d, =- Co, 
E (v+2k\v+2k+24)-l 
Si k>0 d,ı;=0 
Pour j+1<|k|= j<|k|-1=> e e wun = 0 
l k+l,j — 


Dans ces conditions le développement du coefficient d. est de la forme (pour k pair) : 


k =+1,+2 — at gé Sa ka y% = of) 
TLTZ,.... at) e Zb 


j=0 


La récurrence s'écrit donc : 


(v-u+2k-1\v-u=+2k) 
(2v +4k+24-2\(2v+4k+24—4) 
_ V +2k\v +2k+24)+ u? -1+ (24 -1\u+4) 
(2v +4k +24 -2)\(2v +4k +24 +2) 
__(v+u+2k+2A\v+u+2k+24+1) 
(2v +4k +24 +2\2v +4k+24+4) 
v(v+22)+u?-1+(22-1{u+21) 
(2v+24-2)2v+24+2) 


k 


,=0 H=v(v+24) Ll=B,=-2 


La = 4da; + Cd; 
dun =1 dy0 =0<—kZ0 di; =0 j>0 


i=j- i=j-1-k| 
d, =0< j<}k|= ON meer y le |k| = Ze Luz Kli 
i=0 
i=j-1- Jr 


EE 


k,j—1-i Lu 
d, LE i=0 
j 2k(2v +2k +24) 


Ici je ne retranscrit que le développement de la valeur caractéristique w jusqu'à l'ordre 8 : 


= 2 (v+4-1v+4+1) 
epi ioni (~v-u+1fv-u+2\v+u+24v+u+24+1) 


EE 8 ut __1f,, (u+4-2{u+2) 
alc Ian L=vv+2) 1, k | 


2 2Ù (v+24-4/2v+22-2) 0v+22)2v-D)  (2v+24)2v+24+2) (2v+22+42v+3) 
HEET (~v-u+1\v-u+2\v+u+24v+u+24+1) : (~v-u-1\v-uv+u+24-2\v+u+24-1) 
Bi 1 (2v +24-2\2v+24\2v +24 +2}(2v +24 +42 +24+6\2v +3} (2v +24-6\2v +22 -42v +24-2} (2v +22)2v +22 +227 1) 


EE 
Ge (-u-1v-uv+u+22-2(v+u+22) F 
(2v +24-6f(2v +24 -42v +22-2) (2v +24(2v +24 +2} (av -17 
(v-u+I(v-u+2v+u+24v+u+24+1) 
(2v +24-2) (2v +24)(2v +24 +2} (2v +24+42v +24 +6% (2v +3) 
e av -u-v -uy +u+24-4)v +u+24-3\v +u +24-2\v +u +24-1) _ 
(2v +24-8)2v +24 -6% (2v +22 -4) (2v +24 -27 (2v +22 27 -3%2v -17 

(~v-u+1\v-u+2\v-u+3\v-u+4v+u+24\v+u+24+1v+u+24+2v+u+24+3) 

(2v +242 +24+2) (2v +24 +47 (2v +22+6) (2v +24+8%2v +3} (2v +5) 


_(v-u+}Ẹv-u+2 +24) +224) (~v-u-1} (vu) +u+24-2)(v+u+24-1) 
(2v +247 (2v +24+2} (2v +24 +47 (2v +3) (2v +24-47 (2v +24-2} (2v +24} (2v -17 
_(v-u-1\v-ufv-u+1\Ļv-u+2\v+u+24-2v+u+24-1(v + u+24v+u+24+1) 


ECHT EIER (2v +24} (2v +24+2} (2v +24 +4\2v -1F (2v +37 


Développement en série de la valeur caractéristiques w de l'équation des ondes hyper- 
sphéroïdales aplaties lorsque l'hyper-sphéroïdicité y est proche de 0 


Pour l'équation des ondes sphéroïdale aplaties : 


Ga)" (9 x stall ch e) H, jo 


1— x 


xe [- 1,1] hps?” (x, y? )= 5 (— 1} 4 (72 ste (x) 


une solution de la forme: rn aura un 
développement en série autour de la valeur y=0 régit par la récurrence suivante : 


k = +1,+2,.. E St d AT ar} Stier 


Led 
(v-u+2k-1\v-u+2k) 
| (2v+4k+24-2\2v +4k+24-4) 
D) V +2k\v +2k+24)+ u? —1+(22-1{u+2) 
(2v+4k+22-2)2v +4k+22+2) 
_ (v+u+2k+24{v+u+2k+22+1) 
| (2v +4k+24+2\2v+4k+24+4) 
v(v+22)+u?-1+(22-1{u+21) 
(2v+24-2){2v+24+2) 


1,=0 I=v(v+24) L=8,-2 


La = 4da; Fa: 
duu zl d,o =0<—kz0 di; =0 j>0 
i=j-1-|k| 
di; =0< j<k=d,,,,=0<i>j-1- H= an = ka -liu 
i=0 
i=j-1-|k| 


d jail tte +B,d, j Cas) 


k,j-1-i"i+1 
d, e i=0 
S 2k(2v +2k+21) 


Introduction à l'étude du comportement asymptotique des valeurs caractéristiques et des 
solutions de l'équation des ondes hyper-sphéroïdales lorsque y très grand 


Venons en au comportement asymptotique des valeurs caractéristiques et des solutions de 
l'équation des ondes hyper-sphéroïdales lorsque y très grand. 


pr) ne El lech x?) PUD) 0 


y(x)= D Bail > D 7 a" (x)-(1+2u +24) x u'(x)+ Le — u (u+22)+ 7°(1—x2))u(x)= 0 


Avec le changement de variable, l'équation différentielle de la fonction u s'écrit ainsi : 


2 2 2 
st Sit E e ODER 
Y 


27 2y 2y 4 
ikun A 
2 SE 2 2 
>f: à jo 2 TO ‘ Juin 
2y y 2y 4 


Dans l'hypothèse où y tend vers l'infini et dès lors que la limite suivante existe sous la forme : 


omet hpl DOME SCOR 
zs 2y 2 2 4 
L'équation différentielle se transforme en équation de Weber dont les solutions sont les fonctions 
paraboliques cylindriques, où | s'avère être le nombre de zéros de la fonction parabolique 
cylindrique sur l'intervalle (-,+ce), lorsque I à une valeur entière. Tout comme pour les fonctions 
sphéroïdales, l'identification des solutions se réalisent avec les propriétés de la fonction hyper- 
sphéroïdale sur ]-1,1[ de parité et de la permanence du nombre de zéros lorsque la 
« sphéroïdicité » y varie de O à l'infini. Lorsque v et u ne sont pas entiers les fonctions hyper- 
sphéroïdales ne sont ni paires ni impaires quand bien même v-u est entier. C'est uniquement dans 
le cas où v et u sont tous deux entiers que la fonction hyper-sphéroïdale de première espèce 
présente la parité suivante: 


vue N hps!(-x)= (1) aps! (x) = (C1) laps! (x) 
Cette parité est également présente pour la fonction parabolique cylindrique ls 27) lorsque | 


est un entier. Pour déterminer la valeur de l, il suffit donc de comparer le nombre de zéros 
respectifs des deux fonctions. Alors I doit également être le nombre de zéros de la fonction associée 
de Gegenbauer c> (x) dans l'intervalle ]-1,1[ lorsque y=0. Il se trouve que le nombre de zéros 


dans ce cas v, u entiers est exactement V4, il vient donc =Y. Comme on a supposé : 


+7 


Lim I+ 1 v—u+ l 
r> 2y 2 2 alors en première approximation on obtient le comportement 
asymptotique de la valeur caractéristique À : 


o@o+y? 
2y 


be nc mor +y Dir che 


On a donc comme correspondance de la fonction hyper-sphéroïdale de première espèce avec la 
fonction parabolique cylindrique la relation suivante à la constante multiplicative Co près : 


Lorsque 7 — hps (x) = Co D — x? je D. (x 2y ) v,ueN 


Partons maintenant de l'équation des ondes hyper-sphéroïdales aplaties : 


(1 vd 2 er A o sch) geed DEE 


y(x) = (1- x? J2 e70u(x) => (1 -x° Ju” (x)+ DG y D -x?°)-(1 +22 + 2u)x}u'(x)+ Lon ulu + 24) y x(1 +22 + 2u))u(x)= 0 


y(x)= (x het (x) => (1 -x° Ju" (x)+ (-2 y D -x?)-(1 +24 + 2u)x)u'(x)+ Lon — ulu +24)+y x(+22+2u))u(x)=0 


Par les deux changements de variable, on obtient l'équation différentielle suivante : 


el eren 2407 leie Due le PRIE (27 D Ju =0 


Lorsque y->c, il vient : 


2u+2Ââ+l 


Lim 4yt vil Jr O4 (ou + 22)- Qu 2440 7)u()=0 


2 
E recu dt onda D 
e 2 4y 4 
Lim 2 g ZZ) 
Et en supposant que la limite suivante existe: "2 + , il vient l'équation 


différentielle , ONE 42421 vil u'(t)+1u(r)=0- Cette dernière équation est celle des 


fonctions de Laguerre où le paramètre | désigne cette fois le nombre de zéros de la fonction de 


24-1 

A 

2 
Laguerre `“ Om l'intervalle [0,+2°). | reste également le nombre de zéros des 
E) 34-i 
7 2 


combinaisons linéaires paire et impaire ra ) Gr Da reel) ee 


L 


sur l'intervalle ]-1,1[ (pour y tendant vers l'infini). Or la correspondance entre fonctions hyper- 
sphéroïdales et fonctions asymptote se présente comme suit : 


t=2y a = x) t t t 2 ori 
~ Dal T7 «fi | Hi E (2 D) dt E 
Pour la même raison que le nombre de zéros de la fonction hyper-sphéroïdale reste identique 


quelque soit la variation du paramètre y, I est alors le nombre de zéros positifs de la fonction 
associée de Gegenbauer sur l'intervalle [0,1], soit : 


PTE si v—ųu pair 
J= 2 
yen 

e e si V—U impair 


En utilisant la variable t il vient ` — 21 +1 =? SEET es 

si V—uU impair 2 
Pour l'identification de la solution, on utilise également la propriété de parité de la fonction hyper- 
sphéroïdale de première espèce uniquement valable lorsque v, u entiers en construisant deux 
combinaisons linéaires paire et impaire des fonctions de Laguerre. La correspondance s'établit 
donc ainsi à une constante multiplicative Co près : 


EES EES 
2 2 


Lorsque y zm hate Ga fer Jr a-a) CD rer" Toy T 


Développement en série de fonctions paraboliques cylindriques des fonctions d'ondes hyper- 
sphéroïdales allongées 


Les résultats précédents sur la forme limite des solutions suggèrent un développement en fonctions 
paraboliques cylindriques de la forme : 


2 u Lë t? 3 k=+00 
v(x)= D JE > CDS; (x 7) (1-5) ur) u(s)= > CD: 
Kk PEES 


=—00 


En injectant ce développement dans l'équation des ondes sphéroïdales transformées en variable 


22 —1 
+. 2 1+ u + — 7 2 2 
fe Aal iii soit : £ t Jo 2 20H ES u(u+2A)+7? t Juin Cette 
2y D 2y 4 


équation différentielle est équivalente à l'équation différentielle suivante, par la substitution de 
deux paramètres : 


2 ~ ~ tf GI 2 
HR. ee ER Jeo LE eue SEE LT 4 Jain 
2 4 2y D 2y 4 


Il s'agit de la même équation différentielle que celle obtenue lors du développement asymptotique 
des fonctions sphéroïdales allongées. Il est donc facile de déterminer une relation de récurrence 
entre les coefficients C.. Pour cela on utilise les propriétés des fonctions paraboliques cylindriques 
aux plus proches voisins : 


SE , Jp 


dt? 2 4 
t D,(t)= Dir o Di 
D. ln D,(t)+p D. ls 011 Dis D,(t)+(p+1)D, (6) 
tD,.()=D,(6)+(p-1)D,,(0) 
dD, (t) _ 
dt 
CD) E ()+pot D, Al = D,:(t)+(20 +1) Din p(p-1) D) 
PD,()=D,40)+(2p+5)D,,4)+(0+2X0 +1) D, (0) 
PD,()=D,6)+(20-3)D,,0)+(0-2Xp-3)D,.,0) 
Le ()=2D,,(0)+(20+10)7D,(6)+p(0-0rD,,() 
= D,(0)+2(20+3)D,.,()+2(29°+2p+1)D,()+2 o (e -1X2 -1) D,-)+ e (0 -1X -2X -3) D,4(t) 


Loi nu 


(ll Die lt lp Zeie 5. Ar 


DN 2p +1 
#D,"()=7 2,0) 8 FD, ()= 


De Dr E p(e Un, Al Zeie 1Xo-2Xo-3)D, Ji 
Datt -22*1p E Did _D,:()-(CP+92,(0)+p(p-12,:() 
P 2 p 7 


Le terme de coefficient C, lorsque l'on injecte toutes ces expressions dans l'équation différentielle 
donne donc : 


2 A ~ — ox 2 2 
D,a "(0 -é D H ën (42 SIE? Dal) CDt) o 


2y y 2y 4 
d z 
-g7 Dach a Dalle 
7 27 
2(1+2k)+1 4+4 +3-40 +2(22-1{1-2u)+2(1+2k) +2(1+2k) 
t 2 t 8y Ob 


E TE IEN UD. EE 1X2+24-2)1+2k-3) D, ,()=0 
y 


X D'alt)-4ù Daalt- 5 +47° +3-40 +2(1+2k) +2(1+2k)-47(2(1+2k)+1)) D(t) 
+48 (1+2k)(1+2k-1) D, ,(6)+(1+2k)(1+2k—1)1+ 24 —2)1 +24 Ge 4()=0 


Si maintenant Ton considère la contrainte de l'annulation de chaque terme de la fonction p, (t) du 


développement, on obtient une récurrence à 5 termes : 


CRC, -40y +342 + (12H) +2(1+24)-47 012) 0) + ; 
4D (1+2k+2)(1+2k+1)C, a +(1+2k+4)(1+2k+3)1+2k+2(1+2k+1)C,., 
Pour k=+1,k=+2,k= +3... 


La première conséquence de l'existence de cette récurrence est que l'on peut appliquer une 
méthode matricielle pour déterminer, lorsque y est suffisamment grand, et que l'on prend 
suffisamment de termes du développement , une valeur approchée de la valeur caractéristique À 
et une valeur approchée de la fonction d'onde sphéroïdale de première espèce par exemple. Mais 
pour l'instant je laisse ce point de côté, essentiellement numérique 


Ici, je me concentre tout d'abord sur les premiers termes de développement en puissance de 1/y. 
Pour cela supposons que la valeur caractéristique À et les coefficients C, Se développement comme 


suit : 
je ` 2 f 
Ser? Co = coo > Co; = 0 j>0 
sachant que eg , 
je 4 Ee E k#0=0c,,=0 Eet 
E 


C, = De. y” 
0 


Et injectons ces développements dans la récurrence à 5 termes pour d'abord déterminer les deux 
premiers termes en y’ et y du développement de À : 


Zon "A0 Zo fär y +4’ Aale 2k)+1)+3-40° +2(1+ 20 +2(1+2k) Se + 
j=-2 j=0 
=0 
TE EIER H(1+2k+4)(1+2k+3\1+2k+2\(1+2k+1 Doae 


j=0 j=0 


Terme en AT, +) sat, 
Co =0 k0 
k=0=T,=2/+1 


Terme en y {47 T,-4y(20+2k)+1))c,,=0= dë 


Voyons maintenant le terme en 1/y pour j=0,1,....(attention j'inverse le signe de membre gauche) : 


j= +0 (ESA j= =+0 j= 400 
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Connaissant tous les coefficients Cnm au niveau m<=l et les coefficients F, au niveau pel on peut 
alors déterminer les coefficients Canı. Pour la valeur spécifique k=0, il vient : 
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Termeen e EU eh dl an kën Dan Xo y On C2; = 0 


S À Coo et 74 Ci; Gu Co; + Lo Ciy + Oo 2j 
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eT = te a se 


Con Con Don Con Don 


Cela permet donc de compléter le calcul des coefficients F, à partir de la connaissance des 
coefficients c,.. Si l'on normalise tous les coefficients c par le rapport : Gd Coo , alors on arrive aux 
deux relations de récurrence combinées : 


=2(1+2k) +2(1+24)+3-40° 
=2/+21+3-4ÿ 
-a =4k(21+2k+1) 
B,=16k 
y =4à(1+2k+2)(1+2k +1) 
8, =(1+2k+4)(1+2k+3\1+2k+2\1+2k+1) 
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J'ai vérifié cet algorithme de récurrence d'après le programme Mathematica élaboré par P.Falloon 
et j'ai retrouvé exactement les valeurs données par NIST Handbook of Mathematical Functions, 
30.9 Asymptotic Approximations and Expansions, Formules 30.9.1, 30.9.2 et 30.9.3, mais avec les 


nouveaux paramètres Hm: 


I=v-u q=2+1 H=u+À ö=0+2-1 
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2”T, =—527 q -61529 q° -1043961 q° -2241599 q + 32 q (298951127550 q? + 5739 q“) 
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er Et 


2 _10240 q (301+71q°) à +65536 q jl° 


Je vais maintenant démontrer deux propriétés supplémentaires importantes de ce développement : 


is 
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Deux propriétés supplémentaires e 
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) k = +1,+2,.. 


j>1 


=0 tel que ll 


La première propriété à démontrer est celle-ci: Tri 72}, La propriété est vraie pour 


j=0. Pour j=1, il vient : 


x l ~ ~ 
k=+1=c,, = B (4u Ck-1,0 7 Xk čao) 0 
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= (NES 8 e à. SG tre d lys Ge 3 
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Si elle est vrai pour j l'est-elle pour j+1 : 


C; =0 tel que k| >2j 
i=j-1 


2 1 Gi Si _ z D x 
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La première propriété est démontrée. 


=0 Tel oe l+2k<0 


La deuxième propriété est la suivante : "ru . La propriété est vrai pour j=1, en 


effet : 

B=164=B,=-16 B,=-32 
Čo =l sy, =4(+2k+2)(1+2k+1)= y, = 41 (1-1) 
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Si la propriété est vrai pour j, l'est-elle pour j+1 : 
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La deuxième propriété est donc démontrée. 


Cette dernière propriété est fondamentale afin de prouver que le développement ne comporte 
jamais de fonctions paraboliques cylindriques d'indice négatif. Le développement a donc toujours 
la forme suivante : 


k=+0 | 
vele Co SCD 27 ) 


REA YH 
2 


Par le développement en puissance inverse de y ce résultat est démontré, mais une considération 
simple permet également de suggérer un tel développement. En effet la condition de finitude à 
l'infini doit s'appliquer pour la fonction D,a (x 2y ) lorsque y tend vers l'infini. Or ne serait-ce que 


pour p, (x 2y ) cette condition n'est pas remplie puisque : 


Dei, JS. Hi er nein dg, Bes 


Or toutes les fonctions paraboliques d'indice négatif se comporte de la même manière. 


Développement en série de fonctions de Laguerre des fonctions d'ondes hyper-sphéroïdales 
aplaties 


Les résultats précédents sur la forme limite des solutions suggèrent un développement en fonctions 
de Laguerre de la forme : 


wels D fer leen ape nerok rt es) (ea )+ (1) inte 


En injectant ce développement dans l'équation des ondes sphéroïdales transformées en variable : 


LAURE Lal die date ancien: 1) (usa Lal OD 
lara? Lu JO) ami d (urrin) eOr tert 1) [usa Lal OD 


Posons H=u+2 : BET A REI 
= 47 (tu (e)+ (E +1- e)l) -Pu l) EI ettlecles Aalen Aler dl 


S 4y (eu (8 A1 Awe) verk 2a vele gl AR H)-(G ee -))ul)=0 
Saywa- (0) Pa) 2 (7 +0) Du (0 dë 8 a 18 +1) (27 else) 0 


Il se trouve que cette équation différentielle est exactement celle obtenue lors de l'étude du 
développement asymptotique des fonctions d'onde sphéroïdales aplatie, avec la substitution de 
paramètres: H>H @—®, On peut donc facilement déterminer une relation de récurrence 
entre les coefficients C.. Pour cela on utilise les propriétés des fonctions paraboliques cylindriques : 


t IP" + (a +1- OLR + (+4) (0) = 0 
t LE) GU k) 1+ T) LE- U k+ DR a) k aCe) 
t ELO) = UAE E) -0 k + EL ©) 
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t Zë 
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Ainsi, il vient pour le terme C, de la série en utilisant les paramètres #-® 


ay (IR OE- AR O)-r Gët). GG 01. le Hëlt (a — al a)-e H 0-0 
a, =((+k+nX1+k) 

Posons B, = —-2 (2(1+k)+ +1) > ara L aE) Ber + A LE ()+a,1f) Als d 
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Puisque le terme de la fonction de Laguerre lte) doit s'annuler, il vient la relation de récurrence à 
C,=0 k<- 
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trois termes : w()= 0S HESE -(©+ Bıy + Xa)Cı =0 


CILONI -(©+ B, Y + Xr) Cı +a,C,=0 


En introduisant le paramètre q, il vient : 
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Tout comme pour les fonctions d'ondes sphéroïdales allongées, la première conséquence de 
l'existence de cette récurrence est que l'on peut appliquer une méthode matricielle pour 
déterminer, lorsque y est suffisamment grand, et que l'on prend suffisamment de termes du 
développement, une valeur approchée de la valeur caractéristique À et une valeur approchée de la 
fonction d'onde sphéroïdale de première espèce par exemple. 


Mais là encore pour l'instant je laisse ce point de côté, essentiellement numérique, afin de 
déterminer tout d'abord les premiers termes de développement en puissance de 1/y. Pour cela 
supposons que la valeur caractéristique À et les coefficients C, se développement comme suit : 
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Et injectons ces développements dans la récurrence à 3 termes pour d'abord déterminer les deux 
premiers termes en y’ et y du développement de À : 
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Voyons maintenant le terme en 1/y pour j=0,1,.... 
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Con Con Con 


Cela représente donc un système de deux relations de récurrence combinées, et si l'on normalise 
tous les coefficients c par le rapport : Cam/ Coo , alors il vient : 


(t+2k-1+Ar+2k-1- ñ) 
a, = 


4 a, ={1+k+ü1+k) 
B,=-2(r+2k) © p, =-2(2(1+k)+ +1) 
= Et à +1)(r+24-1 DER a = +k+ 4) +) +08 y-k) 
kT T T 
2 
Gels et 1>0 Seat Sch 
GE) i=j-1 
CAE H Du Of H Dr GË Ou H (x ACT CARRE H At GE 
Z nm a= Dm = = At il k>0 
„m Gii 
Ef Fuite [= Eë) GER Q1(1++1)+1+à) 
T = ta; (el 


J'ai vérifié cet algorithme de récurrence d'après le pro 


gramme Mathematica élaboré par P.Falloon 


et j'ai retrouvé exactement les valeurs données par NIST Handbook of Mathematical Functions, 
30.9 Asymptotic Approximations and Expansions, Formules 30.9.4, 30.9.5 et 30.9.6, mais avec les 


nouveaux paramètres H: 
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